Model test of the surface and submerged vehicles with the micro-bubble drag reduction

This study is based on the effective experiment observation and measuring technology to discuss the interaction influence between liquid turbulent boundary layer and a crowded group micro-bubbles. It is in order to understand and quantify the micro-bubbles clouds inside the turbulent boundary layer to eliminate the capacity of skin friction drag. Whenever the micro-bubbles are over supplied, pile up effect happened which makes micro-bubbles to integrate to each other as a large-size air film. Although they still have the drag reduction effect, the efficiency of drag reduction slowed down at this transition period. In the experiment of vertical type circulating water tunnel, when 1 μm porous medium is at 7 m/s flow speed, the C_v value at 0.056 has the best drag reduction efficiency of 26%. While 10 μm porous medium is at the same flow speed, the drag reduction efficiency is only around 23%.     

近岸破碎波水体掺气及气泡输运过程研究

为准确探讨破碎波作用下气体如何卷入以及气泡的形成与输运特性, 文章结合粒子图像测速技术(particle image velocimetry, PIV)、高速相机和气泡测量系统, 以及基于Navier-Stokes方程的三维数值模型对气泡形成及其运动过程进行研究。研究结果表明: 文章建立的数值模型能合理地捕捉到破碎波作用下气体的卷入及其输运过程; 波浪的破碎会形成较大的气腔, 其破裂过程又将产生大量的气体微团; 气泡会增加水体的紊动, 造成水体与空气交界面附近形成大量的漩涡以及水体的飞溅; 气泡的破裂会消耗大量的水体能量, 同时发现较大的紊动动能与气泡的生成有关, 且气泡数随平均紊动动能的增加呈线性增长关系。     

Numerical research on drag reduction by ventilated partial cavity based on two-fluid model

This paper presents a numerical study on the drag reduction mechanism created by a ventilated partial cavity and its associated effects by the downstream dispersed microbubbles. A semi-empirical approach is introduced to model the discrete interface of the ventilated cavity and its complex gas leakage rate induced by the local turbulent shear stress. Based on the Eulerian–Eulerian two-fluid modeling framework, a population balance approach based on MUltiple-SIze-Group (MUSIG) model is incorporated to simulate the dynamical effects of bubbly flow along the test body. Particular attention is also directed to grasp a better understanding of the size evolution of microbubble and its associated effects on drag reduction. Model predictions are validated against three experimental measurements carried out in a high-speed water tunnel by Schauer (2003) and Wosnik et al. (2005). Close examination of the flow structures, gas void fraction distributions and its resultant density ratio provides valuable insights on the complex physical phenomenon, helping to consolidate idea to maximize the drag reduction for ventilated cavitating vehicles.     

长江口湍流剖面的观测与分析

小尺度湍流过程对河口物质输运与能量交换至关重要。受传统观测方法的限制, 河口浅水区域的剖面观测资料至今较为匮乏, 进而限制了湍流过程的研究。为此, 采用新型5波束声学多普勒流速剖面仪(Nortek Signature 1 000 kHz AD2CP)在长江口开展湍流剖面观测, 获取高频、低噪的高质量湍流剖面数据, 并与声学多普勒点式流速仪(acoustic doppler velocimeters, ADV)同步观测的数据进行对比。结果表明, 通过AD2CP与ADV获得的近底部边界层摩阻流速u_*、拖曳系数C_d、雷诺应力S_R等特征参数基本一致, 底摩擦与波浪能量为河口区域湍动能的主要输入源。湍流垂向结构存在显著的非局地平衡, 即温盐等斜压作用引起的浮力通量、对流项以及强波浪作用影响的脉动压力做功、黏性输运等因素可能为长江口湍流非局地平衡的主因。     

Large eddy simulations of a circular cylinder at Reynolds numbers surrounding the drag crisis

Large eddy simulations of the flow around a circular cylinder at high Reynolds numbers are reported. Five Reynolds numbers were chosen, such that the drag crisis was captured. A total of 18 cases were computed to investigate the effect of gridding strategy, turbulence modelling, numerical schemes and domain width on the results. It was found that unstructured grids provide better resolution of key flow features, when a ‘reasonable’ grid size is to be maintained.When using coarse grids for large eddy simulation, the effect of turbulence models and numerical schemes becomes more pronounced. The dynamic mixed Smagorinsky model was found to be superior to the Smagorinsky model, since the model coefficient is allowed to dynamically adjust based on the local flow and grid size. A blended upwind-central convection scheme was also found to provide the best accuracy, since a fully central scheme exhibits artificial wiggles, due to dispersion errors, which pollute the solution.Mean drag, fluctuating lift Strouhal number and base pressure are compared to experiments and empirical estimates for Reynolds numbers ranging from 6.31 × 10⁴ to 5.06 × 10⁵. In terms of the drag coefficient, the drag crisis is well captured by the present simulations, although the other integral quantities (rms lift and Strouhal number) show larger discrepancies. For the lowest Reynolds number, the drag is seen to be more sensitive to the domain width than the spanwise grid spacing, while at the higher Reynolds numbers the grid resolution plays a more important role, due to the larger extent of the turbulent boundary layer.     

