Modeling and Simulation of an Air-Assist Atomizer for Food Sprays

Size: px
Start display at page:

Download "Modeling and Simulation of an Air-Assist Atomizer for Food Sprays"

Transcription

1 ILASS Americas, 21 st Annual Conference on Liquid Atomization and Spray Systems, Orlando FL, May Modeling and Simulation of an Air-Assist Atomizer for Food Sprays Franz X. Tanner 1, Kathleen A. Feigl 1, Tim O. Althaus 2 and Erich J. Windhab 2 1 Department of Mathematical Sciences, Michigan Technological University, Houghton, MI , U.S.A. 2 Laboratory of Food Process Engineering, Swiss Federal Institute of Technology (ETH), CH-892 Zurich, Switzerland ABSTRACT An air-assist atomization model for food sprays has been developed and validated. The model is based on the Cascade Atomization and Drop Breakup (CAB) model. The modifications include a change in the product drop distributions, namely, instead of a uniform distribution, as used in the original CAB model, a χ-squared distribution with the same average drop size has been assumed. The second modification addresses the air-assist atomization process. This process has been modeled by estimating the Weber number due to the increased relative velocity caused by the air jet. This leads to exactly one breakup depending on whether the Weber number is in the catastrophic, stripping or bag breakup regime. The model changes have been validated with a low-pressure atomizer for a cocoa butter melt spray, and for an air-assist atomizer using a highly viscous nutrious liquid. The simulations were performed with a modified version of the KIVA-3 CFD code and they showed good agreement with the experimental data. Corresponding author: phone: (96) , fax: (96) , tanner@mtu.edu

2 Average Diameter [µm] p inj =6 bar p inj =12 bar Experiment Average Diameter [µm] d noz =267 µm d noz =4 µm d noz =6 µm Time [ms] Figure 1. Average drop diameter of the entire spray for different injection pressures Time [ms] Figure 2. Average drop diameter of the entire spray for different nozzle orifice diameters. INTRODUCTION Sprays in food processing are used to produce a powdered substance via spray-drying or spray-freezing. Depending on the application, this powder is required to have specific macroscopic properties such as a desired drop size distribution, and microscopic structures which enclose exact amounts of different substances in each of the powder particles. The liquids involved in food sprays can be very complex, typically consisting of multi-phase/multi-component suspensions or emulsions, which have anisotropic structure and can exhibit non-newtonian and even viscoelastic behavior. Therefore, high-pressure atomizers are not suited for food sprays because the high pressures needed in this process could alter the complex material structures and thus change taste, shelf life and other essential food properties. Using pressure atomizers at relatively low injection pressures is not an option either, as a computational investigation of a highly viscous nutriose liquid spray revealed. In fact, as is seen in Fig. 1, an injection pressure increase from 6 bar to 12 bar had almost no influence on the average drop sizes. The reason for this behavior lies in the fact that the drop Weber numbers 1 were mainly subcritical, i.e., for p inj = 6 bar, We 5.5 < We crit = 12 and therefore, there is only insignificant breakup activity, if there is breakup at all. In addition, the average drop sizes are considerably larger than the ones obtained from the experiment. Moreover, the drop sizes are determined by the nozzle diameters, as is illustrated in the nozzle orifice variation simulations in Fig The Weber number is defined as We = ρ g Dv 2 r /σ, where ρ g is the gas density, D the drop diameter, v r the relative gas-drop velocity, and σ is the surface tension. These findings support the fact that in food-related sprays air-assist atomizers are mainly used. In airassist atomizers, the liquid exiting the nozzle interacts with air jets which leads to the liquid breakup. The interaction between the air and the liquid jets is a complex multi-scale problem which cannot be resolved for a realistic spray by present-day computers. Therefore, the simulation of an air-assist atomizer poses a modeling challenge which is the subject of this study. The starting point of this modeling approach is the Cascade Atomization and Drop Breakup (CAB) model developed in previous studies [1] to describe the atomization process of pressure atomizers. In the CAB model, large droplets of orifice size are injected which subsequently break up into tiny droplets via a cascade of drop decays until the droplets reach a stable condition. Each individual breakup event is caused by aerodynamic instabilities and its nature depends on the drop Weber number. More precisely, the breakups are modeled after the experimentally observed bag, stripping, and catastrophic breakup regimes. In an air-assist nozzle, the interaction of the air streams with the liquid jets increases the relative velocity between the gas and the liquid. This increases the Weber number and hence determines the liquid breakup. This liquid-gas interaction occurs only within a short distance near the nozzle exit, and therefore, the air-assisted breakup occurs exactly once. In the air-assist breakup model introduced in this study, this first breakup is modeled using the breakup criteria of the CAB model, taking the increased Weber number into account. After the initial breakup, the droplets are subject to the usual air-droplet interaction described by the standard CAB model. In addition to the air-assist modeling, since the pri-

3 mary breakup occurs only once, the uniform product droplet assumption results in unrealistic product drop size distributions. Therefore, the product droplets are equipped with a χ-squared distribution which has the same mean value as the uniform distribution. This air-assist breakup model has been validated with experimental drop size distributions obtained from sprays of a cocoa butter melt and a highly viscous nutriose liquid. Both substances were injected into a freezing environment at ambient pressure which resulted in a powder. The corresponding simulations have been performed with a modified version of a Kiva-3-based CFD code. It should be noted that in the simulation results presented in this study, the phase change due to freezing has been neglected. This phenomenon will be considered in a future investigation. BREAKUP MODELING The modeling of the air-assist atomization process is based on the Cascade Atomization and Drop Breakup (CAB) model [1]. A summary of the CAB model is given below, followed by a detailed description of the low-pressure atomization and the air-assist atomization modeling. The CAB Drop Breakup Model The basic idea of the CAB model is the simulation of individual breakup events which are modeled after the experimentally observed bag, stripping or catastrophic breakup mechanism, as reported by Liu and Reitz [3]. The actual breakup criterion of a single drop is determined from the Taylor drop oscillator [4] as introduced by O Rourke and Amsden [5] in the context of spray simulations. In this approach, the drop distortion is described by a forced, damped, harmonic oscillator where the forcing term is given by the aerodynamic droplet-gas interaction, the damping is due to the liquid viscosity and the restoring force is supplied by the surface tension. Breakup occurs when the normalized drop distortion, y(t), exceeds the critical value of one. The behavior of the product droplets is derived from a drop creation rate equation, which, in conjunction with mass conservation and a uniform product drop size distribution assumption, leads to the relation r a = e K but bu, (1) where a and r are the radii of the parent and product drops, respectively, t bu is the breakup time, and K bu is the breakup frequency which depends on the breakup regime. These breakup regimes are classified with respect to increasing gas Weber num- Table 1. Constants used in the CAB model. C λ jet breakup length coefficient 5.5 θ spray angle [deg] 8 & Ref. [2] k 1 breakup regime constant.5 n exponent in IDSD.5 bers into bag breakup (We crit < We We b,s ), stripping breakup (We b,s < We We s,c ) and catastrophic breakup (We > We s,c ), where the regime-dividing Weber numbers are taken to be We crit = 12, We b,s = 8 and We s,c = 35, as suggested in Liu and Reitz [3]. More formally, the breakup frequency can be expressed as k 1 ω if We crit < We We b,s K bu = k 2 ω We if We b,s < We We s,c k 3 ωwe 3/4 if We s,c < We where, according to [5], the drop oscillation frequency, ω, is given by ω 2 = 8σ ρ l a 3 25µ2 l 4ρ 2 l a4. In this equation, σ denotes the surface tension, ρ l the liquid density, a the drop radius and µ l the liquid viscosity. As a consequence of continuity considerations of K bu at the regime-dividing Weber numbers, only one model constant, k 1, is required to characterize the different breakup regimes; its value has been determined to be k 1 =.5 (cf. [1]). An additional property of the product droplets are their initial velocities. The axial velocity is inherited from the parent drop, whereas the transversal (radial) velocity, which is responsible for the radial expansion of the spray, is derived from an energy conservation argument involving the surface, kinetic and deformation energies of the parent and product drops. Jet Breakup Modeling The simulation of a fragmented liquid core at the nozzle exit is achieved by injecting large drops of the size of the nozzle orifice. These drops eventually break up into smaller product droplets until they reach a stable condition, thus forming a breakup cascade where each breakup event is governed by the cascade breakup law given in Eq. (1). In order to avoid the almost immediate breakup of the highly unstable initial drops, each (2)

