
Another important type of atomizer is the pressureswirl atomizer, sometimes referred to by the gasturbine community as a simplex atomizer. This type of atomizer accelerates the liquid through nozzles known as swirl ports into a central swirl chamber. The swirling liquid pushes against the walls of the swirl chamber and develops a hollow air core. It then emerges from the orifice as a thinning sheet, which is unstable, breaking up into ligaments and droplets. The pressureswirl atomizer is very widely used for liquidfuel combustion in gas turbines, oil furnaces, and directinjection sparkignited automobile engines. The transition from internal injector flow to fullydeveloped spray can be divided into three steps: film formation, sheet breakup, and atomization. A sketch of how this process is thought to occur is shown in Figure 15.9.5.
The interaction between the air and the sheet is not well understood. It is generally accepted that an aerodynamic instability causes the sheet to break up. The mathematical analysis below assumes that KelvinHelmholtz waves grow on the sheet and eventually break the liquid into ligaments. It is then assumed that the ligaments break up into droplets due to varicose instability. Once the liquid droplets are formed, the spray evolution is determined by drag, collision, coalescence, and secondary breakup.
The pressureswirl atomizer model used in ANSYS FLUENT is called the Linearized Instability Sheet Atomization (LISA) model of Schmidt et al. [ 308]. The LISA model is divided into two stages:
Both parts of the model are described below.
Film Formation
The centrifugal motion of the liquid within the injector creates an air core surrounded by a liquid film. The thickness of this film, , is related to the mass flow rate by
where is the injector exit diameter, and is the effective mass flow rate, which is defined by Equation 15.95 . The other unknown in Equation 15.923 is , the axial component of velocity at the injector exit. This quantity depends on internal details of the injector and is difficult to calculate from first principles. Instead, the approach of Han et al. [ 122] is used. The total velocity is assumed to be related to the injector pressure by
where is the velocity coefficient. Lefebvre [ 186] has noted that is a function of the injector design and injection pressure. If the swirl ports are treated as nozzles and if it is assumed that the dominant portion of the pressure drop occurs at those ports, is the expression for the discharge coefficient ( ). For singlephase nozzles with sharp inlet corners and ratios of 4, a typical value is 0.78 or less [ 193]. If the nozzles are cavitating, the value of may be as low as 0.61. Hence, 0.78 should be a practical upper bound for . Reducing by 10% to 0.7 approximates the effect of other momentum losses on the discharge coefficient.
Physical limits on require that it be less than unity from conservation of energy, yet be large enough to permit sufficient mass flow. To guarantee that the size of the air core is nonnegative, the following expression is used for :
Assuming that is known, Equation 15.924 can be used to find . Once is determined, is found from
where is the spray angle, which is assumed to be known. At this point, the thickness and axial component of the liquid film are known at the injector exit. The tangential component of velocity ( ) is assumed to be equal to the radial velocity component of the liquid sheet downstream of the nozzle exit. The axial component of velocity is assumed to remain constant.
Sheet Breakup and Atomization
The pressureswirl atomizer model includes the effects of the surrounding gas, liquid viscosity, and surface tension on the breakup of the liquid sheet. Details of the theoretical development of the model are given in Senecal et al. [ 310] and are only briefly presented here. For a more robust implementation, the gasphase velocity is neglected in calculating the relative liquidgas velocity and is instead set by you. This avoids having the injector parameters depend too heavily on the usually underresolved gasphase velocity field very near the injection location.
The model assumes that a twodimensional, viscous, incompressible liquid sheet of thickness moves with velocity through a quiescent, inviscid, incompressible gas medium. The liquid and gas have densities of and , respectively, and the viscosity of the liquid is . A coordinate system is used that moves with the sheet, and a spectrum of infinitesimal wavy disturbances of the form
is imposed on the initially steady motion. The spectrum of disturbances results in fluctuating velocities and pressures for both the liquid and the gas. In Equation 15.927, is the initial wave amplitude, is the wave number, and is the complex growth rate. The most unstable disturbance has the largest value of , denoted here by , and is assumed to be responsible for sheet breakup. Thus, it is desired to obtain a dispersion relation from which the most unstable disturbance can be calculated as a function of wave number.
Squire [ 335], Li and Tankin [ 192], and Hagerty and Shea [ 118] have shown that two solutions, or modes, exist that satisfy the governing equations subject to the boundary conditions at the upper and lower interfaces. The first solution, called the sinuous mode, has waves at the upper and lower interfaces in phase. The second solution is called the varicose mode which has the waves at the upper and lower interfaces radians out of phase. It has been shown by numerous authors (e.g., Senecal et. al. [ 310]) that the sinuous mode dominates the growth of varicose waves for low velocities and low gastoliquid density ratios. In addition, it can be shown that the sinuous and varicose modes become indistinguishable for highvelocity flows. As a result, the atomization model in ANSYS FLUENT is based upon the growth of sinuous waves on the liquid sheet.
As derived in Senecal et al. [ 310], the dispersion relation for the sinuous mode is given by
where and .
Above a critical Weber number of = 27/16 (based on the liquid velocity, gas density, and sheet halfthickness), the fastestgrowing waves are short. For , the wavelengths are long compared to the sheet thickness. The speed of modern high pressure fuel injection systems is high enough to ensure that the film Weber number is well above this critical limit.
An orderofmagnitude analysis using typical values shows that the terms of second order in viscosity can be neglected in comparison to the other terms in Equation 15.928. Using this assumption, Equation 15.928 reduces to
For waves that are long compared with the sheet thickness, a mechanism of sheet disintegration proposed by Dombrowski and Johns [ 74] is adopted. For long waves, ligaments are assumed to form from the sheet breakup process once the unstable waves reach a critical amplitude. If the surface disturbance has reached a value of at breakup, a breakup time, , can be evaluated:
where , the maximum growth rate, is found by numerically maximizing Equation 15.929 as a function of . The maximum is found using a binary search that checks the sign of the derivative. The sheet breaks up and ligaments will be formed at a length given by
where the quantity is an empirical sheet constant that you must specify. The default value of 12 was obtained theoretically by Weber [ 369] for liquid jets. Dombrowski and Hooper [ 73] showed that a value of 12 for the sheet constant agreed favorably with experimental sheet breakup lengths over a range of Weber numbers from 2 to 200.
The diameter of the ligaments formed at the point of breakup can be obtained from a mass balance. If it is assumed that the ligaments are formed from tears in the sheet twice per wavelength, the resulting diameter is given by
where is the wave number corresponding to the maximum growth rate, . The ligament diameter depends on the sheet thickness, which is a function of the breakup length. The film thickness is calculated from the breakup length and the radial distance from the center line to the midline of the sheet at the atomizer exit, :
This mechanism is not used for waves that are short compared to the sheet thickness. For short waves, the ligament diameter is assumed to be linearly proportional to the wavelength that breaks up the sheet,
where , or the ligament constant, is equal to 0.5 by default.
In either the longwave or the shortwave case, the breakup from ligaments to droplets is assumed to behave according to Weber's [ 369] analysis for capillary instability.
Here, is the Ohnesorge number which is a combination of the Reynolds number and the Weber number (see Section 15.10.2 for more details about Oh). Once has been determined from Equation 15.935, it is assumed that this droplet diameter is the most probable droplet size of a RosinRammler distribution with a spread parameter of 3.5 and a default dispersion angle of 6 (which can be modified in the GUI). The choice of spread parameter and dispersion angle is based on past modeling experience [ 307]. It is important to note that the spray cone angle must be specified by you when using this model.