Reynolds number variation in oscillatory boundary layers: Part I. Purely oscillatory motion

A numerical model based upon a low Reynolds number turbulence closure is proposed to study Reynolds number variation in reciprocating oscillatory boundary layers. The model is used to compute the boundary layer for flow regimes ranging from smooth laminar to rough turbulent. Criteria for fully developed turbulence are derived for walls of the smooth and rough types. In particular, a new criterion to identify the rough turbulent regime is determined based on the time-averaged turbulence intensity. The reliability of the present model is assessed through comparisons with detailed experimental data collected by other investigators. The model globally improves upon standard high Reynolds number closures. Variation through the wave cycle of the main flow variables (ensemble-averaged velocity, shear stress, turbulent kinetic energy) is remarkably well-predicted for smooth walls. Predictions are satisfactory for rough walls as well. Yet, the turbulence level in the rough turbulent regime is overpredicted in the vicinity of the bed.     

Numerical investigation on the gas entrainment of ventilated partial cavity based on a multiscale modelling approach

Ventilated cavitation which is acknowledged as an efficient drag reduction technology for underwater vehicle is characterised by the very disparate length and time scales, posing great difficulty in the application of this technology. A multiscale numerical approach which integrates a sub-grid air entrainment model into the two-fluid framework is proposed in this paper to resolve the complex flow field created by ventilated cavity. Simulations have been carried out for the partially ventilated cavity underneath flat plate, with special efforts putting on understanding the gas entrainment at the cavity tail and the bubble dispersion process downstream. The flow parameters including the void fraction, the bubble velocity and the bubble size distributions in and downstream of the ventilated cavity are fully investigated. Comparisons between the numerical results with the experimental data are in satisfactory agreement, demonstrating the potential of the proposed methodology. The ventilation rate effect on the cavity shape and bubbly flow parameters are further investigated, obtaining the law of bubble dispersion and the bubble size evolution. This research not only provide a useful method for the investigation on the multiscale multiphase flow, but also give insight on understanding the combined drag reduction mechanism resulted from large-scale cavity and microbubbles.     

Micro-bubble resistance reduction on a model SES catamaran

The concept of micro-bubble ejection is emerging as an attractive method for viscous drag reduction for high speed craft. Most of the drag reduction literature presents results for micro-bubble ejection above/below horizontal surfaces. Nevertheless, modern high-speed hull forms are slender with vertical surfaces. This paper presents results of model tests with micro-bubble ejection to the vertical sides of a high speed model with a non-wetting coating. The coating yielded drag reduction of 4–6% and micro-bubbles with coating reduced overall drag by 4–11%. The local coefficient of friction with micro-bubbles is analyzed and found comparable to vertical flat plate results of [Madavan et al., 1984].     

The influence of turbulence on drag

The use of the Reynolds number as the only correlating factor for drag force measurements may be inadequate in circumstances involving highly turbulent flows. The results of previous investigations relating to the effects of turbulence scale and intensity are examined. Of special interest is the possibility of a drag minimum, even at low Reynolds number, for a free-stream turbulence intensity of about 5%. This appears to be the result of interaction between the free stream and the boundary layer. As intensity increases beyond 5%, the minimum may be succeeded by an increase in drag to values exceeding the laminar flow values. Further elucidation of the subject is required, particularly because of its importance in various problems related to geophysical flows.     

抛弃式温盐探头阻力系数的数值模拟及其实验验证

采用时变雷诺方程结合κ-ε湍流模型,模拟探头在下沉过程中的运动状态。根据模拟所得探头附近速度与压力分布,积分出阻力系数,由此分析阻力系数随雷诺数变化规律。同时进行物理实验,并与数值结果比较,二者吻合较好。     

Numerical study on the characteristics of flow field and wave propagation near submerged breakwater on slope

In this study,characteristics of flow field and wave propagation near submerged breakwater on a sloping bed are investigated with numerical model. The governing equations of the vertical twodimensional model are Reynolds Averaged Navier Stokes equations. The Reynolds stress terms are closed by a nonlinear k ε turbulence transportation model. The free surface is traced through the PILC-VOF method. The proposed numerical model is verified with experimental results. The numerical result shows that the wave profile may become more asymmetrical when wave propa-gates over breakwater. When wave crest propagates over breakwater,the anticlockwise vortex may generate. On the contrary,when wave hollow propagates over breakwater,the clockwise vortex may generate. Meanwhile,the influenced zone of vortex created by wave crest is larger than that created by wave hollow. All the maximum values of the turbulent kinetic energy,turbulent dissi-pation and eddy viscosity occur on the top of breakwater. Both the turbulent dissipation and eddy viscosity increase as the turbulent kinetic energy increases. Wave energy may rapidly decrease near the breakwater because turbulent dissipation increases and energy in lower harmonics is transferred into higher harmonics.     