4 Table 2. Data used in the model validations for the cocoa butter (CB) melt and the nutriose liquid (NL). Parameter CB NL Atomization type pressure air-assist Gas type air air Gas temperature [K] Gas pressure [bar] 1 1 Liquid temperature [K] Liquid density [g/cm 3 ] Orifice diameter [mm].5.4 Injection duration [ms] 1 1 Injection pressure [bar] 6 6 Injected fuel mass [mg] drop is equipped with a drop deformation velocity to prolong its lifetime such that the jet breakup length matches experimentally determined correlations. In order to account for the droplet surface stripping near the nozzle exit, the initially injected drops have been equipped with an initial drop size distribution (IDSD) such that the small droplets reflect the surface stripping and the large drops yield good penetration and simulate the fragmented liquid core. For highpressure sprays, this (integral) drop size distribution is formally given by the power law H(d) = ( d do ) n+1 if < d < d otherwise where d and d o are the drop and nozzle radii, respectively. The model tuning described in [1] has resulted in a value of n =.5. The initial value of the normalized rate of drop deformation, ẏ o = ẏ(), is chosen such that the first (and only the first) droplet breakup time is considerably delayed, while keeping the initial drop deformation y o = y() =. With this choice of initial conditions, the injected drops have an extended life span; in fact, ẏ o is determined numerically as the largest negative root of the equation y(t bu ;, ẏ o ) = 1, (3) where y(t;, ẏ o ) is the solution of the Taylor drop oscillator with initial conditions (, ẏ o ), and the value 1 is the critical deformation. More specifically, for an inviscid fluid, the largest negative root of Eq. (3) is given by ẏ o /ω o = [1 We 12 (1 cos ω ot bu )]/ sin ω o t bu, (4) where ω 2 o = 8σ/(ρ l a 3 ). The value of the breakup time, t bu, of a high velocity liquid jet injected into a gas is determined from the experimental jet breakup length correlation due to Levich [6] L = v inj t bu = C λ ρl ρ g d, (5) where d is the nozzle diameter and v inj the jet exit velocity. The constant C λ is nozzle dependent and for the computations considered in this study, a value of C λ =5.5 has been used. This value gave good comparisons with the experimental jet breakup lengths. In the CAB model, the spray angle, θ, is prescribed as an initial condition. This approach requires a correction to the radial drop velocity at the first breakup, as described in [1]. The spray angle for the lowpressure cocoa butter sprays were taken to be 8 degrees, and the one for the air-assisted nutriose liquid sprays were determined automatically from the correlation of Naber and Siebers [2]. The CAB model constants are summarized in Table 1. The constants θ and C λ depend on the nozzle and on the injection system specific properties and, in general, need to be adjusted in order to compensate for such influences. Low-Pressure Atomization Modeling In the original CAB model, the product droplet distribution was taken to be uniform, i.e., after a breakup all product droplets have the same diameter determined by the breakup frequency K bu in Eq. (2). The uniform product drop size distribution assumption works well for high-pressure sprays because the initially injected large drops fall into the catastrophic breakup regime and therefore, they undergo several breakups (in different breakup regimes) until they reach a stable condition. Also, it should be noted that the uniform product drop size assumption is realistic for the stripping breakup regime. For low-pressure sprays, however, the initially injected droplets are either stable or they fall into the bag breakup regime, which means that the product droplets are not uniform. Consequently, the assumption of a uniform product drop size distribution leads to unrealistic drop sizes. This is also the case for airassist atomizers, because the air jet causes only one drop breakup at high Weber numbers (usually in the catastrophic regime), whose product droplets are in general not uniform. To describe a more realistic non-uniform product drop behavior, the product droplets are assumed to follow a χ-squared distribution whose mean value corresponds to the product droplet radius determined

5 .1 χ 2 : f(r)=e -r/1 /1 Vol-weighted χ 2 : g(r)= (r/1) 3 e -r/1 /6 5 4 Air jet exit velocity Speed of sound [1/µm].5 Velocity [m/s] Drop Radius [µm] Gas Pressure Ratio [p /p] Figure 3. χ-squared and volume-weighted χ-squared distributions for a hypothetical case with r = 1 µm. via Eq. (1). More precisely, the χ-squared droplet distribution is given by f (r) = exp( r/ r), 1 r where r is the mean drop size determined by Eq. (1). Observe that r = r f (r)dr, i.e., r is also the mean value of the χ-squared distribution. Consequently, Eq. (1) is to be interpreted as determining the mean values of the product droplet distributions. In the actual program implementation, the drop distributions are discretized into parcels, where each parcel is a collection of droplets of equivalent states [7]. In particular, each parcel contains drops of identical mass, and therefore, the χ-squared distribution for the drop radii becomes the volume-weighted density distribution for the parcels g(r) = 1 A r3 f (r). The normalization g(r)dr = 1 yields A = 6 r 3 and, therefore, the χ-squared parcel distribution becomes g(r) = 1 6 ( r r ) 3 1 exp( r/ r). (6) r Eq. (6) is the actual distribution used in the modeling of the product drop sizes. The difference between the χ-squared distribution and the volume-weighted χ-squared distribution is illustrated in Fig. 3 for a hypothetical case with r = 1 µm. Air-Assist Atomization Modeling In the air-assist atomizer used in this study, there are six air jets which intersect the liquid jet transversally Figure 4. Air jet exit velocity as a function of the gas pressure ratio for ambient gas conditions of 1 bar and 273 K. near the nozzle exit. These air jets increase the relative gas-liquid velocity, which increases the Weber number of the liquid, and consequently, leads to the liquid breakup. This air-liquid interaction occurs only on a small region of the liquid jet which determines the primary jet breakup. The air jet exit velocity is obtained from the compressible Bernoulli equation γ p γ 1 ρ + gz + u2 2 = const along streamlines, where γ is the isentropic coefficient, p is the pressure, ρ the density, g the gravitation constant, z the height, and u is the velocity. Note that γ = c p /c v 1.4 for an ideal gas, where c p and c v are the specific heat capacities at constant pressure and constant volume, respectively. Solving this expression for the air exit velocity, assuming an isentropic process (p ρ γ ), neglecting height differences and using the ideal gas law, one obtains 2γ u = γ 1 RT 1 ( p p ) γ 1 γ, (7) where R is the specific gas constant and the subscript indicates the conditions at the nozzle exit. It is interesting to note that for a pressure ratio p /p =.5 and a gas jet temperature of T = K, one obtains an air exit velocity of u = m/s, which is very close to the speed of sound. Since the air jet exit velocity cannot exceed the speed of sound, u sound = γrt, the air jet exit velocity is taken to be u exit = min{u, u sound }, a relationship which is illustrated in Fig. 4 for the ambient gas at 1 bar and

