含微气泡二相湍流边界层减阻理论计算

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

微气泡层减阻的边界层模型计算

以流体中光滑平板为例,建立微气泡层覆盖平板的气液简单边界层模型,采用平板层流理论,对给定的气泡喷射参数和流体特性参数计算了微气泡层状态下平板的摩擦阻力,结果显示微气泡层对平板有明显减阻效果;最后给出了气体边界层的速度分布和剪切应力分布.     

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

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

利用微气泡层减少平板阻力的边界层模型计算

以流体中光滑平板为例，建立微气泡层覆盖平板的气液简单边界层模型，采用平板层流理论，对给定的一系列气泡喷射速度参数和流体特性参数计算了微气泡状态下平板的摩擦阻力，结果显示微气泡对平板有明显减阻效果；最后给出了气体边界层的速度分布和剪切应力分布。     

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.     

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.     

Influence of microbubble diameter and distribution on frictional resistance reduction

Investigations into frictional resistance reduction by microbubbles were carried out using a two-dimensional channel. The flow velocity and the amount of air injected were varied, and the frictional resistance reduction was measured. The frictional resistance reduction increased with increasing mean void ratio. When the bubble diameter was changed, the influence on frictional resistance reduction was negligible. The influence of bubble distribution near a wall was also investigated. Although it is thought that the influence of microbubbles near a wall is large, further investigations are required. Received: February 6, 2002 / Accepted: April 30, 2002 Acknowledgments. This research was supported by a special research program of Toyo University, SR239 Research Committee of the Shipbuilding Research Association of Japan, and the Research Project on Smart Control of Turbulence. The authors express their sincere gratitude for these supports. Address correspondence to: H. Kato (kato@eng.toyo.ac.jp)     

沟槽面减阻效果影响因素及减阻机理的分析

采用雷诺平均N-S方程和RNG k-ε湍流模型计算V型沟槽面的湍流边界层流动和粘性阻力,研究了沟槽尖峰形状和雷诺数对减阻效果的影响规律,初步分析了沟槽面减阻机理.指出:沟槽尖峰处的圆角半径越小其减阻效果越好,沟槽斜面中下部的壁面应力随着圆角半径的减小而降低,但尖峰处的局部壁面应力会随之增大;来流速度对沟槽减阻率的影响很大,对于一种尺度的V型沟槽,存在着一个具有较好减阻效果的来流速度范围,而沟槽面在沿来流方向上的布置位置对减阻效果的影响非常小;沟槽尺度对减阻效果很剧烈;沟槽尖峰处生成的二次涡是产生减阻效果的根本原因,二次涡的强弱与沟槽减阻率的大小紧密相关.     

平面二维k-ε紊流模型不同壁函数的对比及研究

在平面二维紊流模型中壁函数起到连接流固的桥梁作用,克服了边界附近网格过密的困难,也是紊流模型应用中的一个重要环节,因此壁函数的选择对模型的模拟结果起着至关重要的作用。将3种不同形式的壁函数引入平面二维k-ε模型中,对丁坝绕流进行了数值计算,并将计算结果与实测数据进行了对比,从流场、流速、紊动动能分布及紊动耗散率等方面对3种壁函数进行优选验证,结果表明对数率壁函数优于其他两种形式的壁函数,其数值计算结果更接近实测数据。     

Experimental verification of shock sterilization for marine <Emphasis Type="Italic">Vibrio</Emphasis> sp. using microbubbles interacting with underwater shock waves

In this paper, shock sterilization using the motions of microbubbles induced by underwater shock waves is verified experimentally. A bio-experiment is carried out using marine Vibrio sp.. Underwater shock waves are produced by electric discharge in a semi-ellipsoidal discharge reflector. The shock waves are focused to increase the pressure that the generated microbubbles are exposed to. The microbubbles are generated independently. The microbubble generator can produce microbubbles of around 50 μm diameter by means of Kelvin–Helmholz instability and Venturi effect. Propagation behavior of shock waves and generation process of microbubbles are captured by high-speed camera. The experimental results show that the supply of microbubbles increases the potential of shock sterilization. In addition, it is found that shock waves without microbubbles also have the capacity of sterilization, and this means that cavitation bubbles generated behind converging shock waves contribute to inactivating marine bacteria.     

空心四面块体促淤堤近底水动力特性试验研究*