A wave-resolving model for nearshore suspended sediment transport

This paper presents a wave-resolving sediment transport model, which is capable of simulating sediment suspension in the field-scale surf zone. The surf zone hydrodynamics is modeled by the non-hydrostatic model NHWAVE (Ma et al., 2012). The turbulent flow and suspended sediment are simulated in a coupled manner. Three effects of suspended sediment on turbulent flow field are considered: (1) baroclinic forcing effect; (2) turbulence damping effect and (3) bottom boundary layer effect. Through the validation with the laboratory measurements of suspended sediment under nonbreaking skewed waves and surfzone breaking waves, we demonstrate that the model can reasonably predict wave-averaged sediment profiles. The model is then utilized to simulate a rip current field experiment (RCEX) and nearshore suspended sediment transport. The offshore sediment transport by rip currents is captured by the model. The effects of suspended sediment on self-suspension are also investigated. The turbulence damping and bottom boundary layer effects are significant on sediment suspension. The suspended sediment creates a stably stratified water column, damping fluid turbulence and reducing turbulent diffusivity. The suspension of sediment also produces a stably stratified bottom boundary layer. Thus, the drag coefficient and bottom shear stress are reduced, causing less sediment pickup from the bottom. The cross-shore suspended sediment flux is analyzed as well. The mean Eulerian suspended sediment flux is shoreward outside the surf zone, while it is seaward in the surf zone.     

A two-fluid model of turbulent two-phase flow for simulating turbulent stratified flows

In this paper, a two-fluid model of turbulent two-phase flow is used to simulate turbulent stratified flows. This is a unified multi-fluid model for the motion of each phase in the flow, whose turbulent transport is closed by a two-phase k– model. The exchanges of mass, momentum and energy between the two phases are fully accounted for in the simulation. For illustration, a case of turbulent stratified flow with strong buoyancy effects, for which extensive experimental data are available, is selected for examination. It is shown that the numerical results agree well with the experimental data.     

On the application of a turbulence closure modified model to the description of the density jump evolution in a stably stratified medium

The self-similar turbulent density jump evolution has been studied in the scope of a turbulence closure modernized theory which takes into account the anisotropy and mutual transformation of the turbulent fluctuation kinetic and potential energy for a stably stratified fluid. The numerical calculation, performed using the equations for the average density and kinetic and potential energies of turbulent fluctuations, indicates that the vertical profiles of the buoyancy frequency, turbulence scale, and kinetic and potential energies drastically change when the turbulence anisotropy is strong. The vertical profiles of the corresponding energy and spatial discontinuity parameters, calculated at a weaker anisotropy, indicate that similar drastic changes are absent and a qualitative agreement exists with the known analytical solution, which describes the density jump evolution in a freshwater basin and was obtained previously [5, 8] in the scope of a turbulence local-similarity hypothesis applied in combination with the budget equation for the turbulent fluctuation kinetic energy.     

Basic Laws in Mechanics of Turbulent Flows

- Turbulent flow is a basic form of fluid motion widely observed in nature. In hydraulic engineering, especially in the study of sediment movement, turbulence is a key problem. In this paper, based on the stochastic theory of wall turbulence developed by the author and the results by other investigators, fluc-tuation and mean structures and drag coefficient for Newtonian and drag reduction flows in all states (laminar, transitional, turbulent) and in all regions (smooth, transitional, rough) are theoretically discussed in detail. General laws for laminar and turbulent flows obtained by the author are verified by the experimental results obtained by others, and there is good agreement between them.     

Numerical study on evolution of an internal solitary wave across an idealized shelf with different front slopes

Numerical simulations are performed to investigate the influence of variable front slopes on flow evolution and waveform inversion of a depression ISW (internal solitary wave) over an idealized shelf with variable front slopes. A finite volume based on Cartesian grid method is adopted to solve the Reynolds averaged Navier-Stokes equations using a k-ε model for the turbulent closure. Numerical results exhibit the variations of several pertinent properties of the flow field, in the case with or without waveform inversion on the horizontal plateau of an obstacle. The clockwise vortex is stronger than the counterclockwise one, almost throughout the wave-obstacle interaction. Analysis of the turbulent energy budget reveals that the turbulent production term in the governing equations dominates the wave evolution during a wave-obstacle interaction; otherwise the buoyancy production term and the dissipation term due to viscosity within turbulent eddies play a major role in energy dissipation. In addition, the front slope affects mainly the process and reflection of the wave evolution but has less influence than other physical parameters. Moreover, total wave energy of the leading crest is smaller than that of the leading trough even in the cases with waveform inversion on the plateau.     