6 Standard mesh: 5'5 cells Fine mesh: 167'97 cells Weber Number Catastrophic breakup Stripping breakup Bag breakup Nozzle Distance [mm] Gas Pressure Ratio [p /p] Figure 5. Liquid jet Weber number as a function of the gas pressure ratio, and the associated breakup regimes. 273 K. The effect of the air jet on the liquid jet Weber number as a function of the gas pressure ratio p /p is illustrated in Fig. 5 for the injection configuration whose data are listed in Table 2. This figure also relates the gas pressure ratio to the respective breakup regimes; for instance, if p /p =.5 then the resulting liquid jet breakup lies in the catastrophic regime. COMPUTATIONAL DETAILS The computations presented in this study have been performed with a modified version of the KIVA-3 code [7] equipped with many new or improved models. As described above, the CAB atomization and drop breakup model [1] has been utilized. The turbulence has been accounted for via the RNG k-ε turbulence model as implemented by Han and Reitz [8]. Unless stated otherwise, all the standard values of the model parameters, as reported in the respective citations, have been used. The liquid was injected axially downward along the centerline of the cylindrical spray tower by means of a single orifice nozzle. The air in the spray tower was cooled by spraying liquid nitrogen into the chamber via several nozzles located near the top of the tower. Simulations were performed for a cocoa butter melt and a nutriose liquid. More details are given in Table 2. The mesh dependence of the air-assist model has been investigated in terms of the spray penetration for nutriose liquid sprays. The results are presented in Fig. 6, which shows that the penetrations are virtually mesh independent, and the choice of the standard mesh is justified. The meshes used are structured, hexahedral, polar meshes whose cells are concentrated Time [ms] Figure 6. Mesh dependence of the spray penetration for the standard and the fine mesh. radially and vertically around the nozzle exit in the top center of the cylindrical geometry, as can be seen in the vertical cut plane in Fig. 9. The cylinder part of the standard mesh has cells (a total of 5 5 cells) in radial, azimuthal and axial directions. The smallest cell is at the nozzle exit measuring approximately 1.5 mm by 1.5 mm in radial and axial direction. Resolution changes in all directions by a factor of 1.5 have resulted in the fine mesh with (167 97) cells. MODEL VALIDATIONS The two modifications of the CAB model, the volume-weighted χ-squared distribution of the product droplets and the air-assisted primary breakup, have been validated separately. The χ-squared distributions have been compared with experimental data of a cocoa butter melt for a pressure atomizer, and the air-assisted primary breakup has been investigated for a nutriose liquid in a second step. The experimental data have been obtained in the Laboratory of Food Process Engineering at the Swiss Federal Institute of Technology (ETH) Zurich. A summary of the validation conditions is given in Table 2, and the model parameters used are listed in Table 1. Low-Pressure Atomization Model Validation The χ-squared drop size distribution introduced into the breakup model has been validated with experimental data of a cocoa butter melt spray for a simple, solid-cone, pressure atomizer with a nominal spray angle of 8 degrees. The computed volume-weighted drop size distributions of the entire sprays are compared with experimental data in Fig. 7 for the liquid

7 q3(x) [1/µm].1.5 Experiment: p inj =4 bar (subcritical) Experiment: p inj =6 bar Experiment: p inj =8 bar Average Diameter [µm] p /p=1 (no air jet) p /p=.91 (bag BU) p /p=.67 (stripping BU) p /p=.5 (catastrophic BU) Experiment Drop Diameter [µm] Figure 7. Volume-weighted χ-squared drop size distributions for various injection pressures. The filled symbols denote the experiments and the open symbols the simulations Time [ms] Figure 8. Average drop diameter of the entire spray for different air jet pressure ratios of an air-assist atomization process of a nutriose liquid. injection pressures of 4 bar, 6 bar and 8 bar. As seen in Table 3 the injection pressure of 4 bar leads to drop Weber numbers of around 12 which is the critical bag breakup limit. Therefore, the breakup activity is very low, which is reflected in the relatively large drop size distribution in Fig. 7. Note that the agreement between simulation and experiment is reasonably good. The two higher injection pressures of 6 bar and 8 bar lead to droplets at the nozzle exit whose Weber numbers are supercritical; in fact, they are in the bag breakup regime. Here, the breakup activity is large which results in smaller droplets. Again, as is seen in Fig. 7, the agreement between the experimental and simulated drop size distributions are again in good agreement. Air-Assist Atomization Model Validation The air-assist modeling has been investigated for different air jet pressure ratios of nutriose liquid sprays (cf. Table 2). The pressure ratios determine the air jet Table 3. Jet exit velocities and Weber numbers for different injection pressures of the cocoa butter melt at 318 K. p inj [bar] v inj [m/s] We exit velocities according to Eq. (7). The pressure ratio p /p = 1 indicates that there is no air-assist atomization, and the ratios of p /p =.91,.67,.5 result in an air jet-induced bag breakup, stripping breakup or catastrophic breakup, respectively. The average drop sizes of the total sprays, as a function of the injection time, are shown in Fig. 8. This figure shows that the increased air jet pressures, which lead to the different breakup regime of the primary breakup, result in smaller drop sizes. In fact, the average drop sizes obtained with the pressure ratio p /p =.5 are in very good agreement with the experimental results. A comparison between a pure pressure atomization (p /p = 1) and an air-assist atomization in the catastrophic regime (p /p =.5) is shown in Fig. 9 in a vertical cut plane at an injection time of 2 ms. These pictures show that the pure pressure atomizer, apart from larger drop sizes, leads to a narrower spray and deeper penetration. For the air-assist atomizer, the drop sizes do not change much after the primary breakup at the nozzle exit. Both atomization cases show increased drop sizes near the tip of the spray. This phenomenon has been called the tip-clustering effect in a previous computational study by Tanner and Boulouchos [9]; it has been explained in terms of tail wind interactions between successively injected drops. Experimental evidence of this phenomenon has been obtain by means of X-ray measurements and is discussed in [1].

8 although, the actual freezing process has not been taken into account. The phase change due to freezing is the subject of a future investigation. REFERENCES [1] F. X. Tanner, Development and Validation of a Cascade Atomization and Drop Breakup Model for High-Velocity Dense Sprays, Atomization and Sprays, vol. 14, no. 3, pp , 24. [2] J. Naber and D. Siebers, Effects of Gas Density and Vaporization on Penetration and Dispersion of Diesel Sprays, SAE Transactions: Journal of Engines, vol. 15, no. 3, [3] A. B. Liu and R. D. Reitz, Mechanisms of Air- Assisted Liquid Atomization, Atomization and Sprays, vol. 3, pp , Figure 9. Pure pressure atomization (left) and airassist atomization (right) of a nutriose liquid spray at the end of the injection cycle. (The colors represent the drop sizes.) SUMMARY AND CONCLUSIONS The CAB model has been modified to accommodate atomization processes for low-pressure, air-assist atomizers. The modifications include a change in the product drop distributions. Instead of a uniform distribution, as used in the original CAB model, a χ- squared distribution with the same average drop size has been utilized. In the actual implementation, a volume-weighted distribution has been used in order to account for the parcel discretization of the spray probability distribution function. The second modification addresses the air-assist atomization process. This process has been modeled by estimating the Weber number due to the increased relative velocity caused by the air jet. More precisely, the air jet exit velocity has been estimated by means of the compressible Bernoulli equation, assuming an isentropic process. The interaction between the liquid and the air jets leads to exactly one breakup whose characteristics depend on whether the Weber number is in the catastrophic, stripping or bag breakup regime. The model changes have been validated with a low-pressure atomizer for a cocoa butter melt spray, and for an air-assist atomizer using a highly viscous nutriose liquid. The simulations were performed with a KIVA-3-based code which is equipped with well-established spray models. The validations showed good agreement with the experimental data, [4] G. I. Taylor, The Shape and Acceleration of a Drop in a High Speed Air Stream, The Scientific Papers of Sir Geoffrey Ingram Taylor (G. K. Batchelor, ed.), vol. 3, pp , Cambridge University Press, [5] P. J. O Rourke and A. A. Amsden, The TAB Method for Numerical Calculation of Spray Droplet Breakup, SAE Paper 87289, [6] V. G. Levich, Physicochemical Hydrodynamics, Prentice Hall, Englewood Cliffs, N.J., pp , [7] A. A. Amsden, P. J. O Rourke, and T. D. Butler, KIVA II: A Computer Program for Chemically Reactive Flows with Sprays, Tech. Rep. LA MS, Los Alamos National Laboratory, May [8] Z. Y. Han and R. D. Reitz, Turbulence Modeling of Internal Combustion Engines Using RNG k-ε Models, Combust. Sci. and Tech., vol. 16, pp , [9] F. X. Tanner and K. Boulouchos, A Computational Investigation of the Spray-Induced Flow and Its Influence on the Fuel Distribution for Continuous and Intermittent DI-Diesel Sprays, SAE Paper 96631, [1] F. Tanner, K. Feigl, S. Ciatti, C. Powell, S.-K. Cheong, J. Liu, and J. Wang, Structure of High- Velocity Dense Sprays in the Near-Nozzle Region, Atomization and Sprays, vol. 16, pp , 26.

5. SPRAY/WALL IMPINGEMENT

