FLUENT offers two droplet breakup models: the Taylor analogy breakup (TAB) model and the wave model. The TAB model is recommended for low-Weber-number injections and is well suited for low-speed sprays into a standard atmosphere. For Weber numbers greater than 100, the wave model is more applicable. The wave model is popular for use in high-speed fuel-injection applications. Details for each model are provided below.
Taylor Analogy Breakup (TAB) Model
The Taylor analogy breakup (TAB) model is a classic method for calculating droplet breakup, which is applicable to many engineering sprays. This method is based upon Taylor's analogy [ 369] between an oscillating and distorting droplet and a spring mass system. Table 22.7.1 illustrates the analogous components.
The resulting TAB model equation set, which governs the oscillating and distorting droplet, can be solved to determine the droplet oscillation and distortion at any given time. As described in detail below, when the droplet oscillations grow to a critical value the "parent'' droplet will break up into a number of smaller "child'' droplets. As a droplet is distorted from a spherical shape, the drag coefficient changes. A drag model that incorporates the distorting droplet effects is available in FLUENT. See Section 22.6 for details.
Use and Limitations
The TAB model is best for low-Weber-number sprays. Extremely high-Weber-number sprays result in shattering of droplets, which is not described well by the spring-mass analogy.
The equation governing a damped, forced oscillator is [ 270]
where is the displacement of the droplet equator from its spherical (undisturbed) position. The coefficients of this equation are taken from Taylor's analogy:
where and are the discrete phase and continuous phase densities, is the relative velocity of the droplet, is the undisturbed droplet radius, is the droplet surface tension, and is the droplet viscosity. The dimensionless constants , , and will be defined later.
The droplet is assumed to break up if the distortion grows to a critical ratio of the droplet radius. This breakup requirement is given as
where is a constant equal to 0.5 if breakup is assumed to occur when the distortion is equal to the droplet radius, i.e., the north and south poles of the droplet meet at the droplet center. This implicitly assumes that the droplet is undergoing only one (fundamental) oscillation mode. Equation 22.7-8 is nondimensionalized by setting and substituting the relationships in Equations 22.7-9- 22.7-11:
where breakup now occurs for . For under-damped droplets, the equation governing can easily be determined from Equation 22.7-13 if the relative velocity is assumed to be constant:
In Equation 22.7-14, is the relative velocity between the droplet and the gas phase and We is the droplet Weber number, a dimensionless parameter defined as the ratio of aerodynamic forces to surface tension forces. The droplet oscillation frequency is represented by . The default value of is 0, based upon the work of Liu et al. [ 218]. The constants have been chosen to match experiments and theory [ 190]:
If Equation 22.7-14 is solved for all droplets, those with are assumed to break up. The size and velocity of the new child droplets must be determined.
Size of Child Droplets
The size of the child droplets is determined by equating the energy of the parent droplet to the combined energy of the child droplets. The energy of the parent droplet is [ 270]
where is the ratio of the total energy in distortion and oscillation to the energy in the fundamental mode, of the order ( ). The child droplets are assumed to be nondistorted and nonoscillating. Thus, the energy of the child droplets can be shown to be
where is the Sauter mean radius of the droplet size distribution. can be found by equating the energy of the parent and child droplets (i.e., Equations 22.7-21 and 22.7-22), setting , and :
Once the size of the child droplets is determined, the number of child droplets can easily be determined by mass conservation.
Velocity of Child Droplets
The TAB model allows for a velocity component normal to the parent droplet velocity to be imposed upon the child droplets. When breakup occurs, the equator of the parent droplet is traveling at a velocity of . Therefore, the child droplets will have a velocity normal to the parent droplet velocity given by
where is a constant of order (1).
To model droplet breakup, the TAB model first determines the amplitude for an undamped oscillation ) for each droplet at time step using the following:
According to Equation 22.7-25, breakup is possible only if the following condition is satisfied:
This is the limiting case, as damping will only reduce the chance of breakup. If a droplet fails the above criterion, breakup does not occur. The only additional calculations required then, are to update using a discretized form of Equation 22.7-14 and its derivative, which are both based on work done by O'Rourke and Amsden [ 270]:
All of the constants in these expressions are assumed to be constant throughout the time step.
If the criterion of Equation 22.7-26 is met, then breakup is possible. The breakup time, , must be determined to see if breakup occurs within the time step . The value of is set to the time required for oscillations to grow sufficiently large that the magnitude of the droplet distortion, , is equal to unity. The breakup time is determined under the assumption that the droplet oscillation is undamped for its first period. The breakup time is therefore the smallest root greater than of an undamped version of Equation 22.7-14:
If , then breakup will not occur during the current time step, and and are updated by Equations 22.7-27 and 22.7-28. The breakup calculation then continues with the next droplet. Conversely, if , then breakup will occur and the child droplet radii are determined by Equation 22.7-23. The number of child droplets, , is determined by mass conservation:
It is assumed that the child droplets are neither distorted nor oscillating; i.e., . The child droplets are represented by a number of child parcels which are created from the original parcel. These child parcels are distributed equally along the equator of the parent droplet in a plane normal to the parent relative velocity vector. The diameter of each of the child parcels is sampled from a Rosin Rammler distribution based on the Sauter mean radius (Equation 22.7-23) and a spread parameter of 3.5.
A velocity component normal to the relative velocity vector, with magnitude computed by Equation 22.7-24, is imposed upon the child droplets. It is decomposed at the equator into components pointing radially outward.
| A large number of child parcels ensures a smooth distribution of particle diameters and source terms which is needed when simulating, for example, evaporating sprays.
Wave Breakup Model
An alternative to the TAB model that is appropriate for high-Weber-number flows is the wave breakup model of Reitz [ 299], which considers the breakup of the droplets to be induced by the relative velocity between the gas and liquid phases. The model assumes that the time of breakup and the resulting droplet size are related to the fastest-growing Kelvin-Helmholtz instability, derived from the jet stability analysis described below. The wavelength and growth rate of this instability are used to predict details of the newly-formed droplets.
Use and Limitations
The wave model is appropriate for high-speed injections, where the Kelvin-Helmholtz instability is believed to dominate droplet breakup ( ). Because this breakup model can increase the number of computational parcels, you may wish to inject a modest number of droplets initially.
Jet Stability Analysis
The jet stability analysis described in detail by Reitz and Bracco [ 301] is presented briefly here. The analysis considers the stability of a cylindrical, viscous, liquid jet of radius issuing from a circular orifice at a velocity into a stagnant, incompressible, inviscid gas of density . The liquid has a density, , and viscosity, , and a cylindrical polar coordinate system is used which moves with the jet. An arbitrary infinitesimal axisymmetric surface displacement of the form
is imposed on the initially steady motion and it is thus desired to find the dispersion relation which relates the real part of the growth rate, , to its wave number, .
In order to determine the dispersion relation, the linearized equations for the hydrodynamics of the liquid are solved assuming wave solutions of the form
where and are the velocity potential and stream function, respectively, and are integration constants, and are modified Bessel functions of the first kind, , and is the liquid kinematic viscosity [ 299]. The liquid pressure is obtained from the inviscid part of the liquid equations. In addition, the inviscid gas equations can be solved to obtain the fluctuating gas pressure at :
where and are modified Bessel functions of the second kind and is the relative velocity between the liquid and the gas. The linearized boundary conditions are
which are mathematical statements of the liquid kinematic free surface condition, continuity of shear stress, and continuity of normal stress, respectively. Note that is the axial perturbation liquid velocity, is the radial perturbation liquid velocity, and is the surface tension. Also note that Equation 22.7-38 was obtained under the assumption that .
As described by Reitz [ 299], Equations 22.7-37 and 22.7-38 can be used to eliminate the integration constants and in Equations 22.7-34 and 22.7-35. Thus, when the pressure and velocity solutions are substituted into Equation 22.7-39, the desired dispersion relation is obtained:
As shown by Reitz [ 299], Equation 22.7-40 predicts that a maximum growth rate (or most unstable wave) exists for a given set of flow conditions. Curve fits of numerical solutions to Equation 22.7-40 were generated for the maximum growth rate, , and the corresponding wavelength, , and are given by Reitz [ 299]:
where is the Ohnesorge number and is the Taylor number. Furthermore, and are the liquid and gas Weber numbers, respectively, and is the Reynolds number.
In the wave model, breakup of droplet parcels is calculated by assuming that the radius of the newly-formed droplets is proportional to the wavelength of the fastest-growing unstable surface wave on the parent droplet. In other words,
where is a model constant set equal to 0.61 based on the work of Reitz [ 299]. Furthermore, the rate of change of droplet radius in the parent parcel is given by
where the breakup time, , is given by
and and are obtained from Equations 22.7-41 and 22.7-42, respectively. The breakup time constant, , is set to a value of 1.73 as recommended by Liu et al. [ 218]. Values of can range between 1 and 60, depending on the injector characterization.
In the wave model, mass is accumulated from the parent drop at a rate given by Equation 22.7-45 until the shed mass is equal to 5% of the initial parcel mass. At this time, a new parcel is created with a radius given by Equation 22.7-43. The new parcel is given the same properties as the parent parcel (i.e., temperature, material, position, etc.) with the exception of radius and velocity. The new parcel is given a component of velocity randomly selected in the plane orthogonal to the direction vector of the parent parcel, and the momentum of the parent parcel is adjusted so that momentum is conserved. The velocity magnitude of the new parcel is the same as the parent parcel.
You must also specify the model constants which determine how the gas phase interacts with the liquid droplets. For example, the breakup time constant B1 is the constant multiplying the time scale which determines how quickly the parcel will loose mass. Therefore, a larger number means that it takes longer for the particle to loose a given amount. A larger number for B1 in the context of interaction with the gas phase would mean that the interaction with the subgrid is less intense. B0 is the constant for the drop size and is generally taken to be 0.61.