依据长江口某促淤工程的设计条件及物理模型试验方法对空心四面块体促淤堤近底水动力特性开展研究。结果表明：1）近底紊动强度最大值所在位置、近底水平时均流速最小值所在位置以及近底水平时均流动方向的转变位置基本重合;2）对于某种具体堤型而言,在相同单宽流量条件下,随着相对水深的变化（至少在本文试验范围内）,其紊动强度最大值所在位置比较稳定;3）与抛石加护面块体促淤堤相比,透水性较好的堤身结构近底紊动强度最大值所在位置会更加远离堤轴线,可使局部冲刷坑更加远离堤轴线,对堤身结构稳定有益。     

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.     

带附体潜艇绕流场数值计算与验证

采用κ-ε湍流模型和有限体积法用求解RANS方程的方法,对全附体的SUBOFF AFF-8潜艇模型三维粘性绕流场进行了数值模拟计算.一种结合增强型壁面函数的两层模型被采用,以便在近壁区域采用精细网格来模拟边界层内的流场流动特征.给出了雷诺数Re=1.2×107下模型表面纵中剖面线上半部分压力系数与摩擦阻力系数的分布,以及指挥台围壳与尾翼表面不同高度处的压力系数分布,同时计算了桨盘面处的尾流速度分布,并与有关文献给出的实测值进行了比较和验证,分布趋势基本一致,较相关文献给出的计算结果有相当大的改进,此外从指挥台围壳与尾翼后的速度矢量图观察到了与有关文献相同的涡旋流场结构,从而验证了该数值方法的有效性,为潜艇概念设计中评价设计方案优劣提供了重要的技术手段.     

An approach to estimation of near-surface turbulence and CO2 transfer velocity from remote sensing data

The air–sea CO₂ exchange is primarily determined by the boundary-layer processes in the near-surface layer of the ocean since it is a water-side limited gas. As a consequence, the interfacial component of the CO₂ transfer velocity can be linked to parameters of turbulence in the near-surface layer of the ocean. The development of remote sensing techniques provides a possibility to quantify the dissipation of the turbulent kinetic energy in the near-surface layer of the ocean and the air–sea CO₂ transfer velocity on a global scale. In this work, the dissipation rate of the turbulent kinetic energy in the near-surface layer of the ocean and its patchiness has been linked to the air–sea CO₂ transfer velocity with a boundary-layer type model. Field observations of upper ocean turbulence, laboratory studies, and the direct CO₂ flux measurements are used to validate the model. The model is then forced with the TOPEX POSEIDON wind speed and significant wave height to demonstrate its applicability for estimating the distribution of the near-surface turbulence dissipation rate and gas transfer velocity for an extended (decadal) time period. A future version of this remote sensing algorithm will incorporate directional wind/wave data being available from QUIKSCAT, a now-cast wave model, and satellite heat fluxes. The inclusion of microwave imagery from the Special Sensor Microwave Imager (SSM/I) and the Synthetic Aperture Radar (SAR) will provide additional information on the fractional whitecap coverage and sea surface turbulence patchiness.     

Measurement of air–water CO2 transfer at four coastal sites using a chamber method

We measured the air–water CO₂ flux in four coastal regions (two coral reefs, one estuary, and one coastal brackish lake) using a chamber method, which has the highest spatial resolution of the methods available for measuring coastal air–water gas flux. Some of the measurements were considerably higher than expected from reported wind-dependent relationships. The average k₆₀₀ values for Shiraho Reef, Fukido Reef, Fukido River, and Lake Nakaumi were 1.5 ± 0.6, 3.2 ± 0.3, 0.69 ± 0.26, and 2.2 ± 0.9 (mean ± S.D.) times larger than the wind-dependent relationships. Results were compared with current-dependent relationships and vertical turbulence intensity (VTI). VTI is an index of water-surface stirring and is calculated from near-surface vertical velocity. Although some measurements from the reefs and river closely matched those expected from wind-dependent relationships, others were considerably higher. All data were correlated with VTI and were qualitatively explained by bottom macro-roughness enhancement. In Lake Nakaumi, results tended to differ from the wind-dependent relationships, and the difference between the measured and expected gas-transfer velocity was correlated with biological DO changes and/or the intensity of density stratification. We found these factors to have important effects on coastal gas flux. In addition, the chamber method was an effective tool for evaluating coastal gas flux.     

Quantifying copepod kinematics in a laboratory turbulence apparatus