5. SPRAY/WALL IMPINGEMENT 5. SPRAY/WALL IMPINGEMENT 5.1 Wall Interaction Regimes Wachters and Westerling (1966), Akao et al. (1980), Senda et al. (1994) and Nagaoka et al. (1994) describe in detail the phenomena observed when drops

More information

Chapter 2 Spray Formation and Penetration

Chapter 2 Spray Formation and Penetration Chapter 2 Spray Formation and Penetration Abstract The conventional understanding of spray formation when liquid leaves the nozzle is based on the analysis of the following stages: development of a jet,

More information

AN INVESTIGATION OF THE MESH DEPENDENCE OF THE STOCHASTIC DISCRETE DROPLET MODEL APPLIED TO DENSE LIQUID SPRAYS. Simone E.

AN INVESTIGATION OF THE MESH DEPENDENCE OF THE STOCHASTIC DISCRETE DROPLET MODEL APPLIED TO DENSE LIQUID SPRAYS. Simone E. AN INVESTIGATION OF THE MESH DEPENDENCE OF THE STOCHASTIC DISCRETE DROPLET MODEL APPLIED TO DENSE LIQUID SPRAYS By Simone E. Hieber A THESIS Submitted in partial fulfillment of the requirements for the

More information

TRAJECTORY BASED DROPLET COLLISION MODEL FOR SPRAY MODELING

TRAJECTORY BASED DROPLET COLLISION MODEL FOR SPRAY MODELING TRAJECTORY BASED DROPLET COLLISION MODEL FOR SPRAY MODELING Ö. O. TAŞKIRAN a,* and M. ERGENEMAN b a Domestic Research Division, Inventory Control Center, Golcuk, Kocaeli b Automotive Division, Faculty

More information

Liquid Jet Breakup at Low Weber Number: A Survey

Liquid Jet Breakup at Low Weber Number: A Survey International Journal of Engineering Research and Technology. ISSN 0974-3154 Volume 6, Number 6 (2013), pp. 727-732 International Research Publication House http://www.irphouse.com Liquid Jet Breakup at

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

JET AND DROPLET BREAKUP MODELLING APPROACHES

JET AND DROPLET BREAKUP MODELLING APPROACHES Journal of KONES Powertrain and Transport, Vol. 22, No. 3 2015 JET AND DROPLET BREAKUP MODELLING APPROACHES Łukasz Jan Kapusta, Piotr Jaworski Warsaw University of Technology, Institute of Heat Engineering

More information

EMISSION PREDICTIONS IN DIESEL ENGINES USING A NON-EQUILIBRIUM TURBULENCE DISSIPATION CORRECTION FOR THE k-ε MODEL

EMISSION PREDICTIONS IN DIESEL ENGINES USING A NON-EQUILIBRIUM TURBULENCE DISSIPATION CORRECTION FOR THE k-ε MODEL Twelfth International Multidimensional Engine Modeling Meeting at the SAE Congress March 3, 2002, Detroit, Michigan EMISSION PREDICTIONS IN DIESEL ENGINES USING A NON-EQUILIBRIUM TURBULENCE DISSIPATION

More information

MODELING ON THE BREAKUP OF VISCO-ELASTIC LIQUID FOR EFFERVESCENT ATOMIZATION

MODELING ON THE BREAKUP OF VISCO-ELASTIC LIQUID FOR EFFERVESCENT ATOMIZATION 1446 THERMAL SCIENCE, Year 2012, Vol. 16, No. 5, pp. 1446-1450 MODELING ON THE BREAKUP OF VISCO-ELASTIC LIQUID FOR EFFERVESCENT ATOMIZATION by Li-Juan QIAN * China Jiliang University, Hangzhou, China Short

More information

3. ATOMIZATION MODELS

3. ATOMIZATION MODELS 3. ATOMIZATION MODELS The atomization model supplies the initial conditions for spray computations, i.e., the drop sizes, velocities, temperatures etc., at the injector nozzle exit. The fundamental mechanisms

More information

Large Eddy Simulation of High Gas Density Effects in Fuel Sprays

Large Eddy Simulation of High Gas Density Effects in Fuel Sprays Large Eddy Simulation of High Gas Density Effects in Fuel Sprays Kaario, O. 1,*, Vuorinen, V. 1, Hulkkonen, T. 1, Keskinen, K. 1, Nuutinen, M. 1, Larmi, M. 1, and Tanner, F.X. 2 1 Aalto University School

More information

Evaluation of Liquid Fuel Spray Models for Hybrid RANS/LES and DLES Prediction of Turbulent Reactive Flows

Evaluation of Liquid Fuel Spray Models for Hybrid RANS/LES and DLES Prediction of Turbulent Reactive Flows Evaluation of Liquid Fuel Spray Models for Hybrid RANS/LES and DLES Prediction of Turbulent Reactive Flows by Ali Afshar A thesis submitted in conformity with the requirements for the degree of Masters

More information

IV Journeys in Multiphase Flows (JEM 2017) March 27-31, 2017, São Paulo, SP, Brazil Copyright 2017 by ABCM Paper ID: JEM

IV Journeys in Multiphase Flows (JEM 2017) March 27-31, 2017, São Paulo, SP, Brazil Copyright 2017 by ABCM Paper ID: JEM IV Journeys in Multiphase Flows (JEM 2017) March 27-31, 2017, São Paulo, SP, Brazil Copyright 2017 by ABCM Paper ID: JEM-2017-0013 CFD SIMULATION OF TWO-PHASE FLOW IN A VENTURI SCRUBBER: VALIDATION AND

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

Model comparison for single droplet fragmentation under varying accelerations

Model comparison for single droplet fragmentation under varying accelerations , 24th European Conference on Liquid Atomization and Spray Systems, Estoril, Portugal, September 2011 F.-O. Bartz 1, D. R. Guildenbecher 2, R. Schmehl 3, R. Koch 1, H.-J. Bauer 1 and P. E. Sojka 2 1: Institut

More information

Modelling the influence of the nozzle geometry on the primary breakup of diesel jets

Modelling the influence of the nozzle geometry on the primary breakup of diesel jets International journal of spray and combustion dynamics Volume. 6 Number.. 4 pages 3 46 3 Modelling the influence of the nozzle geometry on the primary breakup of diesel jets Oscar J. Soriano-Palao *, Martin

More information

An Improved Fuel Spray Breakup Model for Advanced Internal Combustion Engine Configurations. Chitra Sadanandan

An Improved Fuel Spray Breakup Model for Advanced Internal Combustion Engine Configurations. Chitra Sadanandan An Improved Fuel Spray Breakup Model for Advanced Internal Combustion Engine Configurations by Chitra Sadanandan A thesis submitted to the College of Engineering at Florida Institute of Technology in partial

More information

Effect of Primary Spray Characteristics on the Spray Generated by an Airblast Atomizer under High-Pressure Conditions

Effect of Primary Spray Characteristics on the Spray Generated by an Airblast Atomizer under High-Pressure Conditions ISS mericas, st nnual Conference on iquid tomization and Spray Systems, Orlando, Florida, May 8-8 Effect of Primary Spray Characteristics on the Spray Generated by an irblast tomizer under High-Pressure

More information

The effect of momentum flux ratio and turbulence model on the numerical prediction of atomization characteristics of air assisted liquid jets

The effect of momentum flux ratio and turbulence model on the numerical prediction of atomization characteristics of air assisted liquid jets ILASS Americas, 26 th Annual Conference on Liquid Atomization and Spray Systems, Portland, OR, May 204 The effect of momentum flux ratio and turbulence model on the numerical prediction of atomization

More information

Eulerian Two-Phase Flow CFD Simulation Using a Compressible and Equilibrium Eight- Equation Model. Y. Wang 1 and R. D. Reitz

Eulerian Two-Phase Flow CFD Simulation Using a Compressible and Equilibrium Eight- Equation Model. Y. Wang 1 and R. D. Reitz ILASS Americas 27th Annual Conference on Liquid Atomization and Spray Systems, Raleigh, NC, May 2015 Eulerian Two-Phase Flow CFD Simulation Using a Compressible and Equilibrium Eight- Equation Model Y.

