Reduction of skin friction by microbubbles and its relation with near-wall bubble concentration in a channel

Determination of the flow structure near the wall is essential for a clear insight into the phenomenon of skin friction reduction by microbubbles in a turbulent boundary layer. An important parameter, is the bubble concentration or void fraction in the wall region in drag-reducing conditions. The purpose of this paper is to show drag-reducing effects due to microbubbles in a water channel and, more importantly, to show the dependence of the drag-reduction values on the near-wall void fraction. A two-dimensional channel with an aspect ratio of 10 was specially built for this purpose with provisions for air injection through porous plates. Skin friction was directly measured by a miniature floating element transducer with a 5-mm circular sensing disk mounted flush on the top wall 67 channel-heights downstream of the injector. The wall friction in the presence of air bubbles was found to be reduced under the same bulk velocity when compared with the value without air. Detailed void fraction profiles across the channel were obtained by a sampling probe and a fiber-optic probe. Better collapse of the drag reduction data, independent of different profile shapes, was found when plotted against the near-wall void fraction than against a cross-sectional mean void fraction. While this dependence reconfirms that the phenomena are essentially inner-region dependent, the lack of influence of the bubble distribution patterns away from the wall implies lack of outer region influence.     

采用微气泡缩减高速船拖曳阻力(英文)

Ship hull form of the underwater area strongly influences the resistance of the ship.The major factor in ship resistance is skin friction resistance.Bulbous bows,polymer paint,water repellent paint(highly water-repellent wall),air injection,and specific roughness have been used by researchers as an attempt to obtain the resistance reduction and operation efficiency of ships.Micro-bubble injection is a promising technique for lowering frictional resistance.The injected air bubbles are supposed to somehow modify the energy inside the turbulent boundary layer and thereby lower the skin friction.The purpose of this study was to identify the effect of injected micro bubbles on a navy fast patrol boat(FPB) 57 m type model with the following main dimensions:L=2 450 mm,B=400 mm,and T=190 mm.The influence of the location of micro bubble injection and bubble velocity was also investigated.The ship model was pulled by an electric motor whose speed could be varied and adjusted.The ship model resistance was precisely measured by a load cell transducer.Comparison of ship resistance with and without micro-bubble injection was shown on a graph as a function of the drag coefficient and Froude number.It was shown that micro bubble injection behind the mid-ship is the best location to achieve the most effective drag reduction,and the drag reduction caused by the micro-bubbles can reach 6%-9%.     

Effect of microbubbles on the structure of turbulence in a turbulent boundary layer

The time-averaged velocity and turbulence intensity distributions were measured by a laser Doppler velocimeter in a turbulent boundary layer filled with microbubbles. The void fraction distribution was also measured using a fiber-optic probe. The velocity decreased in the region below 100 wall units with an increase in bubble density. This led to a decrease in the velocity gradient at the wall, which was consistent with a decrease in shearing stress on the wall. The turbulence intensity in the buffer layer increased at a low microbubble density, and then began to decrease with an increasing microbubble density. Based on the present measurements, the mechanism of turbulence reduction by microbubbles is discussed and a model is proposed. Received for publication on Dec. 3, 1999; accepted on April 18, 2000     

微气泡减阻的二相湍流边界层模型计算

利用边界层中注入微气泡来降低固体壁面摩擦阻力是一项有实用价值且有待进一步研究的减阻技术.本文以含微气泡的湍流边界层为研究对象,以气液二相流理论为基础,对二相湍流边界层控制方程进行了详细讨论,并利用差分网格法进行了求解,计算结果表明微气泡具有明显的减阻效果.最后对影响减阻效果的因素进行了分析,指出了今后的研究方向.     

基于CFD技术的滑行艇静水阻力计算

滑行艇作为一种高性能船舶,被广泛地应用于军事、运输、娱乐等方面,而滑行艇的阻力性能备受人们关注。本文基于CFD技术开展了滑行艇阻力计算研究,分析了时间步长与近壁面网格划分对计算精度的影响。数值计算结果表明,时间步长和近壁面网格划分的不同主要影响的是滑行艇的摩擦阻力以及艇底的气液分布情况。选择合理的参数设置可以使数值计算结果与试验较好地吻合,证明了CFD技术计算滑行艇阻力的可行性。     