Turbulent mixing in stably stratified flows of the environment: The current state of the problem (Review)

Specific features of the turbulent transfer of the momentum and heat in stably stratified geophysical flows, as well as possibilities for including them into RANS turbulence models, are analyzed. The momentum (but not heat) transfer by internal gravity waves under conditions of strong stability is, for example, one such feature. Laboratory data and measurements in the atmosphere fix a clear dropping trend of the inverse turbulent Prandtl number with an increasing gradient Richardson number, which must be reproduced by turbulence models. Ignoring this feature can cause a false diffusion of heat under conditions of strong stability and lead, in particular, to noticeable errors in calculations of the temperature in the atmospheric boundary layer. Therefore, models of turbulent transfer must include the effect of the action of buoyancy and internal gravity waves on turbulent flows of the momentum. Such a strategy of modeling the stratified turbulence is presented in the review by a concrete RANS model and original results obtained during the modeling of stratified flows in the environment. Semiempirical turbulence models used for calculations of complex turbulent flows in deep stratified bodies of water are also analyzed. This part of the review is based on the data of investigations within the framework of the large international scientific Comparative Analysis and Rationalization of Second-Moment Turbulence Models (CARTUM) project and other publications of leading specialists. The most economical and effective approach associated with modified two-parameter turbulence models is a real alternative to classical variants of these models. A class of test problems and laboratory and full-scale experiments used by the participants of the CARTUM project for the approbation of numerical models are considered.     

Process-based model for nearshore hydrodynamics, sediment transport and morphological evolution in the surf and swash zones

Hydrodynamics and sediment transport in the nearshore zone were modeled numerically taking into account turbulent unsteady flow. The flow field was computed using the Reynolds Averaged Navier–Stokes equations with a k–ε turbulence closure model, while the free surface was tracked using the Volume-Of-Fluid technique. This hydrodynamical model was supplemented with a cross-shore sediment transport formula to calculate profile changes and sediment transport in the surf and swash zones. Based on the numerical solutions, flow characteristics and the effects of breaking waves on sediment transport were studied. The main characteristic of breaking waves, i.e. the instantaneous sediment transport rate, was investigated numerically, as was the spatial distribution of time-averaged sediment transport rates for different grain sizes. The analysis included an evaluation of different values of the wave friction factor and an empirical constant characterizing the uprush and backwash. It was found that the uprush induces a larger instantaneous transport rate than the backwash, indicating that the uprush is more important for sediment transport than the backwash. The results of the present model are in reasonable agreement with other numerical and physical models of nearshore hydrodynamics. The model was found to predict well cross-shore sediment transport and thus it provides a tool for predicting beach morphology change.     

The influence of bubble clouds on acoustic propagation in the surfzone

High-frequency propagation close to an active surf line is explored with 12and 100-kHz propagation paths together with measurements of bubble clouds, bubble size distributions, and waves. Breaking waves inject massive bubble plumes that are mixed downwards from the roller region by intense turbulence. If these injections follow one another at intervals less than the time taken for the bubbles to rise to the surface, acoustic signals will be continuously blocked, forming an acoustical barrier that effectively inhibits any propagation. Occasionally, waves break seaward of this barrier. In this case, dense bubble clouds are mixed down beneath the air entrainment zone, but there is sufficient time for them to disappear before succeeding breakers, allowing intermittent high-frequency propagation recharge the bubble field. The duration and shape of signal dropouts are then determined by the selective removal of bubbles by buoyancy and dissolution. In addition to turbulence created by the air entrainment process, a lower level of continuous background turbulence may be generated by interaction of residual currents with the wave boundary layer. Our observations illustrate the variable character of acoustic blocking by bubble clouds and serve as a basis for quantitative analysis of these effects with a 2D propagation model coupled to 2D models of bubble cloud evolution and background turbulence     

Study of the Neutral Ekman Flow Using an Algebraic Reynolds Stress Model

A recently developed fully explicit algebraic model of Reynolds stress and turbulent heat flux in a thermally stratified planetary atmospheric boundary layer without stratification has been used for a numerical study of the Ekman turbulent boundary layer over a homogeneous rough surface for different dimensionless surface Rossby numbers. A comparative analysis has been conducted for a closure model of the transport term in the prognostic equation of turbulent kinetic energy dissipation including third-order moments. Dependences of the total wind rotation angle on the Rossby number have been obtained. The calculated vertical profiles of mean velocity, turbulent stress, turbulent kinetic energy, surface-friction velocity, and boundary-layer height agree satisfactorily with observational and earlier obtained LES data.