More information

Liquid-Rocket Transverse Triggered Combustion Instability: Deterministic and Stochastic Analyses

Liquid-Rocket Transverse Triggered Combustion Instability: Deterministic and Stochastic Analyses Liquid-Rocket Transverse Triggered Combustion Instability: Deterministic and Stochastic Analyses by W. A. Sirignano Mechanical and Aerospace Engineering University of California, Irvine Collaborators:

More information

Modeling the influence of the nozzle flow on diesel spray atomization under high pressure injection conditions

Modeling the influence of the nozzle flow on diesel spray atomization under high pressure injection conditions Modeling the influence of the nozzle flow on diesel spray atomization under high pressure injection conditions J. O. Soriano 1*, M. Sommerfeld, A. Burkhardt 1 and U. Leuteritz 1 1 Siemens VDO Automotive,

More information

A Model to Predict Spray-tip Penetration for Time-varying Injection Profiles

A Model to Predict Spray-tip Penetration for Time-varying Injection Profiles ILASS Americas, 0 th Annual Conference on Liquid Atomization and Spray Systems, Chicago, IL, May 007 A Model to Predict Spray-tip Penetration for Time-varying Injection Profiles Neerav Abani* and Rolf

More information

Contents. Microfluidics - Jens Ducrée Physics: Laminar and Turbulent Flow 1

Contents. Microfluidics - Jens Ducrée Physics: Laminar and Turbulent Flow 1 Contents 1. Introduction 2. Fluids 3. Physics of Microfluidic Systems 4. Microfabrication Technologies 5. Flow Control 6. Micropumps 7. Sensors 8. Ink-Jet Technology 9. Liquid Handling 10.Microarrays 11.Microreactors

More information

An extension of dynamic droplet deformation models to secondary atomization

An extension of dynamic droplet deformation models to secondary atomization , 23rd Annual Conference on Liquid Atomization and Spray Systems, Brno, Czech Republic, September 2 F.-O. Bartz, R. Schmehl, R. Koch and H.-J. Bauer Institut für Thermische Strömungsmaschinen Karlsruher

More information

Detailed Numerical Simulation of Liquid Jet in Cross Flow Atomization: Impact of Nozzle Geometry and Boundary Condition

Detailed Numerical Simulation of Liquid Jet in Cross Flow Atomization: Impact of Nozzle Geometry and Boundary Condition ILASS-Americas 25th Annual Conference on Liquid Atomization and Spray Systems, Pittsburgh, PA, May 23 Detailed Numerical Simulation of Liquid Jet in Cross Flow Atomization: Impact of Nozzle Geometry and

More information

Numerical Investigation of the Fluid Flow around and Past a Circular Cylinder by Ansys Simulation

Numerical Investigation of the Fluid Flow around and Past a Circular Cylinder by Ansys Simulation , pp.49-58 http://dx.doi.org/10.1457/ijast.016.9.06 Numerical Investigation of the Fluid Flow around and Past a Circular Cylinder by Ansys Simulation Mojtaba Daneshi Department of Mechanical Engineering,

More information

Predicting Breakup Characteristics of Liquid Jets Disturbed by Practical Piezoelectric Devices

Predicting Breakup Characteristics of Liquid Jets Disturbed by Practical Piezoelectric Devices ILASS Americas 2th Annual Conference on Liquid Atomization and Spray Systems, Chicago, IL, May 27 Predicting Breakup Characteristics of Liquid Jets Disturbed by Practical Piezoelectric Devices M. Rohani,

More information

Reciprocating Internal Combustion Engines

Reciprocating Internal Combustion Engines Reciprocating Internal Combustion Engines Prof. Rolf D. Reitz Engine Research Center University of Wisconsin-Madison 014 Princeton-CEFRC Summer School on Combustion Course Length: 15 hrs (Mon.- Fri., June

More information

CFD Analysis of Gas Cooling

CFD Analysis of Gas Cooling As presented at: ILASS Americas, 20th Annual Conference on Liquid Atomization and Spray Systems, Chicago, IL, May 2007 J. S. Markus* Spraying Systems Deutschland GmbH Paul-Strähle Str. 10 73614 Schorndorf,

More information

Figure 11.1: A fluid jet extruded where we define the dimensionless groups

Figure 11.1: A fluid jet extruded where we define the dimensionless groups 11. Fluid Jets 11.1 The shape of a falling fluid jet Consider a circular orifice of a radius a ejecting a flux Q of fluid density ρ and kinematic viscosity ν (see Fig. 11.1). The resulting jet accelerates

More information

NUMERICAL INVESTIGATION ON THE EFFECT OF COOLING WATER SPRAY ON HOT SUPERSONIC JET

NUMERICAL INVESTIGATION ON THE EFFECT OF COOLING WATER SPRAY ON HOT SUPERSONIC JET Volume 119 No. 12 2018, 59-63 ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu ijpam.eu NUMERICAL INVESTIGATION ON THE EFFECT OF COOLING WATER SPRAY ON HOT SUPERSONIC JET Ramprasad T and Jayakumar

More information

LES Investigation of Fuel Effects on Lean Blow off (LBO) for a Realistic Two-Phase Flow Combustor

LES Investigation of Fuel Effects on Lean Blow off (LBO) for a Realistic Two-Phase Flow Combustor LES Investigation of Fuel Effects on Lean Blow off (LBO) for a Realistic Two-Phase Flow Combustor J.W. LABAHN, P. C. M A, L. E SCLAPE, M. I HME S TANFORD U NIVERSITY 2017 SIAM International Conference

More information

A KIVA-based Model for Liquid Jet in Cross Flow

A KIVA-based Model for Liquid Jet in Cross Flow ILASS-Americas nd Annual Conference on Liquid Atomization and Spray Systems, Cincinnati, OH, May 010 A KIVA-based Model for Liquid Jet in Cross Flow M. Behzad 1, A. Mashayek and N. Ashgriz Multiphase Flow

More information

NUMERICAL MODELING FOR ATOMIZATION OF COAXIAL LIQUID/GAS JETS

NUMERICAL MODELING FOR ATOMIZATION OF COAXIAL LIQUID/GAS JETS 90 Journal of Marine Science and Technology, Vol. 1, No. 4, pp. 90-99 (004) NUMERICAL MODELING FOR ATOMIZATION OF COAXIAL LIQUID/GAS JETS Stephen Gen-Ken Chuech*, Andrzej J. Przekwas**, and Chih-Yuan Wang***

More information

Cluster fusion in a high magnetic field

Cluster fusion in a high magnetic field Santa Fe July 28, 2009 Cluster fusion in a high magnetic field Roger Bengtson, Boris Breizman Institute for Fusion Studies, Fusion Research Center The University of Texas at Austin In collaboration with:

More information

Dynamics of Transient Liquid Injection:

Dynamics of Transient Liquid Injection: Dynamics of Transient Liquid Injection: K-H instability, vorticity dynamics, R-T instability, capillary action, and cavitation William A. Sirignano University of California, Irvine -- Round liquid columns

More information

ATOMISATION TECHNOLOGY TO MINIMIZE THE EFFECTS OF NOZZLE WEAR ON THE DROPLET SIZE

ATOMISATION TECHNOLOGY TO MINIMIZE THE EFFECTS OF NOZZLE WEAR ON THE DROPLET SIZE ILASS-Europe Darmstadt 11-13 September ATOMISATION TECHNOLOGY TO MINIMIZE THE EFFECTS OF NOZZLE WEAR ON THE DROPLET SIZE J. Kohlmann*, M. Schmidt*, G. Slowik** and M.Bürgermeister*** * Fachhochschule für

More information

Prototype Instabilities

Prototype Instabilities Prototype Instabilities David Randall Introduction Broadly speaking, a growing atmospheric disturbance can draw its kinetic energy from two possible sources: the kinetic and available potential energies

More information

Simulation of soot formation and analysis of the "Advanced soot model" parameters in an internal combustion engine

Simulation of soot formation and analysis of the Advanced soot model parameters in an internal combustion engine Simulation of soot formation and analysis of the "Advanced soot model" parameters in an internal combustion engine Marko Ban, dipl. ing.* Power Engineering Department Faculty of Mechanical Engineering

More information

Simulation of a Pressure Driven Droplet Generator

Simulation of a Pressure Driven Droplet Generator Simulation of a Pressure Driven Droplet Generator V. Mamet* 1, P. Namy 2, N. Berri 1, L. Tatoulian 1, P. Ehouarn 1, V. Briday 1, P. Clémenceau 1 and B. Dupont 1 1 DBV Technologies, 2 SIMTEC *84 rue des

More information

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t) IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common