气泡减阻技术研究进展

总结气泡减阻技术研究的发展历程,介绍微气泡减阻和大气泡减阻技术研究进展,介绍气泡减阻技术的实船应用情况,指出存在的问题和重点技术研究方向.     

Characteristics of a bubble jet near a vertical wall

A numerical model of a coupled bubble jet and wall was built on the assumption of potential flow and calculated by the boundary integral method. A three-dimensional computing program was then developed. Starting with the basic phenomenon of the interaction between a bubble and a wall, the dynamics of bubbles near rigid walls were studied systematically with the program. Calculated results agreed well with experimental results. The relationship between the Bjerknes effect of a wall and characteristic parameters was then studied and the calculated results of various cases were compared and discussed with the Blake criterion based on the Kelvin-impulse theory. Our analyses show that the angle of the jet＇s direction and the pressure on the rigid wall have a close relationship with collapse force and the bubble＇s characteristic parameters. From this, the application range of Blake criterion can be determined. This paper aims to provide a basis for future research on the dynamics of bubbles near a wall.     

减压条件下深水圆柱附近气泡的运动特性

为深入探究气泡与附近潜艇的相互作用规律,采用电火花气泡发生装置与高速摄像分析系统,在减压环境中对近似潜艇的深水圆柱周围不同位置的气泡进行研究。并以气泡通过圆柱结构中心的铅垂线为界的分类情况讨论深水圆柱周围气泡的动态变化。实验结果表明：圆柱结构附近的气泡的运动主要取决于气泡与圆柱结构间相对的攻角以及气泡质心到就近圆柱壁面的无量纲距离这两种因素的影响。在此基础上同时考虑减压环境中浮力的因素分析了浮力参数及攻角和爆距半径比对气泡第一周期运动时间和位移的影响,并总结圆柱下侧气泡不同爆距半径比下射流宽度的作用范围变化规律。研究结论为圆柱状的潜艇结构在水下受到攻击时的毁伤状态提供了参考。     

Simple Lagrangian formulation of bubbly flow in a turbulent boundary layer (bubbly boundary layer flow)

For the theoretical consideration of a system for reducing skin friction, a mathematical model was derived to represent, in a two-phase field, the effect on skin friction of the injection of micro air bubbles into the turbulent boundary layer of a liquid stream. Based on the Lagrangian method, the equation of motion governing a single bubble was derived. The random motion of bubbles in a field initially devoid of bubbles was then traced in three dimensions to estimate void fraction distributions across sections of the flow channel, and to determine local bubble behavior. The liquid phase was modeled on the principle of mixing length. Assuming that the force exerted on the liquid phase was equal to the fluid drag generated by bubble slip, an equation was derived to express the reduction in turbulent shear stress. Corroborating experimental data were obtained from tests using a cavitation tunnel equipped with a slit in the ceiling from which bubbly water was injected. The measurement data provided qualitative substantiation of the trend shown by the calculated results with regard to the skin friction ratio between cases with and without bubble injection as function of the distance downstream from the point of bubble injection.List of symbols B law of wall constant - C _f local coefficient of skin friction - C _f0 local coefficient of skin friction in the absence of bubbles - d _b bubble diameter [m] - g acceleration of gravity [m/s²] - k ₁ k₄ proportional coefficient - k _L turbulent energy of the liquid phase [m²/s²] - L representative length [m] - l _b mean free path of a bubble [m] - m _A added mass of a single bubble [kg] - m _b mass of a single bubble [kg] - N _x ,N _y ,N _z force perpendicular to the wall or ceiling exerted on a bubble adhering to that wall or ceiling [N] - P absolute pressure [Pa] - Q _G rate of air supply [/min] - q _L ⁽ⁱ⁾ turbulent velocity at the ith time increment [m/s] - R> _ex Reynolds number defined by Eq. 32 - T _*L integral time scale of the liquid phase [s] - U velocity of the main stream [m/s] - ,¯v,¯w time-averaged velocity components [m/s] - u,v,w turbulent velocity components [m/s] - û _L ,v_L root mean square values of liquid phase turbulence components in thex- and y-directions [m/s] - V volume of a single bubble [m³] - X,Y,Z components of bubble displacement [m] - x _s ,y _s ,z _s coordinate of a random point on a sphere of unit diameter centered at the coordinate origin - root mean square of bubble displacement in they-direction in reference to the turbulent liquid phase velocity [m] - local void fraction - _m mean void fraction in a turbulent region - regular random number - R _v increment of the horizontal component of the force acting on a single bubble, defined by Eq. 22 [N] - t time increment [s] - ₁ reduction of turbulent stress [N/m²] - _L rate of liquid energy dissipation [m²/s³] - _m coefficient defined by Eq. 30 - law of wall constant in the turbulent region in absence of bubbles - ₁ law of wall constant in the turbulent region in presence of bubbles     