We describe application of a new apparatus that permits simultaneous detailed observations of plankton behavior and turbulent velocities. We are able to acquire 3D trajectories amenable to statistical analyses for comparisons of copepod responses to well-quantified turbulence intensities that match those found in the coastal ocean environment. The turbulence characteristics consist of nearly isotropic and homogeneous velocity fluctuation statistics in the observation region. In the apparatus, three species of copepods, Acartia hudsonica, Temora longicornis, and Calanus finmarchicus were exposed separately to stagnant water plus four sequentially increasing levels of turbulence intensity. Copepod kinematics were quantified via several measures, including transport speed, motility number, net-to-gross displacement ratio, number of escape events, and number of animals phototactically aggregating per minute. The results suggest that these copepods could control their position and movements at low turbulence intensity. At higher turbulence intensity, the copepods movement was dominated by the water motion, although species-specific modifications due to size and swimming mode of the copepod influenced the results. Several trends support a dome-shaped variation of copepod kinematics with increasing turbulence. These species-specific trends and threshold quantities provide a data set for future comparative analyses of copepod responses to turbulence of varying duration as well as intensity.     

航道整治工程软体排压载体的水力特性

软体排是航道整治工程中的基础结构,其水力特性直接影响航道整治工程的护滩效果。通过建立水槽模型,研究软体排压载体附近流场结构、回流区长度、近底流速、紊动流场、切应力等水力特性。结果表明,设置软体排压载体后,近底流速和回流区流速减小,上层流速和水流纵向紊动强度增加,近底层的切应力则明显减小。压载体的设置有利于调节近底水流流态、促进泥沙落淤,提高工程边滩守护效果。     

半潜式钻井平台风载荷CFD预报

在平台设计过程中,海洋平台受到的风载荷是稳性分析与结构设计需要考虑的重要因素。本文采用CFD技术,对某型半潜式钻井平台在作业工况下受到的风载荷进行了预报,分析了湍流模型、风速剖面形式、湍流强度分布等因素对半潜式钻井平台风载荷的影响。结果表明：不同湍流模型对海洋平台风载荷的预报偏差在2%以内,可以忽略湍流模型的影响；风速剖面对风载荷有较大影响,随着风速梯度的增加而增大；湍流强度分布对海洋平台风载荷的影响在6%以内。研究结论可为建立高精度的海洋平台风载荷数值预报方法提供技术支持。     

航道清礁对西江鱼类保护区紊流体的影响

河底裸露的礁石增加水体紊动,有助于鱼卵正常形成,大量存在礁石的河段往往是鱼类栖息和产卵的重要场所。礁石也存在碍航的可能,航道清礁一定程度上影响水体紊动强度,改变鱼类栖息环境。为此,在重要鱼类栖息河段进行航道清礁需充分论证。以紊动动能作为评价水体紊动强度指标,开展西江中游珍稀鱼类保护区河段三维紊流数值模拟,计算并统计航道清礁工程前后水体紊动动能及不同强度的紊流水体体积,分析其变化规律及原因。研究结果表明,断面中紊动动能强度与水深成正相关关系;清礁工程导致河段中弱紊动动能的紊流水体体积增加,强紊动动能的紊流体积减小;底层水体紊动动能的变化量高于表层;强紊动动能水体体积减少的变化率大于弱紊动动能水体体积增加的变化率。     

Late summer stratification, internal waves, and turbulence in the Yellow Sea

Microstructure profiling measurements at two locations in the Yellow Sea (a deeper central basin and a local shelf break) were analyzed focusing on tidal and internal-wave induced turbulence near the bottom and in the pycnocline. A classical three-layer density structure consisting of weakly stratified surface and bottom boundary layers and a narrow sharp pycnocline is developed by the end of warm season. Turbulence in the surface layer was not influenced by the tidal forcing but by the diurnal cycle of buoyancy flux and wind forcing at the sea surface. The enhanced dissipation and diffusivity generated by the shear stress at the seafloor was found in the water interior at heights 10–15 m above the bottom with a phase shift of ~ 5–6 m/h. No internal waves, turbulence, or mixing were detected in the pycnocline in the central basin, in contrast to the pycnocline near the local shelf break wherein internal waves of various frequencies were observed all the time. The thickness of the surface layer near the local shelf break slightly exceeded that of the bottom layer (20 vs. 18 m). A 5–6 m high vertical displacement of the pycnocline, which emerged during the low tide, was arguably caused by the passage of an internal soliton of elevation. During this episode, the gradient Richardson number decreased below 0.25 due to enhanced vertical shear, leading to local generation of turbulence with dissipation rates exceeding the background level by an order of magnitude.