More information

Viscous potential flow analysis of stress induced cavitation in an aperture flow

Viscous potential flow analysis of stress induced cavitation in an aperture flow cavitation-july9.tex Viscous potential flow analysis of stress induced cavitation in an aperture flow T. Funada, J. Wang and D. D. Joseph Department of Digital Engineering, Numazu College of Technology,

More information

The Bernoulli Equation

The Bernoulli Equation The Bernoulli Equation The most used and the most abused equation in fluid mechanics. Newton s Second Law: F = ma In general, most real flows are 3-D, unsteady (x, y, z, t; r,θ, z, t; etc) Let consider

More information

INTERNAL FLOW IN A Y-JET ATOMISER ---NUMERICAL MODELLING---

INTERNAL FLOW IN A Y-JET ATOMISER ---NUMERICAL MODELLING--- ILASS-Europe 2002 Zaragoza 9 11 September 2002 INTERNAL FLOW IN A Y-JET ATOMISER ---NUMERICAL MODELLING--- Z. Tapia, A. Chávez e-mail: ztapia@imp.mx Instituto Mexicano del Petróleo Blvd. Adolfo Ruiz Cortines

More information

CST Investigation on High Speed Liquid Jet using Computational Fluid Dynamics Technique

CST Investigation on High Speed Liquid Jet using Computational Fluid Dynamics Technique The 23 rd Conference of the Mechanical Engineering Network of Thailand November 4 7, 2009, Chiang Mai Investigation on High Speed Liquid Jet using Computational Fluid Dynamics Technique Wirapan Seehanam*,

More information

Visualization of high-speed gas jets and their airblast sprays of cross-injected liquid

Visualization of high-speed gas jets and their airblast sprays of cross-injected liquid Short communications Experiments in Fluids 27 (1999) 102 106 Springer-Verlag 1999 Visualization of high-speed gas jets and their airblast sprays of cross-injected liquid K. D. Kihm, T. K. Kim, S. Y. Son

More information

Dynamics of Secondary Breakup of Emulsified Fuel Drop

Dynamics of Secondary Breakup of Emulsified Fuel Drop 25 th ICDERS August 2 7, 2015 Leeds, UK Dynamics of Secondary Breakup of Emulsified Fuel Drop Oleksandr G. Girin Odessa National Maritime University Odessa, Ukraine 1 Introduction Diesel engines are used

More information

Mixing and Combustion in Dense Mixtures by William A. Sirignano and Derek Dunn-Rankin

Mixing and Combustion in Dense Mixtures by William A. Sirignano and Derek Dunn-Rankin Mixing and Combustion in Dense Mixtures by William A. Sirignano and Derek Dunn-Rankin At very high pressures and densities, what is different and what is similar about the processes of Injection and Atomization,

More information

CHARACTERISTICS OF ELLIPTIC CO-AXIAL JETS

CHARACTERISTICS OF ELLIPTIC CO-AXIAL JETS ELECTRIC POWER 2003 March 4-6, 2003 George R Brown Convention Center, Houston, TX EP 03 Session 07C: Fuels, Combustion and Advanced Cycles - Part II ASME - FACT Division CHARACTERISTICS OF ELLIPTIC CO-AXIAL

More information

Simulation analysis using CFD on vibration behaviors of circular cylinders subjected to free jets through narrow gaps in the vicinity of walls

Simulation analysis using CFD on vibration behaviors of circular cylinders subjected to free jets through narrow gaps in the vicinity of walls Fluid Structure Interaction V 85 Simulation analysis using CFD on vibration behaviors of circular cylinders subjected to free jets through narrow gaps in the vicinity of walls K. Fujita Osaka City University,

More information

Reciprocating Internal Combustion Engines

