CHARACTERISTICS OF BUBBLE SPLITTING IN A TIP VORTEX FLOW

Transcription

1 CAV3-GS-7- Fifth International Symposium on Cavitation (CAV3 Osaka, Japan, November -4, 3 CHARACTERISTICS OF BUBBLE SPLITTING IN A TIP VORTEX FLOW Jin-Keun Choi DYNAFLOW, INC. jkchoi@dynaflow-inc.com Georges L. Chahine DYNAFLOW, INC. glchahine@dynaflow-inc.com ABSTRACT The splitting behavior of a bubble in a tip vortex flow is studied using an axisymmetric boundary element method accounting for the effect of surrounding viscous flow. Bubble splitting behavior, bubble sizes before and after splitting, and emitted pressure signals are characterized in a wide range of cavitation numbers. The effect of the Reynolds number and the effect of initial bubble nucleus size are also investigated. These observations lead us to a simple bubble splitting model which can be applied within a practical spherical model to include bubble subdivision into smaller sub-bubbles. INTRODUCTION Recently, our axisymmetric free surface code, DYNAFS has been extended to simulate extreme bubble deformations including bubble splitting and subsequent behavior []. When the method is applied to a bubble in a tip vortex, it is found that, under defined conditions, the elongated bubble splits, after which violent reentrant jets develop [-3]. These reentrant jets result in pressure pulse emissions which are orders of magnitude higher than those due to spherical bubble dynamics. The method is based on coupling bubble dynamics models and an unsteady Reynolds Averaged Navier Stokes (RANS solver. In the coupled approach, surface averaged pressure (SAP spherical models [4,5] and non-spherical axisymmetric [-3] or fully three-dimensional models [4] using boundary element methods can be utilized to describe the bubble dynamics. Bubble behavior in a vortex occur in three phases; (a bubble capture by the vortex, (b interaction between the vortex and an initially quasi-spherical bubble, and (c the dynamics of elongated bubbles on the vortex axis [6]. In this study, the SAP spherical bubble dynamics model is used for the two phases while the axisymmetric method is utilized in the third phase. DYNAFS, is a boundary element potential flow solver. The code has been verified successfully for diverse types of fluid dynamic problems [7]. Recently, it has been extended to accommodate the ambient vortex flow field and to model the extreme deformation of bubbles including splitting [-3]. A summary of the numerical method and its validation against laboratory experiments can be found in Choi and Chahine [,]. A few past investigations considered bubble splitting between two parallel plates [8-]. These showed the hourglasslike bubble deformation experimentally and numerically. Ishida et al. [] reported the existence of a small peak in the pressure signal at the splitting. However, the numerical studies while able to simulate the bubble behavior until spitting, could not predict the dynamics beyond that point. In this paper, we study the characteristics of bubble splitting in a tip vortex flow. The bubble behavior is simulated by DYNAFS with the effect of surrounding viscous flow taken into account using one-way coupling []. Bubble size before and after the split and the emitted pressure signal peaks are characterized in a range of cavitation numbers varying from.5 to.8 and are related to the various bubble dynamic phases. The effect of the Reynolds number and the effect of initial bubble nucleus size are also investigated. Finally, a predictive model for the bubble splitting is suggested so that the SAP spherical model could be extended to provide simulations including bubble subdivision into smaller sub-bubbles. TIP VORTEX FLOW FIELDS The tip vortex flow field used in this study is computed with DF_UNCLE, a Reynolds Averaged Navier-Stokes (RANS flow solver. Two different scales are considered as follows: Hydrofoils : elliptic plan form; (a chord = m, half span =.5 m, (b chord =.5 m, half span =.75 m Inflow: U =.88m/s; (a Re=.44x 6, (b Re=.88x 6 After obtaining the three dimensional tip vortex flow fields, these are cast by averaging azimuthally into axisymmetric flow fields and input to the axisymmetric solver DYNAFS. This is achieved by finding the vortex center first and taking circumferential average of the field variable around various radii from the vortex center. The locus of the tip vortex center is found as the lowest pressure point in each grid plane roughly perpendicular to the stream direction. However, the circumferential average is taken on each plane exactly perpendicular to the inflow direction. One of the axisymmetric flow fields obtained by this procedure is shown in Figure.

2 r x p r x u Figure. Bubble behavior at σ=.5 as predicted by DYNAFS in time sequence from top left to bottom right. The tip vortex flows from upper left to lower right. Figure. The axisymmetric tip vortex flow field for Re=.88x6. (x,r are made dimensionless by the chord, p and u, by ρu and U, respectively. REPRESENTATIVE RESULTS In this section, the behavior of a bubble nucleus of radius 5 µm placed at. m upstream of the tip is used as an example. The cavitation number studied for this bubble condition varies from. to.8. The initial pressure in the nucleus bubble is set to be in equilibrium with the ambient pressure for the corresponding cavitation number σ defined as p p σ = v. ( U ρ A representative bubble splitting behavior from the simulations is shown in Figure for σ =.5. The figure shows the final stage of the simulation when the elongated bubble splits into two sub-bubbles which then develop reentrant jets. In general, bubbles do not always split. The same size nuclei at a slightly higher σ =.54 shows no splitting. Instead a single jet from the downstream end of the elongated bubble forms. The bubble size vs. time can be compared between the results of the SAP spherical model [4.5] and DYNAFS by using a volume equivalent bubble radius for the non-spherical bubble. The equivalent radius of the bubble as it flows downstream is shown in Figure 3 for two selected cavitation numbers. For the cavitation number σ =.5, the bubble experiences splitting shortly after reaching its maximum size, while a single reentrant jet forms for σ =.54 at the downstream end of the bubble. In both cases, the SAP spherical model predicts collapse and rebound sequences with decreasing maximum radii. The general trend of decreasing maximum bubble size as the cavitation number increases can be observed in both of the predictions by SAP spherical model and by DYNAFS. Equivalent Radius Equivalent Radius Re =.88e6, R o =5µm, σ =.5 Re =.88e6, R o =5µm, σ =.54 single jet in DynaFS DynaFS SAP DynaFS SAP Figure 3. Equivalent radius of the bubble as it flows downstream for cavitation numbers,.5 and.54. Figure 4, shows the maximum radius the bubble reached for a given σ, as well as its radius at the moment of splitting and the radii of the sub-bubbles generated. The size of the larger sub-bubble after splitting is only slightly smaller than the size of the bubble just before splitting. This difference between the two sizes is even smaller as the cavitation number increases and approaches the inception number. When the cavitation number is very close to inception no splitting occurs.

3 Figure 5 shows the relative sizes of the bubbles before and after splitting. We can see that the ratio of the maximum equivalent radii to the radius at splitting increases from.6 to.5 as the cavitation number increases from. to.5. The ratio of the larger sub-bubble generated to the equivalent radius of the bubble before splitting increases from.95 to.99, while the relative size of the smaller sub-bubble decreases from.55 to.. The curves of these ratios flatten as the cavitation number decreases. R EQ Equivalent Radius for Re=.88E6, R o =5µm Max. DynaFS, 8 seg DynaFS before splitting no splitting Figure 4. Equivalent radii of the bubble at the moment of splitting and radii of the sub-bubbles generated. minimum pressures during the initial sudden growth phase. Similarly, the collapse peak-to-peak pressure and the jet peakto-peak pressure are defined respectively for the collapse/rebound phase and for the jet development phase. 5 5 Field Point Pressure at (r,z = (.597,. m DynaFS SAP Spherical Model st peak to peak Figure 6. Pressure signals at a field point predicted for σ=.5. The peak-to-peak pressures are defined as the maximum pressure minus the minimum pressure within the indicated portion of the pressure signal. 5 Jet peak to peak collapse peak to peak Re =.88e6, Ro = 5e-6 m R aft,r aft,r max Maximum Before splitting After splitting # After splitting # 4 3 Field Point Pressure for Re=.88E6, R o =5µm Peak-to-peak pressure at (r,z=(.597,. m Split + Jets Single Jet Collapse/ Rebounds (Jets cannot be detected with 3 segments Sigma Figure 5. The ratios of the equivalent radii of the bubble and the sub-bubbles, relative to the equivalent radius just before the splitting. During each simulation, the pressure at a field point was computed and recorded. The field point was located.53 m away from the vortex center in the radial direction and at the same axial location of the tip vortex. As shown in Figure 6, the field point pressure predicted by SAP spherical model usually has multiple peaks with the first peak near the explosive growth of the bubble and the largest peak at the first collapse/rebound of the bubble. DYNAFS predicts similar behavior during the initial growth but develops much higher pressure peaks when a reentrant jet is formed. The first peak-to-peak pressure is defined as the difference between the maximum and the - SAP (st SAP (collapse DynaFS (st 8 seg DynaFS (jet 8 seg Figure 7. Comparison of the pressure peaks predicted by the SAP spherical model and by DYNAFS. The peak-to-peak values of the pressure at the field point for a range of cavitation numbers are compared in Figure 7. The first peak-to-peak pressure signals predicted by the two methods agree very well with each other for the range of cavitation number studied. However, the maximum jet peak-to-peak pressures from the DYNAFS prediction are about two orders of magnitude higher than the maximum collapse peak-to-peak pressures predicted by the SAP spherical model. This is due to the high pressures created at the development of the jet in the 3

4 DYNAFS predictions. Overall the flattening trend toward the lower cavitation numbers is observed for all pressure peaks. Axial Location of C.G Figure 8. Axial locations of the center of the bubble and the formed sub-bubbles for three cavitation number cases;.,.4, and.5 (from left to right The trajectories of the center of the initial bubble and the resulting sub-bubbles are shown in Figure 8. In all three cavitation number cases shown in the figure, the bubble nucleus slows down as it flows downstream and it almost stops traveling downstream at the moment of splitting. Just after the splitting, the center of the larger sub-bubble is located slightly upstream of the center of the pre-split bubble whereas the center of the smaller sub-bubble is located much downstream of the pre-split bubble. The centers of the sub-bubbles for the three cases shown in Figure 8 can be described in terms of the equivalent radius of each case. That is, -.7 and 4.9 for σ =., -.5 and 4.9 for σ =.4, and -.5 and 3.87 for σ =.5. maximum bubble size is found to be very similar for the two Reynolds numbers. The bubbles do not split for cavitation numbers just below the inception, however, the relative magnitude of the bubble size is smaller for the lower Reynolds number cases. Normalized Equivalent Radius (R EQ /a c.5.5 Equivalent Radius, Re=.44x 6 vs..88x 6,R o =5µm Re =.44x 6 -Cp min =.7 a c =.339 m at Cp min Max. SAP Max. DynaFS Re =.88x 6 -Cp min =.73 a c =.558 m at Cp min no splitting Normalized (σ/-cp min Figure. Comparison of the maximum equivalent radii for the two Reynolds numbers with both axes normalized..5 Equivalent Radius, Re=.44x 6 vs..88x 6,R o =5µm Equivalent Radius, Re=.44x 6 vs..88x 6,R o =5µm.75 Max. SAP Max. DynaFS Re=.88x 6.5 Re=.44x 6 no splitting R EQ.5 no splitting Normalized Equivalent Radius (R EQ /a c.5 Re =.44x 6 -Cp min =.7 a c =.339 m at Cp min Re =.88x 6 -Cp min =.73 a c =.558 m at Cp min Max. DynaFS DynaFS before splitting Normalized (σ/-cp min Figure 9. Comparison of the maximum equivalent radii for two Reynolds numbers,.44x 6 and.88x 6. REYNOLDS NUMBER EFFECT In order to observe the Reynolds number effect on the bubble size and the pressure peaks, the results in the previous section are compared with the corresponding results from the simulations using the tip vortex flow field at Reynolds number.44x 6. The maximum equivalent radii of the bubbles for the two Reynolds numbers are compared in Figure 9. The cavitation inception based on the bubble size occurs at σ =.57 for Re=.88x 6, and at σ =. for Re=.44x 6. The trend of the Figure. Comparison of the equivalent radii at the splitting for two for the two Reynolds numbers with both axes normalized. Because the cavitation inception is known to be scaled for different Reynolds numbers by the minimum pressure in the tip vortex core, the abscissa of Figure 9 can be normalized by using the minimum pressures, C p,min =.7 and.73, respectively for Re=.44x 6 and.88x 6. Also from our previous results [4] the bubble maximum size is scaled by the core radius of the tip vortex. Thus, the ordinate of Figure 9 can be normalized by the corresponding core radius. The tip vortex core radius at the axial location of minimum pressure is found from the RANS equation solutions as a c =.339 m and.558 m respectively for Re=.44x 6 and.88x 6. After normalizing 4

5 both axes, the cavitation inception curves for the two Reynolds numbers nicely collapse near σ / C p,min =.94, and the equivalent radius curves from the two Reynolds number cases also collapse nicely as shown in Figure. Based on the normalized bubble size, the bubble splitting happens only when Req / a c >.55~.6. This threshold is based on the 5 µm bubble nucleus here, but this value does not vary much for other nucleus sizes as will be discussed later. The equivalent radii just before and after the bubble splitting are compared for the two Reynolds numbers in Figure with the normalization described above. A strong resemblance between the two Reynolds numbers is observed in the bubble and sub-bubble sizes. Notice that the equivalent radii curves also collapse nicely. In Figure, the ratios of the equivalent radii are compared for the two Reynolds numbers using normalized cavitation numbers. The ratios look very similar for the two Reynolds numbers. The radius ratios flatten as the cavitation number decreases much below the cavitation inception number. R aft,r aft,r max Ratio of Radii, Re =.44x 6 vs..88x 6,Ro=5e-6m Re=.44x 6 Re=.88x 6 Maximum Before splitting After splitting # After splitting # Normalized (σ/-cp min Figure. Comparison of the equivalent radii ratios for the two Reynolds numbers with normalized cavitation number on the abscissa. The peak magnitudes of the pressure signals are compared in Figure 3. The pressure is computed at a field point.53 m away from the tip of the hydrofoil for the case of Re=.88x 6, but at a field point.8 m away for the smaller hydrofoil case of Re=.44 x 6. The normalization of the cavitation number also works for the pressure signal of the two Reynolds numbers. The highest pressure peak predicted by the SAP spherical model for Re=.44 x 6 has an up and down near the flat region of the lower cavitation numbers, which is due to the fact that the highest pressure is observed at the subsequent collapses/rebounds rather than at the usual first one. The magnitude of the pressure does not depend on the Reynolds numbers as expected, but depends on how far the cavitation number deviates from the cavitation inception number. The axial location of the bubble centers (such as shown earlier in Figure 8 is studied also for the Re=.44x 6. From this numerical study, the observed centers of the bubble and the sub-bubbles can be compared for the two Reynolds numbers as shown in Figure 4. The dimensional bubble center values of the two Reynolds number cases seem to follow roughly the length scale ratio,, of the two cases. It can be concluded that the upstream sub-bubble center just after the splitting approaches -.8 R eq while the downstream sub-bubble center approaches about 4. R eq as the cavitation number decreases for both Reynolds numbers Field Point Pressure with R o =5µm, Re=.44x 6 and.88x 6 Peak-to-peak pressure at (r,z=(.597,. m SAP (st min/max SAP (max. at collapse DynaFS (st min/max DynaFS (jet min/max Re=.44x 6 Re=.88x Normalized (σ/-cp min Figure 3. Comparison of the pressure signals for the two Reynolds numbers with normalized cavitation number on the abscissa. CG /R eq,cg /R eq 6 4 Center of Gravity of Sub-bubbles, R o =5µm Re=.44x 6 Re=.88x 6 CGb CGa CGa CG/Req CG/Req σ Figure 4. Axial location of the bubble center and the subbubble centers at the splitting for two Reynolds numbers,.44x 6 and.88x 6. Ro=5µm, CGb: center just before the splitting, CGa and CGa: centers of the sub-bubbles just after the splitting, Req: equivalent radius just before the splitting. BUBBLE NUCLEUS SIZE EFFECT The effect of the initial nucleus size is studied by comparing results from the simulations with four initial sizes:,, 5, and µm, in the flow field for Re=.88x 6. The radii just before and after the splitting as well as the maximum radii from the SAP spherical model and the DYNAFS are compared in Figure 5. The equivalent radii of the smaller CG b,cg a,cg a 5

6 initial bubble sizes seem to lie on straight lines extended from the lower cavitation number portion of the radius curves for the larger initial bubble sizes. The threshold cavitation number between bubble splitting and single jet formation decreases as the nucleus size decreases from µm to 5 µm. The threshold in terms of bubble size is R eq.8 m, which can be normalized as Req / ac.5. This value of the threshold is not very different from the value Req / ac.55~.6 found earlier for the two Reynolds number cases. Note that single jet behaviors cannot be observed for the smaller initial bubble radii of and µm. In Figure 6, the ratios of the equivalent radii are compared for the four initial bubble radii. The flattening trend prevails for the lower cavitation numbers regardless of the initial bubble sizes. As the cavitation number decreases, the ratio for the maximum equivalent radius approaches.6, the ratio for the larger sub-bubble just after the split approaches.95, that for the smaller sub-bubble approaches.55 regardless of the initial bubble size...5 Equivalent Radius, Re =.88x 6,R o =,, 5, µm Ro = µm Max. SAP Max. DynaFS DynaFS before splitting cavitation numbers. This might be due to the fact that the jet developed from µm initial bubble is much longer (through a long tubular part of the sub-bubble than the other jets found in larger size bubble simulations. Rmax, Rbef, Raft, Raft Ratio of Radii, Re =.88x 6, Ro =,, 5, µm Ro = µm Ro = µm Ro = 5 µm Ro = µm Maximum Before splitting After splitting # After splitting # Sigma Figure 6. Comparison of the ratios of the equivalent radii for four initial bubble sizes,,, 5, and µm. Field Point Pressure at Re=.88x 6, with Ro=,, 5, and µm Peak-to-peak pressure at (r,z=(.597,. m 6 µm SAP (st min/max SAP (max. at collapse DynaFS (st min/max DynaFS (jet min/max µm 4 5 µm R EQ. Ro = µm Ro = 5 µm no splitting (5 µm µm.5 no splitting ( µm - Ro = µm Figure 5. Comparison of the equivalent radii for four initial bubble sizes,,, 5, and µm. The peak-to-peak values of the pressure signal peaks are compared in Figure 7. The first peak pressures predicted by the SAP spherical model and by DYNAFS agree throughout the studied cavitation numbers. Moreover these pressure data form a common curve regardless of the initial bubble sizes. This curve of the first peak pressure flattens out as the cavitation number decreases. The maximum peak pressure predicted by the SAP spherical model forms a hump just below the cavitation inception but becomes closer to the first peak pressure as the initial bubble size decreases and the cavitation number decreases. This is due to the fact that the first peak becomes the largest peak as the cavitation number and the initial bubble size decrease. The maximum pressure peaks predicted by DYNAFS are always observed at the development of the jet that follows the bubble splitting. The maximum pressure peaks predicted with 5 and µm show more or less similar behavior with a small hump just below the cavitation inception number. This curve continues smoothly to the data obtained from and µm initial bubbles, curving up slightly toward the lower Figure 7. Comparison of pressure signals for four initial bubble sizes:,, 5, and µm. The axial location of the bubble center obtained from several computations with different initial bubble nucleus sizes is shown in Figure 8. It is interesting to observe that the subbubble centers normalized by the equivalent radius form common curves throughout the cavitation number ranges studied. The common curve for the downstream sub-bubbles is observed to be more or less flat near 4.4 R eq for the cavitation numbers between.7 and.. The normalized axial location of the upstream sub-bubble varies from approximately -. in the low cavitation numbers to. at the cavitation number where the splitting does not happen. 6

7 6 Center of Gravity of Sub-bubbles, Re=.88x 6 R o =µm. Bubble splitting criterion (for low cavitation numbers: When the bubble reaches it maximum radius and the radius decreases to.95 of the maximum, the bubble splits. CG /R eq,cg /R eq 4 R o =µm R o =5µm CGb CGa CGa CG/Req CG/Req R o =µm.5..5 CG b,cg a,cg a Equivalent Radius Re =.88e6, R o =5µm, σ = σ Figure 8. Axial location of the bubble center and the subbubble centers at the splitting for four initial bubble sizes,,, 5, and µm. Re=.88x6, CGb: center just before the splitting, CGa and CGa: centers of the sub-bubbles just after the splitting, Req: equivalent radius just before the splitting. BUBBLE SPLITTING MODEL The axisymmetric simulations using DYNAFS have the advantage of predicting non-spherical bubble behaviors, such as elongation, splitting, and reentrant jets, and the resulting pressure signals. However, these simulations are more costly than simulations using the SAP spherical model. In order to perform massive simulations with multiple bubble nuclei, a simpler model is required so that the bubble splitting can be handled within the spherical model. In the upgraded spherical model, once a bubble reaches a condition that it should split, the bubble is replaced with two sub-bubbles and the simulation is continued with the new sub-bubbles. The splitting model will be more useful for cavitation numbers much lower than the inception because the bubble usually does not split very near the cavitation inception. The splitting prediction model should be able to answer the next two questions. Criteria of splitting: When does a bubble split? Initial condition of the sub-bubbles: What is the size of each sub-bubble? What is the pressure inside each subbubble? From the studies described previously [3], a bubble splits when its elongation reaches a value between.8 and 4.8 depending on the cavitation number and the initial bubble size. Since the elongation cannot be calculated within the spherical model, it cannot be used in the spherical model to determine if a bubble should split. The ratio of the maximum equivalent bubble radius relative to the equivalent radius just before the splitting can be used as a criterion. This ratio is flat at.5 for lower cavitation numbers and increases up to.5 as the cavitation number approaches the inception. In other words, the following criterion can be used to determine when a bubble should split Figure 9. Equivalent radius of the bubble as a function of time predicted with 5 µm nucleus at σ=.. Equivalent Radius Re =.88e6, R o =µm, σ =.54 DynaFS SAP Figure. Equivalent radius of the bubble as a function of time predicted with a µm nucleus at σ=.54. Once the splitting is detected from the above criteria, the size of sub-bubbles can be determined also from the ratios shown in Figure 6. The equivalent radii of the larger and the smaller sub-bubbles are, respectively,.95 and.55 of the equivalent radius before the splitting. The gas pressure inside the sub-bubbles can be modeled with the introduction of a pressure reduction factor α such that the pressure after the splitting can be written as pg, after = α pg, before. ( In an ideal case of lossless splitting, the factor α is.. The radial velocity of the sub-bubbles can be observed indirectly from equivalent radius curves such as shown in Figure 3. A magnified view near the splitting is shown in Figure 9 for the case of σ =.. The slope of the equivalent radius in the figure can be interpreted as the average radial velocity of the bubble. Just after the split, the larger sub-bubble experiences a very slow radial velocity (nearly zero slope while the smaller sub-bubble has a slightly faster radial velocity (steeper slope than before the split. This trend is also found in most of the bubble behaviors simulated with 5 and µm bubble nuclei at relatively higher cavitation numbers. However, in simulations 7

8 with and µm bubble nuclei at relatively lower cavitation numbers, the slope just after the split tends to be zero for both sub-bubbles as shown in Figure. Therefore, zero radial velocity as the initial condition of the sub-bubbles seems to be appropriate for the lower cavitation number simulations. The initial location of the sub-bubbles depends weakly on the Reynolds numbers (Figure 4 and more strongly on the cavitation numbers (Figure 8. The location can be read directly from these two figures. For cavitation numbers in the mid-low range (.7 σ., the locations of the subbubbles relative to the bubble center just before the splitting are fairly constant and are -. and 4.4, respectively for the upstream and the downstream ones. The initial condition of the new sub-bubbles can be summarized as follows. Initial condition of the sub-bubbles (for low cavitation numbers: The radii of the larger and smaller sub-bubbles are respectively,.95 and.55 of the radius just before splitting. The initial gas pressure is determined from equation (, and the initial radial velocity is zero. The initial locations of the subbubbles are determined from Figure 4 or Figure 8. If a simpler approach is desired, it can be assumed that the initial location of the larger sub-bubbles is. radius upstream of the pre-split bubble center and that of the smaller one is 4.4 radii downstream of the pre-split bubble center. CONCLUSION Bubble behaviors in two tip vortex flow fields of Reynolds number.44x 6 and.88 x 6 were studied by using the SAP spherical model and DYNAFS, a non-spherical free surface flow solver. The equivalent radius of the non-spherical bubble and the radius predicted by the spherical model are compared in a range of cavitation numbers varying from.5 to.8. In addition, the peak values of the pressure signals detected from the bubble collapse/rebound and from the formation of a jet or jets are compared. The effect of the Reynolds number is found to be mainly a shifting of the cavitation inception number both in the bubble radii and in the pressure peaks. With cavitation numbers normalized by the minimum pressure in the tip vortex core, the curves from the two Reynolds numbers collapse very close to each other nicely as expected. The effect of initial bubble nucleus size is found to be such that, if observed from higher to lower cavitation numbers, a smaller initial bubble nucleus brings an abrupt transition from a non-cavitating status to the cavitating status represented by the common curves of the radii (or of the pressure peaks much later than a larger bubble nucleus case. These common curves seem to be independent of the nucleus size within the range of to µm for Reynolds number.88 million. The observation of the ratios of the equivalent bubble maximum radii, before and after splitting, leads to a simple model of the bubble splitting for low cavitation numbers. This suggested model is based on the ratios of the equivalent radii of the bubble and sub-bubbles. Studies to expand and apply the model are currently underway. ACKNOWLEDGMENTS This work was conducted at DYNAFLOW, INC. ( The support of several colleagues including Dr. Chao-Tsung Hsiao is appreciated. The work has been supported by the Office of Naval Research under contract No. N4-99-C-369 monitored by Dr. Ki-Han Kim. This support is greatly appreciated. DYNAFLOW s RANS code DF_UNCLE is based on the UNCLE code originally developed by Mississippi State University and graciously provided to DYNAFLOW. REFERENCES [] Choi, J.-K., Chahine, G. L., Non-spherical bubble behavior in vortex flow fields, Computational Mechanics, (in print 3. [] Choi, J.-K., Chahine, G. L., Noise due to extreme bubble deformation near inception of tip vortex cavitation, Proc. FEDSM 3, International Symposium on Cavitation Inception, 4th ASME/JSME Joint Fluids Engineering Conference, Honolulu, Hawaii, July 6-, 3. [3] Choi, J.-K., Chahine, G. L., A numerical study on the bubble noise and the tip vortex cavitation inception, 8 th International Conference on Numerical Ship Hydrodynamics, Busan, Korea, September -5, 3. [4] Hsiao, C.-T., Chahine, G. L., Prediction of vortex cavitation inception using coupled spherical and nonspherical models, Proc. 4th Symposium on Naval Hydrodynamics, ONR, Fukuoka, Japan, July 8-3,. [5] Hsiao, C.-T., Chahine G. L., Liu, H., Scaling effects on prediction of cavitation inception in a line vortex flow, Journal of Fluids Engineering, Vol.5, 3, pp [6] Chahine, Georges L., Fluid Vortices, S. Green (ed., Kluwer Academic, Chapter 8, 995. [7] Chahine, G. L., Duraiswami, R., Kalumuck, K. H., Boundary element method for calculating -D and 3-D underwater explosion bubble loading on nearby structures including fluid structure interaction effects, Technical Report NSWC-DD/TR-93/46, Dynaflow, Inc., 996. [8] Chahine, G. L., Experimental and asymptotic study of nonspherical bubble collapse, Applied Scientific Research, Vol.38, 98, pp [9] Kucherenko, V. V., Shamko, V. V., 986, Dynamics of electric-explosion cavities between two solid parallel walls, Journal of Applied Mechanics and Technical Physics, Vol.7, 986, pp.-5. [] Ishida, H., Nuntadusit, C., Kimoto, H., Nakagawa, T., Yamamoto, T., Cavitation bubble behavior near solid boundaries, Proc. CAV, 4th International Symposium on Cavitation, session A3, California Inst. of Tech., Pasadena, CA.. [] Chahine, G. L., Sarkar, K., Duraiswami, R., Strong bubble/flow interaction and cavitation inception, Technical Report 943-ONR, DYNAFLOW, INC.,