Long-term prediction method for the green water load and volume for an assessment of the load line

 To develop a practical prediction method for the green water load and volume on the bow deck in irregular waves, model tests were conducted using a tanker and a cargo ship on a domestic Japanese voyage. The relation between green water load and relative water height at the stem was considered. Based on the finding that the maximum value of the green water load is proportional to the square of the maximum value of the water elevation over the bow top, the probability density functions of the green water load and volume in short-term predictions were proposed. It was verified that the proposed functions show good agreement with the measured distributions, especially in the tails, and were better than conventional functions. Using these functions, long-term predictions of the green water load were carried out. It was confirmed that the present method is more rational than the conventional one for estimating the long-term probability of the green water load. An assessment of the bow height of a domestic Japanese ship from the viewpoint of deck wetness was carried out using these prediction methods. This research was used as the technical background for the revision of domestic rules on load lines, which was enforced in October 2001. Received: July 19, 2002 / Accepted: October 30, 2002 Acknowledgment. Some of the present study was carried out as part of a cooperation project (RR45) with the Shipbuilding Research Association of Japan, supported by the Nippon Foundation. Address correspondence to: Y. Ogawa (e-mail: ogawa@nmri.go.jp)     

中速三体船型PCC的可行性研究(英文)

In the present paper, a new trimaran Pure Car Carrier (PCC) is proposed and a feasibility study on the ship is carried out. In this study, first, the effective horse power (EHP)/car of the PCC running in still water is predicted. By comparing the predicted EHP/car with that of a conventional mono-hull PCC, it is found that the trimaran PCC is superior to the conventional mono-hull PCC at rather higher speed. As ship speed increases, the reduction of the resistance of the trimaran is bigger. It is also found that at common service speed of PCCs, the EHP/car of a small PCC is lower than that of a conventional PCC. Secondly, the optimal L/B of a main-hull of the trimaran PCC in still water is determined. The optimal L/B of the main-hull varies with ship speed and size because the wave resistance decreases but the frictional resistance increases as L/B of the hull increases. As ship size increases, the optimal L/B of the main-hull of the trimaran PCC decreases.     

气泡减阻实验研究

为研究气泡对航行体侧壁的减阻效果，设计一套试验装置。利用该装置，采用壁面通气法研究不同通气量、不同来流速度和不同通气角度对气泡减阻效果的影响。实验结果表明：在通气量和通气角度不变的条件下，减阻率随来流速度增加而增加；在来流速度和通气角度不变的条件下，随通气量增加，减阻率先增加然后变小；在来流速度和通气量不变的条件下，减阻率随通气角度的增大而减少；在一定的条件下，取得了20．34％的减阻率。     

