The boundary mixture model is derived to predict the performance of the microbubble drag reduction technique for a flat plate. The flat plate with a porous material microbubble injecting system and resistance-measuring system are set up to measure the frictional resistance of the flat plate without and with injected microbubbles. The tests are conducted in a water tunnel and a towing tank. The test results show that the boundary mixture model predicts the drag reduction well for the flat plate when testing with injected microbubbles in the water tunnel. However, the boundary mixture model overestimates the drag reduction effect for the flat plate tested in the towing tank. The possible mechanism for the overestimation of drag reduction effect in the towing tank may be due to the different behaviors of microbubbles in the velocity gradient of boundary layer.
1. Introduction
The microbubble drag reduction technique has been one of the major research topics both in fluid mechanics and in ship hydrodynamics over the last two decades [1]. Since Madavan et al. [2] achieved an impressive 80% drag reduction in a water tunnel test, many researchers devoted themselves to the study of the microbubble drag reduction technique. The major research directions of the microbubble drag reduction technique include (1) developing the microbubble injection methods and to validate the drag reduction effect in a water tunnel, such as by using the electrolysis method [1], the porous method [2, 3], or the WAIP (winged air induction pipe) method [4], (2) experimentally studying the mechanisms and the key parameters of the microbubble drag reduction technique, such as the bubble size effect [5], the bubble concentration effect [6], the gas type effect [7], or the effect of microbubbles on the structure of turbulence [8], (3) modeling the mechanism of the microbubble drag reduction technique by numerical methods [9] or by DNS simulation [10], and (4) applying the microbubble drag reduction technique to the ship model and conducting the ship model resistance test to study the drag reduction effect in a towing tank [11–13].
Some important research results have emerged for the microbubble drag reduction technique over the last two decades. First, the drag reduction effect of the microbubble in the flat plate was confirmed and an 80% maximum drag reduction effect was also validated in water tunnel tests [2, 3, 14, 15]. Second, it was confirmed that the effect of bubble size (ranging from 0.05 mm to 5 mm) on drag reduction could be neglected [5, 16]. Third, the most important parameter for microbubble drag reduction is the volume fraction of the injected air in the water boundary layer. Skudarnov and Lin [17] also validated the density effect on the drag reduction effect of microbubble drag reduction by using a 2D RANS simulation. But the prediction model for the microbubble drag reduction technique is still not well developed.
In this study, a simple boundary layer mixture model is derived to predict the drag reduction effect of the microbubble drag reduction technique for a flat plate. The flat plate resistance measuring system with a porous material microbubble injecting system is set up to measure the total resistances of the plate without and with injecting the microbubbles. The tests are conducted in a water tunnel and a towing tank to verify the drag reduction effect predicted by the boundary layer mixture model.
2. Boundary Mixture Model
Figure 1 is a schematic drawing of a microbubble injecting system for a plate on the top. The injected microbubbles are assumed to be distributed uniformly across the boundary layer. The air volume fraction Cv is defined as the ratio of the injected air flow rate divided by the summation of the air flow rate and the water flow rate within the boundary layer,
Cv=QaQa+Qw,
where Qa is the injected air flow rate, and Qw is the water flow rate within the boundary layer of the plate. Based on the turbulent boundary theory [18], the water flow rate within the boundary layer of the plate can be calculated by
Qw=U0(δ-δ*)b,
where b is the width of the plate, U0 is the inflow velocity, and δ is the boundary layer thickness, which is defined as the distance from the wall where the velocity is 0.99U0. A seventh power velocity distribution [18] is assumed for the velocity distribution across the boundary layer
UU0=(yδ)1/7.
And the displacement thickness δ* is defined as
δ*=∫0δ(1-UU0)dy.
The Schlichting boundary thickness formula [18] is used to estimate the thickness of the boundary layer
δ=0.37xRex0.2,
where Rex=U0x/ν and x is the distance from the origin of the plate.
Schematic drawing of a microbubble injecting system (plate on top).
Using (4) and (5) to solve (2), the water flow rate can be calculated by
Qw=0.3238U0bxRex-0.2.
Then the air volume fraction Cv in (1) can subsequently be calculated with measured injected air flow rate and the water flow rate which was estimated by using (6).
For a flat plate without injected bubbles, the frictional resistance Df of a flat plate with length l and width b can be derived as [18]
Df=0.036ρwU02l2bRel-0.2,
where Rel=U0l/ν and l is the total length of the plate, and ρw is the density of the water.
The resistance coefficient Cf is defined as
Cf=Df0.5ρwU02lb=0.074lRel-0.2.
For the water-bubble mixture of the boundary layer with injected bubbles, the mixture density ρb can be calculated by the linear combination of the density of air and the density of water according to the air volume fraction Cv, and is given by
ρb=ρaQa+ρwQwQa+Qw=ρaCv+ρw(1-Cv),
where ρa is the density of the injected air.
The dynamic viscosity of the water-bubble mixture can also be calculated by using the same approach as
μb=μaQa+μwQwQa+Qw=μaCv+μw(1-Cv).
The kinematical viscosity of the water-bubble mixture νb is then defined as
νb=μbρb.
The frictional resistance of a flat plate with a water-bubble mixture boundary layer Dfb can be calculated by using the same approach [18] as
Dfb=0.036ρbU02l2bRebl-0.2,
where the Reynolds number of the water-bubble mixture Rebl is defined as
Rebl=U0lνb.
The ratio of the frictional resistance of water-bubble mixture to the frictional resistance of the water is then expressed as
DfbDf=ρbRebl-0.2ρwRel-0.2.
Equation (14) predicts the frictional resistance of water-bubble mixture of the plate from the frictional resistance of flat plate in pure water.
If the frictional resistances Dfb and Df are normalized by the dynamic pressure of the inflow (1/2)ρwU02 and the surface area of the flat plate b*l, then the same result will be had for the nondimensional resistance coefficients as will be
CfbCf=ρbRebl-0.2ρwRel-0.2,
where Cf=Df/(1/2)ρwU02lb and Cfb=Dfb/(1/2)ρwU02lb.
Equation (15) predicts the nondimensional frictional resistance of water-bubble mixture of the plate from the nondimensional frictional resistance of flat plate in water. The drag reduction ratio DR predicted by the boundary layer mixture model can be calculated by
DR=1-CfbCf.
Figure 2 shows the ratio of nondimensional resistance coefficients of the bubble-water mixture to the pure water with the parameter of the air volume fraction and the drag reduction ratio. The effect of the density of the mixture and Reynolds number on the microbubble drag reduction technique is shown in Table 1. The effect of the Reynolds number is very small when compared with the effect of density of the mixture. The density of the bubble mixture becomes the key parameter for the microbubble drag reduction technique. The ratio of the frictional resistance of the water-bubble mixture boundary layer to the water boundary layer is almost directly proportional to the density ratio.
Density and Reynolds number effect on the microbubble drag reduction.
Cv
ρb/ρw
Rebl-0.2/Rel-0.2
Cfb/Cf
DR
0.000
1.000
1.000
1.000
0.000
0.100
0.900
1.000
0.900
0.100
0.200
0.800
1.001
0.801
0.199
0.300
0.700
1.001
0.701
0.299
0.400
0.600
1.002
0.602
0.398
0.500
0.501
1.003
0.502
0.498
0.600
0.401
1.004
0.402
0.598
0.700
0.301
1.006
0.303
0.797
0.800
0.201
1.011
0.203
0.897
0.900
0.101
1.024
0.103
0.997
0.990
0.011
1.174
0.013
0.987
The drag reduction effect predicted by the boundary layer mixture model.
3. Experimental Apparatus and Test Procedures
Figure 3 shows a flat plate resistance measuring system. The length of the flat plate is 700 mm and the width is 190 mm. The flat plate is suspended by four steel strips at the four corners. The force gage is located at the center of the plate to measure the resistance of the suspended flat plate. Figure 4 shows a porous plate microbubble injecting system, which is designed to inject the microbubbles through the porous plate. The length of the porous plate is 120 mm and the width is 190 mm. The pore sizes of porous plate used in the water tunnel test are 10 μm and 100 μm. The microbubble injecting system is placed in front of the flat plate resistance measuring system, and they are combined together as a whole unit. Figure 5 shows the whole system installed in the water tunnel. The water tunnel is the K23 cavitation tunnel from K&R Germany.
The flat plate resistance measuring system.
Air chamber and porous plate (the yellow one).
The whole microbubble resistance measuring system in water tunnel.
The resistances of the flat plate without injecting microbubbles are first measured for different velocities varying from 2.45 m/s to 4.2 m/s. The resistances of the flat plate while injecting microbubbles are then measured for different air flow rates for each velocity. The air flow rates range from 12 L/min to 75 L/min.
The same system is also installed at the towing carriage to conduct the resistance test in the towing tank. The dimensions of the towing tank are 130 m in length, 8 m in width, and 4 m in depth. The maximum speed of the towing carriage is 5 m/s. Figure 6 shows the flat plate resistance measuring system with the towing carriage in the towing tank. The test procedure and conditions in the towing tank are the same as in the water tunnel.
The whole microbubble resistance measuring system in towing tank.
4. Test Results and Validation of the Boundary Layer Mixture Model
Figure 7 shows the measured resistance coefficients of the flat plate without injecting microbubbles tested in the water tunnel and in the towing tank. The measured resistance coefficients are about the same, both in the water tunnel and in the towing tank. The resistance coefficient calculated by the formula (8) is also shown in Figure 7. The measured resistance coefficients of the flat plate both from the water tunnel and the towing tank are a little higher than the resistance coefficients predicted by the formula (8). It is reasonably common to achieve a higher resistance in the test environment than when using the formula because of the imperfection of the flat plate test setup and the three-dimensional effects.
Resistance of flat plate without injected microbubbles.
Figure 8 shows the resistance coefficient ratio of the flat plate with the injected microbubbles by a 10 μm porous plate and without the injected microbubbles tested in the water tunnel. The measured resistance coefficient ratio is in good agreement with the drag reduction line predicted by the boundary mixture model. Figure 9 shows the resistance coefficient ratio of the flat plate with the injected microbubbles by a 100 μm porous plate and without the injected microbubbles tested in the water tunnel. The measured resistance coefficient ratio of the 100 μm porous plate is also in good agreement with the frictional resistance line predicted by the boundary mixture model.
The drag reduction ratio of the flat plate with a 10 μm porous plate injecting microbubbles in the water tunnel.
Drag reduction ratio of the flat plate with a 100 μm porous plate injecting microbubbles in the water tunnel.
Figure 10 shows the resistance coefficient ratio of the flat plate with the injected microbubbles by a 10 μm porous plate and without the injected microbubbles tested in the towing tank. The drag reduction effect of the 10 μm porous plate in towing tank is much smaller then that predicted by the boundary mixture model. Figure 11 shows the resistance coefficient ratio of the flat plate with the injected microbubbles by a 100 μm porous plate and without the injected microbubbles tested in the towing tank. The drag reduction effect of the 100 μm porous plate in towing tank is also much smaller then that predicted by the boundary mixture model. In addition, the drag reduction effect is decreased at a higher air flow rate for velocities 2.45 m/s and 2.9 m/s. For example, at a velocity 2.45 m/s, the resistance ratio is 0.77 when the air volume fraction is 0.65, but the resistance ratio is 0.85 when the air volume fraction is 0.71. The drag reduction effect is decreased from 23% to 15% when the air volume fraction is increased from 0.65 to 0.71. This implies that there exits an optimal air flow rate when using the microbubble drag reduction technique in the towing tank for each velocity. This trend is different from the air flow rate increasing monotonically with increasing the drag reduction effect in the water tunnel.
Drag reduction ratio of the flat plate with a 10 μm porous plate injecting microbubbles in the towing tank.
Drag reduction ratio of the flat plate with a 100 μm porous plate injecting microbubbles in the towing tank.
5. Discussions
The test results show some differences for the microbubble drag reduction technique tested in the water tunnel and in the towing tank. The first is that the maximum drag reduction effect of the microbubbles in the water tunnel is about 80%, and the drag reduction effect is in good agreement with the value predicted by the boundary mixture model. However, in this study, the maximum drag reduction effect of microbubble in the towing tank is only about 30%. Takahashi et al. [19] only achieved a 22% drag reduction for a 50-meter flat plate ship in the towing tank. The second is that the drag reduction effect increases monotonically with an increasing the air flow rate in the water tunnel. However, an optimal air flow rate exists in the towing tank. The drag reduction effect will be decreased when the air flow rate is above a critical air flow rate. The third is that the drag reduction effect is only a function of the air volume fraction Cvand is independent of velocity in the water tunnel. Madavan et al. [2] showed the same trend in which the drag reduction is only a function of nondimensional air flow rate Cq=Qa/(U0A)in water tunnel. However, in the towing tank, the drag reduction effect is slightly dependent on velocity as shown in Figures 10 and 11.
The mechanism that creates the difference in the drag reduction effect of the flat plate when it is tested between the water tunnel and the towing tank plays a crucial role for the application of microbubble drag reduction technique in ships. The most likely mechanism may be the different behaviors of microbubbles between in the water tunnel and in the towing tank. A transparent acrylic plate is used instead of the steel plate of the resistance measuring system to observe the microbubble behaviors in the water tunnel and in the towing tank. Figure 12 shows the distribution of microbubbles in the water tunnel. The microbubbles are distributed uniformly across the whole plate, and the individual bubbles can be seen clearly. Figure 13 shows the bubble distribution across the flat plate in the towing tank. The bubbles are deformed irregularly, and many bubbles are attached to the plate. The irregular black holes in Figure 13 are large deformed bubbles attached to the plate. There are more like an air packets for the white bubbles as shown in Figure 13.
Microbubble distribution across the flat plate in the water tunnel.
Microbubble distribution across the flat plate in the towing tank.
The different behaviors of bubbles in the water tunnel and in the towing tank may be due to the different velocity gradient in the boundary layer. The bubble will suffer a lift force when it remains in the velocity gradient. Figure 14 shows a schematic diagram for a bubble of radius R being lifted by a velocity gradient. The lift force can be calculated by
FL=ρwVCLU⃗×ϖ0,
where V is the volume of the bubble, ϖ0=∇×U⃗ is calculated at the bubble center, and CL is the lift coefficient. For an inviscid fluid satisfying the nondimensional shear rate Sr=2Rα/Um≪1, Auton et al. [20] obtained the lift force as
FL=23πρwR3αUm,
where α is the linear slope of the velocity gradient across the bubble and Um is the minimum velocity of the bubble encountered.
Bubble’s lift force produced by a velocity gradient.
The direction of the lift force is from the low velocity to the high velocity.
Besides the lift force, the bubble also has a buoyancy force in the upward direction. Figure 15 shows the velocity gradients of the boundary layer both in the water tunnel and in the towing tank. For the water tunnel, the flat plate is fixed at the top of the water tunnel wall, the velocity of the flat plate is zero, and the velocity outside the boundary of the flat plate is U0. Because the microbubble is injected into the boundary layer, it will suffer a velocity gradient and produce a lift force with a downward direction and an upward buoyancy force. The two forces could be balanced in some situations, so the bubble could stay in the buffer layer [21] and will not attach to the surface of the flat plate as shown in Figure 15(a).
Forces of a bubble suffered in the boundary layer of the water tunnel and the towing tank.
For the flat plate installed in the towing tank, the flat plate moves with the speed U0 of the towing carriage. The speed of the water outside the boundary layer of the plate is zero. The velocity gradient of the boundary layer in the towing tank is reversed when compared with the velocity gradient in the water tunnel. The direction of the bubble’s lift force produced by the velocity gradient as well as the direction of the bubble’s buoyancy is upward in the towing tank. The two forces will push the bubble toward the surface of the flat plate, and the bubbles will attach to the surface of the flat plate as shown in Figure 15(b). The bubbles will be accumulated at the surface of the plate. In addition, the coalescence of bubbles occurs at the surface of the plate and forms several large flat bubbles attached to the surface of the moving flat plate, which can be seen in Figure 13. The resistance of the flat plate is increased when a large deformed bubble is attached to its surface. Also, the drag reduction effect is decreased. Future study is required to examine the detailed process of the injected bubble to forming an irregular large bubble and attaching to the flat plate surface in the towing tank.
6. Conclusions
A boundary layer mixture model was derived to predict the drag reduction effect of microbubble drag reduction technique in the flat plate. Also, a flat plate total drag measurement with microbubble injecting system was designed to conduct the resistance test in the water tunnel as well as in the towing tank. From the prediction of the boundary mixture model, the test results, and the discussions, the following conclusions can be drawn.
The drag reduction effect predicted by the boundary mixture model is almost directly proportional to the density ratio of mixture and water.
The maximum drag reduction effect of the microbubbles in water tunnel is about 80%, and the frictional resistance coefficient is in good agreement with the value predicted by the boundary mixture model. The maximum drag reduction effect of the microbubbles in the towing tank is only about 30% and drag reduction is much smaller than that predicted by the boundary mixture model.
The drag reduction effect increases monotonically with increasing the air flow rate in the water tunnel. However, an optimal air flow rate exists for each velocity in the towing tank.
The different drag reduction effect in the water tunnel and in the towing tank may be due to the different bubble behaviors produced by the different velocity gradient. Future study is required to examine the details of the developing process of the injected bubble in the towing tank.
NomenclaturesA:
Area of the porous plate
b:
Breadth of the flat plate
Cf:
Frictional resistance coefficient of the flat plate in water
Cfb:
Frictional resistance coefficient of the flat plate in bubble-water mixture
CL:
Lift coefficient of bubble in velocity gradient
Cq:
Nondimensional air flow rate
Cv:
Air volume fraction
Df:
Frictional resistance of the flat plate in water
Dfb:
Frictional resistance of the flat plate in bubble-water mixture
DR:
Drag reduction ratio
FL:
Lift of a bubble in velocity gradient
l:
Length of the flat plate
R:
Radius of a bubble
Rel:
Reynolds number of flat plate with length l in water
Rebl:
Reynolds number of flat plate with length l in bubble-water mixture
Rex:
Reynolds number at any x coordinate of flat plate in water
Sr:
Nondimensional shear rate of a bubble
U:
Velocity component in x axes within the boundary layer of the flat plate
U0:
Uniform inflow velocity outside boundary layer of the flat plate
Um:
Minimum velocity of a bubble experienced within the diameter of the bubble
V:
Volume of a bubble
x:
Coordinate of horizon axes along the flat plate
y:
Coordinate of vertical axes perpendicular to the flat plate
α:
Slope of velocity within the range of bubble’s diameter
ρa:
Density of air
ρb:
Density of bubble-water mixture
ρw:
Density of water
μa:
Dynamic viscosity of air
μb:
Dynamic viscosity of bubble-water mixture
μw:
Dynamic viscosity of water
νa:
Kinematical viscosity of air
νb:
Kinematical viscosity of bubble-water mixture
νw:
Kinematical viscosity of water
δ:
Boundary layer thickness
δ*:
Boundary layer displacement thickness
ϖ0:
Vorticity of a bubble within a velocity gradient.
Acknowledgment
The authors are grateful for the funding support of the National Science Council of Taiwan (NSC92-2611-E-002-022).
McCormickM. E.BhattacharyyaR.Drag reduction of a submersible hull by electrolysis197385211162-s2.0-0015615228MadavanN. K.DeutschS.MerkleC. L.Reduction of turbulent friction by microbubbles1983TM 83-23The Pennsylvania State UniversitySandersW. C.CeccioS. L.IvyE. M.PerlinM.DowlingD. R.Microbubble drag reduction at high Reynolds numberProceedings of the 4th ASME/JSME Joint Fluids Engineering ConferenceJuly 2003Honolulu, Hawaii, USA6836952-s2.0-0346903931KumagaiI.NakamuraN.MuraiY.TasakaY.TakedaY.TakahashiY.A new power-saving device for air bubble generation: hydrofoil air pump for ship drag reductionProceedings of the International Conference on Ship Drag Reduction SMOOTH-SHIPSMay 2010Istanbul, TurkeyShenX.CeccioS. L.PerlinM.Influence of bubble size on micro-bubble drag reduction20064134154242-s2.0-3374789203810.1007/s00348-006-0169-yGuinM. M.KatoH.YamaguchiH.MaedaM.MiyanagaM.Reduction of skin friction by microbubbles and its relation with near-wall bubble concentration in a channel1996152412542-s2.0-0030354046FontaineA. A.DeutschS.The influence of the type of gas on the reduction of skin friction drag by microbubble injection1992132-31281362-s2.0-000856647010.1007/BF00218158KatoH.IwashinaT.MiyanagaM.YamaguchiH.Effect of microbubbles on the structure of turbulence in a turbulent boundary layer2000441551622-s2.0-23044520408MadavanN. K.MerkleC. L.DeutschS.Numerical investigation into the mechanisms of microbubble drag reduction198510733703772-s2.0-0022130249SugiyamaK.KawamuraT.TakagiS.MatsumotoY.Numerical simulations on drag reduction mechanism by microbubblesProceeding of the 3rd Symposium on Smart Control of Turbulence2002University of Tokyo, Japan129139LatorreR.Ship hull drag reduction using bottom air injection19972421611752-s2.0-0031080980KodamaY.TakahashiT.MakinoM.HoriT.UedaT.Practical application of microbubbles to ships—large scale model experiments and a new full scale experimentProceedings of the 6th International Symposium on Smart Control of TurbulenceMarch 2005Tokyo, JapanThillC.ToxopeusS.van WalreeF.Project Energy-saving air-Lubricated Ships (PELS)Proceedings of the 2nd International Symposium on Seawater Drag ReductionMay 2005Busan, KoreaKatoH.MiyanagaM.YamaguchiH.Frictional drag reduction by injecting bubbly water into a turbulent boundary layer and the effect of plate orientation1994190ASME185194TsaiJ. F.ChenC. C.Experimental study on the micro-bubble drag reduction effect in water tunnel and towing tankProceedings of the 3rd Asia-Pacific Workshop on Marine Hydrodynamics (APHydro '06)June 2006Shanghai, ChinaKawamuraT.KakugawaA.KodamaY.MoriguchiY.KatoH.Controlling the size of microbubbles for drag reductionProceeding on 3rd Symposium on Smart Control of Turbulence2002University of Tokyo, Japan129139SkudarnovP. V.LinC. X.Drag reduction by gas injection into turbulent boundary layer: density ratio effect20062734364442-s2.0-3364593392210.1016/j.ijheatfluidflow.2005.12.002SchlichtingH.19686thMcGraw-HillTakahashiT.KakugawaA.MakinoM.KodamaY.Experimental study on scale effect of drag reduction by microbubbles using very large flat plate ships20032391120AutonT. R.HuntJ. C. R.Prud'hommeM.Force exerted on a body in inviscid unsteady non-uniform rotational flow19881972412572-s2.0-0024142881Ortiz-VillafuerteJ.HassanY. A.Investigation of microbubble boundary layer using particle tracking velocimetry2006128350751910.1115/1.2174062