
The plainorifice is the most common type of atomizer and the most simply made. However there is nothing simple about the physics of the internal nozzle flow and the external atomization. In the plainorifice atomizer model in FLUENT, the liquid is accelerated through a nozzle, forms a liquid jet and then breaks up to form droplets. This apparently simple process is dauntingly complex. The plain orifice may operate in three different regimes: singlephase, cavitating and flipped [ 348]. The transition between regimes is abrupt, producing dramatically different sprays. The internal regime determines the velocity at the orifice exit, as well as the initial droplet size and the angle of droplet dispersion. Diagrams of each case are shown in Figures 22.8.1, 22.8.2, and 22.8.3.
Internal Nozzle State
To accurately predict the spray characteristics, the plainorifice model in FLUENT must identify the correct state of the internal nozzle flow because the nozzle state has a tremendous effect on the external spray. Unfortunately, there is no established theory for determining the nozzle state. One must rely on empirical models obtained from experimental data. FLUENT uses several dimensionless parameters to determine the internal flow regime for the plainorifice atomizer model. These parameters and the decisionmaking process are summarized below.
A list of the parameters that control internal nozzle flow is given in Table 22.8.1. These parameters may be combined to form nondimensional characteristic lengths such as and , as well as nondimensional groups like the Reynolds number based on hydraulic "head'' ( ) and the cavitation parameter ( ).
The liquid flow often contracts in the nozzle, as can be seen in Figures 22.8.2 and 22.8.3. Nurick [ 266] found it helpful to use a coefficient of contraction ( ) to represent the reduction in the crosssectional area of the liquid jet. The coefficient of contraction is defined as the area of the stream of contracting liquid over the total crosssectional area of the nozzle. FLUENT uses Nurick's fit for the coefficient of contraction:
Here, is a theoretical constant equal to 0.611, which comes from potential flow analysis of flipped nozzles.
Coefficient of Discharge
Another important parameter for describing the performance of nozzles is the coefficient of discharge ( ). The coefficient of discharge is the ratio of the mass flow rate through the nozzle to the theoretical maximum mass flow rate:
where
is the
effective mass flow rate of the nozzle, defined by
Here, is the mass flow rate specified in the user interface, and is the difference between the azimuthal stop angle and the azimuthal start angle
as input by the user (see Section 22.12.1). Note that the mass flow rate input by the user should be for the appropriate start and stop angles, in other words the correct mass flow rate for the sector being modeled. Note also that for of , the effective mass flow rate is identical to the mass flow rate in the interface.
The cavitation number ( in Equation 22.82) is an essential parameter for predicting the inception of cavitation. The inception of cavitation is known to occur at a value of 1.9 for short, sharpedged nozzles. However, to include the effects of inlet rounding and viscosity, an empirical relationship is used:
Similarly, a critical value of where flip occurs is given by
If is greater than 0.05, then flip is deemed impossible and is set to 1.0.
The cavitation number, , is compared to the values of and to identify the nozzle state. The decision tree is shown in Figure 22.8.4. Depending on the state of the nozzle, a unique closure is chosen for the above equations.
For a singlephase nozzle ( ) [ 207], the coefficient of discharge is given by
where is the ultimate discharge coefficient, and is defined as
For a cavitating nozzle ( ) [ 266] the coefficient of discharge is determined from
For a flipped nozzle ( ) [ 266],
All of the nozzle flow equations are solved iteratively, along with the appropriate relationship for coefficient of discharge as given by the nozzle state. The nozzle state may change as the upstream or downstream pressures change. Once the nozzle state is determined, the exit velocity is calculated, and appropriate correlations for spray angle and initial droplet size distribution are determined.
Exit Velocity
For a singlephase nozzle, the estimate of exit velocity ( ) comes from the conservation of mass and the assumption of a uniform exit velocity:
For the cavitating nozzle, Schmidt and Corradini [ 323] have shown that the uniform exit velocity is not accurate. Instead, they derived an expression for a higher velocity over a reduced area:
This analytical relation is used for cavitating nozzles in FLUENT. For the case of flipped nozzles, the exit velocity is found from the conservation of mass and the value of the reduced flow area:
Spray Angle
The correlation for the spray angle ( ) comes from the work of Ranz [ 295]:
The spray angle for both singlephase and cavitating nozzles depends on the ratio of the gas and liquid densities and also the parameter . For flipped nozzles, the spray angle has a constant value.
The parameter , which you must specify, is thought to be a constant for a given nozzle geometry. The larger the value, the narrower the spray. Reitz [ 300] suggests the following correlation for :
The spray angle is sensitive to the internal flow regime of the nozzle. Hence, you may wish to choose smaller values of for cavitating nozzles than for singlephase nozzles. Typical values range from 4.0 to 6.0. The spray angle for flipped nozzles is a small, arbitrary value that represents the lack of any turbulence or initial disturbance from the nozzle.
Droplet Diameter Distribution
One of the basic characteristics of an injection is the distribution of drop size. For an atomizer, the droplet diameter distribution is closely related to the nozzle state. FLUENT's spray models use a twoparameter RosinRammler distribution, characterized by the most probable droplet size and a spread parameter. The most probable droplet size, is obtained in FLUENT from the Sauter mean diameter, [ 200]. For more information about the RosinRammler size distribution, see Section 22.12.1.
For singlephase nozzle flows, the correlation of Wu et al. [ 407] is used to calculate and relate the initial drop size to the estimated turbulence quantities of the liquid jet:
where , is the radial integral length scale at the jet exit based upon fullydeveloped turbulent pipe flow, and We is the Weber number, defined as
Here, is the droplet surface tension. For a more detailed discussion of droplet surface tension and the Weber number, see Section 22.7.2. For more information about mean particle diameters, see Section 22.16.8.
For cavitating nozzles, FLUENT uses a slight modification of Equation 22.818. The initial jet diameter used in Wu's correlation, , is calculated from the effective area of the cavitating orifice exit, and thus represents the effective diameter of the exiting liquid jet, . For an explanation of effective area of cavitating nozzles, please see Schmidt and Corradini [ 323].
The length scale for a cavitating nozzle is , where
For the case of the flipped nozzle, the initial droplet diameter is set to the diameter of the liquid jet:
where is defined as the most probable diameter.
The second parameter required to specify the droplet size distribution is the spread parameter, . The values for the spread parameter are chosen from past modeling experience and from a review of experimental observations. Table 22.8.2 lists the values of for the three nozzle states. The larger the value of the spread parameter, the narrower the droplet size distribution.
Since the correlations of Wu et al. provide the Sauter mean diameter, , these are converted to the most probable diameter, . Lefebvre [ 200] gives the most general relationship between the Sauter mean diameter and most probable diameter for a RosinRammler distribution. The simplified version for =3.5 is as follows:
At this point, the droplet size, velocity, and spray angle have been determined and the initialization of the injections is complete.