Broaching prediction in the light of an enhanced mathematical model,with higher-order terms taken into account

 We have attempted to develop a more consistent mathematical model for capsizing associated with surf-riding in following and quartering waves by taking most of the second-order terms of the waves into account. The wave effects on the hull maneuvring coefficients were estimated, together with the hydrodynamic lift due to wave fluid velocity, and the change in added mass due to relative wave elevations. The wave effects on the hydrodynamic derivatives with respect to rudder angles were estimated by using the Mathematical Modelling Group (MMG) model. Then captive ship model experiments were conducted, and these showed reasonably good agreements between the experiments and the calculations for the wave effects on the hull and the rudder maneuvring forces. It was also found that the wave effects on restoring moments are much smaller than the Froude–Krylov prediction, and the minimum restoring arm appears on a wave downslope but not on a wave crest amidship. Thus, an experimental formula of the lift force due to the heel angle of the ship is provided for numerical modelling. Numerical simulations were then carried out with these second-order terms of waves, and the results were compared with the results of free-running model experiments. An improved prediction accuracy for ship motions in following and quartering seas was demonstrated. Although the boundaries of the ship motion modes were also obtained with both the original model and the present one, the second-order terms for waves are not so crucial for predicting the capsizing boundaries themselves. Received: June 20, 2002 / Accepted: October 10, 2002 Acknowledgments. This research was supported by a Grant-in-Aid for Scientific Research of the Ministry of Education, Culture, Sports, Science and Technology of Japan (No. 13555270). The authors thank Prof. N. Rakhmanin of the Krylov Ship Research Institute for providing the Russian literature, as well as Mr. H. Murata of NHK (Japan Broadcasting Corporation) for translating it into Japanese. Address correspondence to: N. Umeda (e-mail: umeda@naoe.eng.osaka-u.ac.jp)     

两相流相互作用下的船模微气泡减阻性能数值仿真

基于两相流理论,对集装箱船船模进行微气泡减阻性能数值仿真研究。分析船模在不同喷气口位置、流体含气率、不同喷射速度和角度下的流场分布和减阻效果。结果表明：当流体微气泡含量过高时,会增大船体摩擦阻力,流体含气率为10%~20%时,有较好的减阻效果,喷气口位置在距离球鼻艏后缘约三分之一船长附近具有较高的空气覆盖率,喷气口喷气方向垂直向下所产生的减阻效果比喷气方向偏向船侧效果要好,且船速较高时可允许较高含气率流体注入流场并增大减阻效果。     

Hydrodynamic optimization of ship hull forms in shallow water

A computational method for improving hull form in shallow water with respect to wave resistance is presented. The method involves coupling ideas from two distinct research fields: numerical ship hydrodynamics and nonlinear programming techniques. The wave resistance is estimated by means of Morinos panel method, which is extended to free surface flow and considers the influence of finite depth on the wave resistance of ships. This is linked to the optimization procedure of the sequential quadratic programming (SQP) technique, and an optimum hull form can be obtained through a series of iterations giving some design constraints. Sinkage is an important factor in shallow water, and this method considers sinkage as a hydrodynamic design constraint. The optimization procedure developed is demonstrated by selecting a Wigley (C _B = 0.444) hull and the Series 60 (C _B = 0.60) hull, and new hull forms are obtained at Froude number 0.316. The Froude number specified corresponds to a lower than critical speed since most of the ships operating in shallow water move below their critical speed. The numerical results of the optimization procedure indicate that the optimized hull forms yields a reduction in wave resistance.     

Elastic response characteristics of a very large floating structure in waves moored inside a reef

 This article describes the results of hydraulic model tests of the elastic response of a very large floating structure (VLFS) moored inside a reef in an isolated island. The distributions of strains and vertical displacements due to the elastic response of the VLFS were measured. The response characteristics were strongly affected by deformed nonlinear waves inside the reef. A two-step analytical method to compute the elastic response of a VLFS is proposed, and its validity is verified using the results of the hydraulic model tests. Received: May 2, 2002 / Accepted: March 17, 2003 RID="*" ID="*" Address correspondence to: S. Shiraishi (shiraishi@pari.go.jp) Acknowledgment. This study was supported by the Program for Promotions of Fundamental Transport Technology Research from the Corporation for Advanced Transport and Technology (CATT).     

Comparison of Microbubble and Air Layer Injection with Porous Media for Drag Reduction on a Self-propelled Barge Ship Model