Reciprocating Internal Combustion Engines Reciprocating Internal Combustion Engines Prof. Rolf D. Reitz, Engine Research Center, University of Wisconsin-Madison 01 Princeton-CEFRC Summer Program on Combustion Course Length: 9 hrs (Wed., Thur.,

More information

A numerical study on the effects of cavitation on orifice flow

A numerical study on the effects of cavitation on orifice flow PHSICS OF FLUIDS, A numerical study on the effects of cavitation on orifice flow S. Dabiri, W. A. Sirignano, and D. D. Joseph, University of California, Irvine, California 9697, USA University of Minnesota,

More information

Computational Study of Sprays for the Development of a Monte Carlo Model

Computational Study of Sprays for the Development of a Monte Carlo Model 38th Dayton-Cincinnati Aerospace Sciences Symposium Computational Study of Sprays for the Development of a Monte Carlo Model Presenter: Murat Dinc West Virginia University Donald D. Gray West Virginia

More information

LINKING NOZZLE FLOW AND PRIMARY BREAKUP OF HIGH PRESSURE DIESEL JETS USING CFD. Oscar J. Soriano, Rossella Rotondi

LINKING NOZZLE FLOW AND PRIMARY BREAKUP OF HIGH PRESSURE DIESEL JETS USING CFD. Oscar J. Soriano, Rossella Rotondi Paper ID ILASS08-2-6 ILASS 2008 Sep. 8-10,2008, Como Lake, Italy LINKING NOZZLE FLOW AND PRIMARY BREAKUP OF HIGH PRESSURE DIESEL JETS USING CFD Oscar J. Soriano, Rossella Rotondi Diesel Systems, Hydraulic

More information

Numerical Investigation of Ignition Delay in Methane-Air Mixtures using Conditional Moment Closure

Numerical Investigation of Ignition Delay in Methane-Air Mixtures using Conditional Moment Closure 21 st ICDERS July 23-27, 27 Poitiers, France Numerical Investigation of Ignition Delay in Methane-Air Mixtures using Conditional Moment Closure Ahmad S. El Sayed, Cécile B. Devaud Department of Mechanical

More information

Numerical Simulation of Unsteady Nozzle Flow and Spray Formation under Diesel Engine Conditions

Numerical Simulation of Unsteady Nozzle Flow and Spray Formation under Diesel Engine Conditions Numerical Simulation of Unsteady Nozzle Flow and Spray Formation under Diesel Engine Conditions Mikhail Konstantinov * & Claus Wagner German Aerospace Center (DLR) Institute for Aerodynamics and Flow Technology

More information

LASER BASED DIAGNOSTIC SYSTEM FOR SPRAY MEASUREMENTS

LASER BASED DIAGNOSTIC SYSTEM FOR SPRAY MEASUREMENTS Journal of KONES Powertrain and Transport, Vol. 22, No. 3 2015 LASER BASED DIAGNOSTIC SYSTEM FOR SPRAY MEASUREMENTS Łukasz Jan Kapusta, Piotr Jaworski, Andrzej Teodorczyk Warsaw University of Technology,

More information

Modeling of Direct Gas Injection in Internal Combustion Engines

Modeling of Direct Gas Injection in Internal Combustion Engines Modeling of Direct Gas Injection in Internal Combustion Engines A. Y. Deshmukh 1, C. Giefer 1, M. Khosravi 2, H. Pitsch 1 1 Institute for Combustion Technology, RWTH Aachen University 2 Ford Research and

More information

Nonlinear shape evolution of immiscible two-phase interface

Nonlinear shape evolution of immiscible two-phase interface Nonlinear shape evolution of immiscible two-phase interface Francesco Capuano 1,2,*, Gennaro Coppola 1, Luigi de Luca 1 1 Dipartimento di Ingegneria Industriale (DII), Università di Napoli Federico II,

More information

EULERIAN LAGRANGIAN RANS MODEL SIMULATIONS OF THE NIST TURBULENT METHANOL SPRAY FLAME

EULERIAN LAGRANGIAN RANS MODEL SIMULATIONS OF THE NIST TURBULENT METHANOL SPRAY FLAME Combust. Sci. Technol., 187: 111 1138, 215 Copyright Taylor & Francis Group, LLC ISSN: 1-222 print / 1563-521X online DOI: 1.18/1222.215.119616 EULERIAN LAGRANGIAN RANS MODEL SIMULATIONS OF THE NIST TURBULENT

More information

Emulsion Processing - Homogenization -

Emulsion Processing - Homogenization - Emulsion Processing - Homogenization - Jochen Weiss *Food Structure and Functionality Laboratories Department of Food Science & Biotechnology University of Hohenheim Garbenstrasse 21, 70599 Stuttgart,

More information

Effect of Viscosity on the Breakup Length of Liquid Sheets Formed by Splash Plate Nozzles. University of Toronto, Toronto, ON M5S 3G8

Effect of Viscosity on the Breakup Length of Liquid Sheets Formed by Splash Plate Nozzles. University of Toronto, Toronto, ON M5S 3G8 ILASS Americas, 19 th Annual Conference on Liquid Atomization and Spray Systems, Toronto, Canada, May 2006 Effect of Viscosity on the Breakup Length of Liquid Sheets Formed by Splash Plate Nozzles M. Ahmed

More information

Size-velocity pdfs for Drop Fragments Formed via Bag Breakup

Size-velocity pdfs for Drop Fragments Formed via Bag Breakup ILASS-Americas 29th Annual Conference on Liquid Atomization and Spray Systems, Atlanta, GA, May 2017 Size-velocity pdfs for Drop Fragments Formed via Bag Breakup C.M.L. White,* G. Sondgeroth,+ W. Shang,+

More information

Linear analysis of three-dimensional instability of non-newtonian liquid jets

Linear analysis of three-dimensional instability of non-newtonian liquid jets J. Fluid Mech. (2006), vol. 559, pp. 451 459. c 2006 Cambridge University Press doi:10.1017/s0022112006000413 Printed in the United Kingdom 451 Linear analysis of three-dimensional instability of non-newtonian

More information

PAGE. 4.e' ce.u W -e1- eo- oeo-otwcata V'_2'e ~~ ~~~~~~ 1,A,.,,..Alt 3 REPORT TYPE AND DATES COVERED

PAGE. 4.e' ce.u W -e1- eo- oeo-otwcata V'_2'e ~~ ~~~~~~ 1,A,.,,..Alt 3 REPORT TYPE AND DATES COVERED MASTER COPY KI-I'l 'rills COPY FOR 10.'PRODUCTION ILUII'O;si:S AD -.A263 6.50 Il IiIIUIa-leIslsh!Iile!II!IEI 11,110 PAGE ~~ ~~~~~~ Form mon Approved....M. No07088 4.e' ce.u W -e1- eo- oeo-otwcata V'_2'e

More information

Simulation of Liquid Jet Breakup Process by Three-Dimensional Incompressible SPH Method

Simulation of Liquid Jet Breakup Process by Three-Dimensional Incompressible SPH Method Seventh International Conference on Computational Fluid Dynamics (ICCFD7), Big Island, Hawaii, July 9-13, 212 ICCFD7-291 Simulation of Liquid Jet Breakup Process by Three-Dimensional Incompressible SPH

More information

William A. Sirignano Mechanical and Aerospace Engineering University of California, Irvine

William A. Sirignano Mechanical and Aerospace Engineering University of California, Irvine Combustion Instability: Liquid-Propellant Rockets and Liquid-Fueled Ramjets William A. Sirignano Mechanical and Aerospace Engineering University of California, Irvine Linear Theory Nonlinear Theory Nozzle

More information

V. MODELING, SIMILARITY, AND DIMENSIONAL ANALYSIS To this point, we have concentrated on analytical methods of solution for fluids problems.

V. MODELING, SIMILARITY, AND DIMENSIONAL ANALYSIS To this point, we have concentrated on analytical methods of solution for fluids problems. V. MODELING, SIMILARITY, AND DIMENSIONAL ANALYSIS To this point, we have concentrated on analytical methods of solution for fluids problems. However, analytical methods are not always satisfactory due

More information

Chapter 3 Bernoulli Equation

Chapter 3 Bernoulli Equation 1 Bernoulli Equation 3.1 Flow Patterns: Streamlines, Pathlines, Streaklines 1) A streamline, is a line that is everywhere tangent to the velocity vector at a given instant. Examples of streamlines around

More information

Fluid Mechanics Theory I

Fluid Mechanics Theory I Fluid Mechanics Theory I Last Class: 1. Introduction 2. MicroTAS or Lab on a Chip 3. Microfluidics Length Scale 4. Fundamentals 5. Different Aspects of Microfluidcs Today s Contents: 1. Introduction to

More information

Flow Focusing Droplet Generation Using Linear Vibration

Flow Focusing Droplet Generation Using Linear Vibration Flow Focusing Droplet Generation Using Linear Vibration A. Salari, C. Dalton Department of Electrical & Computer Engineering, University of Calgary, Calgary, AB, Canada Abstract: Flow focusing microchannels

More information

Numerical Study of Laminar Annular Two-Phase Flow in Effervescent Atomizers

Numerical Study of Laminar Annular Two-Phase Flow in Effervescent Atomizers ILASS Americas 28th Annual Conference on Liquid Atomization and Spray Systems, Dearborn, MI, May 2016 Numerical Study of Laminar Annular Two-Phase Flow in Effervescent Atomizers C.K. Mohapatra and M.A.

More information

The Effect of Nozzle Height on Cooling Heat Transfer from a Hot Steel Plate by an Impinging Liquid Jet

The Effect of Nozzle Height on Cooling Heat Transfer from a Hot Steel Plate by an Impinging Liquid Jet , pp. 704 709 The Effect of Nozzle Height on Cooling Heat Transfer from a Hot Steel Plate by an Impinging Liquid Jet Piljong LEE, Haewon CHOI 1) and Sunghong LEE 2) Technical Research Center, POSCO, Pohang

More information

Review of Fluid Mechanics

Review of Fluid Mechanics Chapter 3 Review of Fluid Mechanics 3.1 Units and Basic Definitions Newton s Second law forms the basis of all units of measurement. For a particle of mass m subjected to a resultant force F the law may

More information

Spray formation with complex fluids

Spray formation with complex fluids Journal of Physics: Conference Series Spray formation with complex fluids To cite this article: S Lustig and M Rosen 2011 J. Phys.: Conf. Ser. 296 012019 View the article online for updates and enhancements.

More information

Fluid Dynamics Exercises and questions for the course

Fluid Dynamics Exercises and questions for the course Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r

More information

Separating Cloud Forming Nuclei from Interstitial Aerosol

Separating Cloud Forming Nuclei from Interstitial Aerosol Chapter 16 Separating Cloud Forming Nuclei from Interstitial Aerosol Gourihar Kulkarni Additional information is available at the end of the chapter http://dx.doi.org/10.5772/50589 1. Introduction Our

More information

Simulating the effect of in-nozzle cavitation on liquid atomisation using a three-phase model

Simulating the effect of in-nozzle cavitation on liquid atomisation using a three-phase model Simulating the effect of in-nozzle cavitation on liquid atomisation using a three-phase model M.G. Mithun*,1, P. Koukouvnis 1, I. K. Karathanassis 1, M. Gavaises 1 1 City, University of London, UK Abstract

More information

LIQUID FILM THICKNESS OF OSCILLATING FLOW IN A MICRO TUBE

LIQUID FILM THICKNESS OF OSCILLATING FLOW IN A MICRO TUBE Proceedings of the ASME/JSME 2011 8th Thermal Engineering Joint Conference AJTEC2011 March 13-17, 2011, Honolulu, Hawaii, USA AJTEC2011-44190 LIQUID FILM THICKNESS OF OSCILLATING FLOW IN A MICRO TUBE Youngbae

More information

ANSYS Advanced Solutions for Gas Turbine Combustion. Gilles Eggenspieler 2011 ANSYS, Inc.

