
Introduction
An alternative to the TAB model that is appropriate for highWebernumber flows is the wave breakup model of Reitz [ 287], 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 fastestgrowing KelvinHelmholtz instability, derived from the jet stability analysis described below. The wavelength and growth rate of this instability are used to predict details of the newlyformed droplets.
Use and Limitations
The wave model is appropriate for highspeed injections, where the KelvinHelmholtz 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 [ 289] 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 [ 287]. 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
and
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 15.1031 was obtained under the assumption that .
As described by Reitz [ 287], Equations 15.1030 and 15.1031 can be used to eliminate the integration constants and in Equations 15.1027 and 15.1028. Thus, when the pressure and velocity solutions are substituted into Equation 15.1032, the desired dispersion relation is obtained:
As shown by Reitz [ 287], Equation 15.1033 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 15.1033 were generated for the maximum growth rate, , and the corresponding wavelength, , and are given by Reitz [ 287]:
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.
Droplet Breakup
In the wave model, breakup of droplet parcels is calculated by assuming that the radius of the newlyformed droplets is proportional to the wavelength of the fastestgrowing 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 [ 287]. 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 15.1034 and 15.1035, respectively. The breakup time constant, , is set to a value of 1.73 as recommended by Liu et al. [ 205]. 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 15.1038 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 15.1036. 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.