Ship resistance issues are related to fuel economy, speed, and cost efficiency. Air lubrication is a promising technique for lowering hull frictional resistance as it is supposed to modify the energy in the turbulent boundary layer and thereby reduce hull friction. In this paper, the objective is to identify the optimum type of air lubrication using microbubble drag reduction (MBDR) and air layer drag reduction (ALDR) techniques to reduce the resistance of a 56-m Indonesian self-propelled barge (SPB). A model with the following dimensions was constructed: length L?=?2000 mm, breadth B?=?521.60 mm, and draft T?=?52.50 mm. The ship model was towed using standard towing tank experimental parameters. The speed was varied over the Froude number range 0.11–0.31. The air layer flow rate was varied at 80, 85, and 90 standard liters per minute (SLPM) and the microbubble injection coefficient over the range 0.20–0.60. The results show that the ship model using the air layer had the highest drag reduction up to a maximum of 90%. Based on the characteristics of the SPB, which operates at low speed, the optimum air lubrication type to reduce resistance in this instance is ALDR.     

A hybrid analytical model of sterilization effect on marine bacteria using microbubbles interacting with shock wave

The present paper reports on a hybrid analytical model consisting of a biological probability model and a physical impact model proposed to predict the cell viability ratio of a sterilization method for marine bacteria using microbubbles interacting with a shock wave. The physical impact model of interaction between a microbubble and a shock wave is developed on the basis of the bubble motion analysis with experimental pressure data and a one-dimensional numerical simulation. An underwater shock wave produced by microbubble motion is simulated using a second-order finite differential scheme by means of bubble surface velocity variation obtained from Herring’s bubble motion equation, and the radius of the sterilized space around a microbubble is estimated by the critical pressure that just causes cell wall damage of marine bacteria. The sterilization effect predicted by the present hybrid analysis shows a good agreement with the bio-experimental result.     

Structural safety assessment of a pontoon-type VLFS considering damage to the breakwater

 A structural safety assessment of a pontoon-type very large floating structure (VLFS) surrounded by a gravity-type breakwater was carried out for extreme wave conditions by considering the damage to the breakwater. Bending and shear collapses are considered to be a failure mode of the floating structure, while overturning damages the breakwater. The probability of the breakwater overturning, and the transmitted wave height before and after damage to the breakwater, are evaluated using design formulae for port and harbor facilities in Japan. The ultimate bending and shear strengths of the floating structure are calculated by the idealized structural unit method (ISUM) and FEM, respectively. The calculated failure probability for the floating structure is compared with the specified target safety level. It was found that the floating structure under consideration is most likely to fail by bending in transverse waves, and that the corresponding failure probability satisfies the target level. Received: September 12, 2002 / Accepted: October 4, 2002 Acknowledgment. The authors are grateful to Dr. Shigeo Ohmatsu, National Maritime Research Institute, Japan, for allowing us to use the program of hydroelastic response analysis. Address correspondence to: M. Fujikubo (e-mail: fujikubo@naoe.hiroshima-u.ac.jp) Updated from the Japanese original, which won the 2002 SNAJ prize (J Soc Arthit Jpn 2002;190:337–345)     

Experimental study on microbubble ejection method for frictional drag reduction

The formation of air bubbles ejected through a single hole in a flat plate was observed in uniform flow of 2–10 m/s It was confirmed that the size of the air bubbles was governed by main flow velocity and air flow rate. According to previous experiments, the size of the bubbles is an important factor in frictional drag reduction by microbubble ejection. Usually bubbles larger than a certain diameter (for example 1 mm) have no effect on frictional drag reduction. Three different methods were proposed and tested to generate smaller bubbles. Among them, a 2D convex (half body of an NACA 64-021 section) with ejection holes at the top was the best and most promising. The diameter of the bubbles became about one-third the size of the reference ejection on a flat plate. Moreover, the bubble size did not increase with increasing flow rate. This is a favorable characteristic for practical purposes. The skin friction force was measured directly with a miniature floating element transducer, and decreased drastically by microbubble ejection from the top of the 2D convex shape.