ANSYS Advanced Solutions for Gas Turbine Combustion. Gilles Eggenspieler 2011 ANSYS, Inc. ANSYS Advanced Solutions for Gas Turbine Combustion Gilles Eggenspieler ANSYS, Inc. 1 Agenda Steady State: New and Existing Capabilities Reduced Order Combustion Models Finite-Rate Chemistry Models Chemistry

More information

Experimental Study of Injection of Oil-Water Emulsions into a Crossflow

Experimental Study of Injection of Oil-Water Emulsions into a Crossflow ILASS Americas, 24 th Annual Conference on Liquid Atomization and Spray Systems, San Antonio, TX, May 2012 Experimental Study of Injection of Oil-Water Emulsions into a Crossflow C.D. Bolszo, G.A. Gomez

More information

ME3250 Fluid Dynamics I

ME3250 Fluid Dynamics I ME3250 Fluid Dynamics I Section I, Fall 2012 Instructor: Prof. Zhuyin Ren Department of Mechanical Engineering University of Connecticut Course Information Website: http://www.engr.uconn.edu/~rzr11001/me3250_f12/

More information

A Novel FEM Method for Predicting Thermoacoustic Combustion Instability

A Novel FEM Method for Predicting Thermoacoustic Combustion Instability Excerpt from the Proceedings of the COMSOL Conference 009 Milan A Novel FEM Method for Predicting Thermoacoustic Combustion Instability G. Campa *, S.M. Camporeale Politecnico di Bari * campa@imedado.poliba.it,

More information

Application of the immersed boundary method to simulate flows inside and outside the nozzles

Application of the immersed boundary method to simulate flows inside and outside the nozzles Application of the immersed boundary method to simulate flows inside and outside the nozzles E. Noël, A. Berlemont, J. Cousin 1, T. Ménard UMR 6614 - CORIA, Université et INSA de Rouen, France emeline.noel@coria.fr,

More information

XXXVIII Meeting of the Italian Section of the Combustion Institute

XXXVIII Meeting of the Italian Section of the Combustion Institute Coupling a Helmholtz solver with a Distributed Flame Transfer Function (DFTF) to study combustion instability of a longitudinal combustor equipped with a full-scale burner D. Laera*, S.M. Camporeale* davide.laera@poliba.it

More information

A Study of Grid Resolution and SGS Models for LES under Non-reacting Spray Conditions

A Study of Grid Resolution and SGS Models for LES under Non-reacting Spray Conditions ILASS Americas, 25 th Annual Conference on Liquid Atomization and Spray Systems, Pittsburgh, PA, May 2013 A Study of Grid Resolution and SGS Models for LES under Non-reacting Spray Conditions Q. Xue 1*,

More information

FLUID MECHANICS. Chapter 3 Elementary Fluid Dynamics - The Bernoulli Equation

FLUID MECHANICS. Chapter 3 Elementary Fluid Dynamics - The Bernoulli Equation FLUID MECHANICS Chapter 3 Elementary Fluid Dynamics - The Bernoulli Equation CHAP 3. ELEMENTARY FLUID DYNAMICS - THE BERNOULLI EQUATION CONTENTS 3. Newton s Second Law 3. F = ma along a Streamline 3.3

More information

R13 SET - 1 '' ''' '' ' '''' Code No RT21033

R13 SET - 1 '' ''' '' ' '''' Code No RT21033 SET - 1 II B. Tech I Semester Supplementary Examinations, June - 2015 THERMODYNAMICS (Com. to ME, AE, AME) Time: 3 hours Max. Marks: 70 Note: 1. Question Paper consists of two parts (Part-A and Part-B)

More information

Parametric Investigation of the Common Geometry Shapes for Added Mass Calculation

Parametric Investigation of the Common Geometry Shapes for Added Mass Calculation Parametric Investigation of the Common Geometry Shapes for Added Mass Calculation Afsoun Koushesh* and Jin Lee Department of Mechanical Engineering, American University of Sharjah, Sharjah, UAE *Corresponding

More information

Efficient Numerical Calculation of Evaporating Sprays in Combustion Chamber Flows

Efficient Numerical Calculation of Evaporating Sprays in Combustion Chamber Flows 51-1 Efficient Numerical Calculation of Evaporating Sprays in Combustion Chamber Flows R. Schmehl, G. Klose, G. Maier and S. Wittig Lehrstuhl und Institut für Thermische Strömungsmaschinen Universität

More information

Temperature distribution and heat flow across the combustion chamber wall.

Temperature distribution and heat flow across the combustion chamber wall. ΜΕΤΑΔΟΣΗ ΘΕΡΜΟΤΗΤΑΣ ΣΤΟΝ ΚΥΛΙΝΔΡΟ (J.B. Heywood: Internal Combustion Engine Fundamentals McGraw Hill 1988) Temperature distribution and heat flow across the combustion chamber wall. Throughout each engine

More information

Correlating Results from Numerical Simulation to SLIPI-based Measurements for a non-combusting Diesel Spray

Correlating Results from Numerical Simulation to SLIPI-based Measurements for a non-combusting Diesel Spray ICLASS 2012, 12 th Triennial International Conference on Liquid Atomization and Spray Systems, Heidelberg, Germany, September 2-6, 2012 Correlating Results from Numerical Simulation to SLIPI-based Measurements

More information

COMPUTATIONAL METHODS FOR THE INVESTIGATION OF LIQUID DROP PHENOMENA IN EXTERNAL GAS FLOWS

COMPUTATIONAL METHODS FOR THE INVESTIGATION OF LIQUID DROP PHENOMENA IN EXTERNAL GAS FLOWS Michigan Technological University Digital Commons @ Michigan Tech Dissertations, Master's Theses and Master's Reports 2016 COMPUTATIONAL METHODS FOR THE INVESTIGATION OF LIQUID DROP PHENOMENA IN EXTERNAL

More information

Implementing Fundamental Pharmaceutical Science and Materials/Engineer Expertise in Scale-up

Implementing Fundamental Pharmaceutical Science and Materials/Engineer Expertise in Scale-up Implementing Fundamental Pharmaceutical Science and Materials/Engineer Expertise in Scale-up 2 nd FDA/PQRI Conference on Advancing Product Quality Session: The Science of Tech Transfer/Scale-up North Bethesda,

More information

A high resolution collision algorithm for anisotropic particle populations

A high resolution collision algorithm for anisotropic particle populations A high resolution collision algorithm for anisotropic particle populations Philipp Pischke 1 and Reinhold Kneer1 1 Institute of Heat and Mass Transfer, RWTH Aachen University July 2014 Collision algorithm,

More information

Chapter 7 Separation of Particles from a Gas

Chapter 7 Separation of Particles from a Gas Chapter 7 Separation of Particles from a Gas For either gas cleaning (removal of dusts) or recovery of particulate products Separation Mechanisms Sedimentation : Settling chamber, centrifuge Migration

More information

NUMERICAL SIMULATION OF LDI COMBUSTOR WITH DISCRETE-JET SWIRLERS USING RE-STRESS MODEL IN THE KIVA CODE

NUMERICAL SIMULATION OF LDI COMBUSTOR WITH DISCRETE-JET SWIRLERS USING RE-STRESS MODEL IN THE KIVA CODE NUMERICAL SIMULATION OF LDI COMBUSTOR WITH DISCRETE-JET SWIRLERS USING RE-STRESS MODEL IN THE KIVA CODE S. L. Yang, C. Y. Teo, and Y. K. Siow Department of Mechanical Engineering Engineering Mechanics

More information

For example an empty bucket weighs 2.0kg. After 7 seconds of collecting water the bucket weighs 8.0kg, then:

For example an empty bucket weighs 2.0kg. After 7 seconds of collecting water the bucket weighs 8.0kg, then: Hydraulic Coefficient & Flow Measurements ELEMENTARY HYDRAULICS National Certificate in Technology (Civil Engineering) Chapter 3 1. Mass flow rate If we want to measure the rate at which water is flowing

More information

Simulation of Condensing Compressible Flows

Simulation of Condensing Compressible Flows Simulation of Condensing Compressible Flows Maximilian Wendenburg Outline Physical Aspects Transonic Flows and Experiments Condensation Fundamentals Practical Effects Modeling and Simulation Equations,

More information