Inleakage of a vortex-ring bunch onto a plane screen

A method is suggested for simulating axisymmetric laminar or turbulent flows formed during the motion of a vortex-ring bunch of given geometry and circulation toward a plane screen. Earlier, similar problems were simulated with the numerical solution of the Navier-Stokes equations for laminar flows. Turbulent flows have remained unconsidered until now. When a vortex ring approaches the screen, the secondary nonstationary flow is induced near the screen’s surface and this secondary flow causes the formation of the radial boundary layer (provided that air viscosity is taken into account). First, the medium spreads out from the critical point at the screen’s center with the negative pressure gradient along the radial coordinate and then detaches in the region of the positive pressure gradient. This radial wall flow and the corresponding boundary layer are considered in the quasi-stationary approximation. When the boundary layer detaches at successive instances, the flow is replenished with the radially moving secondary vortex rings whose circulations have the sign opposite to that of the circulation of the primary vortex ring. It is the interaction of the primary and secondary vortices that governs process dynamics, which differs substantially from that in the case when the formation of secondary vortices is disregarded. The suggested method is based on the method of discrete vortices (a perfect liquid) and the boundary-layer (laminar or turbulent) theory. During the development of the flow under investigation, the nonstationary ascending flow in the direction perpendicular to the screen’s plane is formed and then this flow decays and dissipates. Simulations for large Reynolds numbers corresponding to the formation of the turbulent boundary layer show that the velocity of ascending vortices in the plane of the initial vortex bunch is less than one-tenth of the initial velocity of the descending vortex ring. The boundary layer is introduced into calculations with the sole goal of determining the parameters of the secondary vortex rings formed during boundary-layer detachments. The interaction of the primary and secondary vortices is then considered within the framework of a perfect medium. Simulations for large Reynolds numbers corresponding to the formation of the turbulent boundary layer on the screen were correlated with the available data obtained in laboratory experiments for small Reynolds numbers. Qualitative agreement between the simulations and experiments is fairly satisfactory. The simulation for one combination of the circulation and vortex-ring geometry takes at most 10–15 min with the use of an average PC.     

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.     

海底滑坡对海洋单桩冲击力试验研究

为了研究海底滑坡对海洋单桩的冲击力大小,首先通过调整高岭土、粉砂的不同含量,得到不同流变特性、不同密度的碎屑流,采用Herschel-Bulkley模型和幂率模型对流体流变性质进行描述;随后利用自制海底滑坡模型槽,模拟碎屑流在不同流速和黏度下对模型桩的冲击;并结合流体力学理论,建立阻力系数与非牛顿流体雷诺数之间关系表达式。试验数据表明:碎屑流黏度和流速是影响海底滑坡冲击力的主要因素,海底滑坡冲击力随着泥浆黏度和流速的增加而增大。同时,考虑碎屑流剪切稀释特性,得到管桩阻力系数随雷诺数变化的拟合公式,为海洋桩基础设计提供参考。     

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.     

An Analytical Solution for the Run-Out of Submarine Debris Flows

Submarine debris flows have a significant impact on offshore and coastal facilities. The unique characteristics of submarine debris flows involve large mass movements and long travel distances over very gentle slopes. To improve our insight and knowledge of the basic mechanism behind submarine debris flows, an analytical model was derived for the mobility of submarine debris flows. This model takes into account the mass change of debris flows induced by deposition, stagnation pressure, and the topography of the depositional area. One case study on the Palos Verdes debris flow proves its ability to predict the run-out distance of a submarine debris flow to a reasonable level of accuracy. On the gentle slopes, the submarine debris flow progressively loses mass due to deposition, which in turn influences the flow velocity. In addition, the results show that the slope angle and spreading angle of the debris depositional zone play important roles in the sliding process.     

高雷诺数下近壁圆柱绕流数值模拟与分析

近壁圆柱绕流问题在海底悬跨管道的研究中具有重要的意义。在绕流阻力、升力以及海底土壤的耦合作用下,海底管道所发生的移位、悬跨等现象对于海底管道的安全运行构成了很大的威胁。正确预测各种绕流条件下管流之间的作用力是保证油气管道安全的首要任务。海底管道在极端海洋环境条件下的管、流相互作用为高雷诺数绕流问题,处于高雷诺数下的绕流模拟比处于低雷诺数下的绕流模拟要复杂很多,它需要更精细的网格以及合适的湍流模型。此文对处于悬跨状态下的海底管道进行数值研究,给出不同间隙比下海流绕流海底管道的流场结构形态,分析了间隙比对绕流阻力和绕流升力的影响,为进一步研究海底悬跨管道的受力和变形提供载荷边界数据。     

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.     

Wave-induced pore pressures around submarine pipelines

For the response of an undisturbed seabed to water waves there exists a number of solutions. Based on the concept of linear consolidation theory they all lead to the result that the pore pressure satisfies the Laplace equation, at least when incompressibility of the pore fluid is assumed. This suggests that the flow around an embedded object, e.g. a submarine pipeline, is also governed by the Laplace equation, either exactly or approximately. In this paper a number of such problems is solved using complex function theory and the technique of conformal mapping. The hydrodynamical force is calculated from the solution for the pore pressure. A good agreement is obtained with existing solutions. A comparison is also made with a numerical solution based on the finite-element method. It is concluded that the hydrodynamical force acting upon a submarine pipeline is about ten to thirty percent of the buoyancy of the pipe depending on the maximum wave load and the burial depth.     

Sedimentology, stratigraphic occurrence and origin of linked debrites in the West Crocker Formation (Oligo-Miocene), Sabah, NW Borneo

The West Crocker Formation (Oligocene–Early Miocene), NW Borneo, consists of a large (>20 000 km²) submarine fan deposited as part of an accretionary complex. A range of gravity-flow deposits are observed, the most significant of which are mud-poor, massive sandstones interpreted as turbidites and clast-rich, muddy sandstones and sandy mudstones interpreted as debrites. An upward transition from turbidite to debrite is commonly observed, with the contact being either gradational and planar, or sharp and highly erosive. Based on their repeated vertical relationship and the nature of the contact between them, these intervals are interpreted as being deposited from one flow event which consisted of two distinct flow phases: fully turbulent turbidity current and weakly turbulent to laminar debris flow. The associated bed is called a co-genetic turbidite–debrite, with the upper debrite interval termed a linked debrite. Linked debrites are best developed in the non-channellised parts of the fan system, and are absent to poorly-developed in the proximal channel-levee and distal basin floor environments. Due to outcrop limitations, the genesis of linked debrites within the West Crocker Formation is unclear. Based on clast size and type, it seems likely that a weakly turbulent to laminar debris-flow flow phase was present when the flow event entered the basin. A change in flow behaviour may have led to deposition of a sand-rich unit with ‘turbidite’ characteristics, which was subsequently overlain by a mud-rich unit with ‘debrite’ characteristics. Flow transformation may have been enhanced by the disintegration and incorporation into the flow of muddy clasts derived from the upstream channel floor, channel mouth or from channel-levee collapse. Lack of preservation of this debrite in proximal areas may indicate either bypass of this flow phase or that the available outcrops fail to capture the debris flow entry point. Establishing robust sedimentological criteria from a variety of datasets may lead to the increasing recognition of co-genetic turbidite-debrite beds, and an increased appreciation of the importance of bipartite flows in the transport and deposition of sediments in deepwater environments.     

The numerical study of wave-induced pore water pressure response in highly permeable seabed

The coupling numerical model of wave interaction with porous medium is used to study waveinduced pore water pressure in high permeability seabed.In the model,the wave field solver is based on the two dimensional Reynolds-averaged Navier-Stokes(RANS) equations with a k-ε closure,and Forchheimer equations are adopted for flow within the porous media.By introducing a Velocity-Pressure Correction equation for the wave flow and porous flow,a highly efficient coupling between the two flows is implemented.The numerical tests are conducted to study the effects of seabed thickness,porosity,particle size and intrinsic permeability coefficient on regular wave and solitary wave-induced pore water pressure response.The results indicate that,as compared with regular wave-induced,solitary wave-induced pore water pressure has larger values and stronger action on seabed with different parameters.The results also clearly show the flow characteristics of pore water flow within seabed and water wave flow on seabed.The maximum pore water flow velocities within seabed under solitary wave action are higher than those under regular wave action.     

Relevance of transition turbulent model for hydrodynamic characteristics of low Reynolds number propeller

The propeller of an Autonomous Underwater Vehicle (AUV) operates at low Reynolds number in laminar to turbulent transition region. The performance of these propellers can be calculated accurately using RANSE solver with γ − Re_θ transition model. In this study, the global and local hydrodynamic characteristics of open and ducted propeller are investigated using the γ − Re_θ transition model. The capability of the γ − Re_θ transition model to capture laminar to turbulent transition on the surface of the open propeller is demonstrated by comparison with published experimental results. The application of transition model for the propeller K_a-4-70 inside the duct 19A shows that the centrifugal forces are dominant at low Reynolds number and the flow is mainly directed in the radial direction. The transition model is able to predict complex flow physics such as leading-edge separation, tip leakage vortex, and the separation bubble on outer surface of the duct. The accurate prediction of these flow phenomenon can lead to correct calculation of global hydrodynamic forces and moments acting on the propeller at low Reynolds number.     

Numerical Investigation of Submarine Hydrodynamics and Flow Field in Steady Turn

This paper presents numerical simulations of viscous flow past a submarine model in steady turn by solving the Reynolds-Averaged Navier?Stokes Equations (RANSE) for incompressible, steady flows. The rotating coordinate system was adopted to deal with the rotation problem. The Coriolis force and centrifugal force due to the computation in a body-fixed rotating frame of reference were treated explicitly and added to momentum equations as source terms. Furthermore, velocities of entrances were coded to give the correct magnitude and direction needed. Two turbulence closure models (TCMs), the RNG model with wall functions and curvature correction and the Shear Stress Transport (SST) model without the use of wall functions, but with curvature correction and low-Re correction were introduced, respectively. Take DARPA SUBOFF model as the test case, a series of drift angle varying between 0° and 16° at a Reynolds number of 6.53×106 undergoing rotating arm test simulations were conducted. The computed forces and moment as a function of drift angle during the steady turn are mostly in close agreement with available experimental data. Though the difference between the pressure coefficients around the hull form was observed, they always show the same trend. It was demonstrated that using sufficiently fine grids and advanced turbulence models will lead to accurate prediction of the flow field as well as the forces and moments on the hull.     

Centrifuge study on the runout distance of submarine debris flows

ABSTRACT

In this paper, a series of mini-drum centrifuge experiments on motion of submarine debris flow, which are able to correctly reproduce the self-weight stresses and gravity-dependent processes, are presented. These tests were performed using artificial submarine clay with high water content, from 93 to 149%. The extremely low shear strength made the debris material behave as idealistic lubricating material when it was deposited, resulting in a linear relationship between water content and runout distance of strongly coherent debris flow. On the other hand, the dilation of the flow body and hydroplaning was observed for weakly coherent debris flow, which further increased the mobility of flow body. A densimetric Froude number Fr_d was used to indicate the threshold of hydroplaning, which occurs if the Fr_d is greater than 0.2. Finally, two simple analytical models based on prototype debris flow under 1?g condition was used to validate the experiment results, which further prove the effects of soft marine clay on the high mobility of submarine debris flow. On the other hand, when the water content exceeded 120%, the experiment results deviated from the analytical solution due to the effects of hydroplaning.     

Simulation of the nonlinear dynamic interactions between waves, a submerged breakwater and the seabed

This study employed direct numerical simulation to simulate the fully nonlinear interaction between the water waves, the submerged breakwater, and the seabed under differing wave conditions. In the numerical simulation, the laminar flow condition in the seabed was applied to evaluate the more exact fluid resistance acting on the porous media. Varying incident wave conditions were applied to the flow field resulting from the wave–structure–seabed interaction, and the variation in the pore water pressure beneath the submerged breakwater was investigated along the cross-section of the submerged breakwater. Structural safety and scouring were also considered on the basis of the numerical results for the flow field around the structure and the variation of the pore water pressure.     

Simulation of vortex-induced vibrations of a cylinder using ANSYS CFX

In this paper, vortex-induced vibrations of a cylinder are simulated by use of ANSYS CFX simulation code. The cylinder is treated as a rigid body and transverse displacements are obtained by use of a one degree of freedom spring damper system. 2-D as well as 3-D analysis is performed using air as the fluid. Reynolds number is varied from 40 to 16000 approx., covering the laminar and turbulent regimes of flow. The experimental results of （Khalak and Williamson, 1997） and other researchers are used for validation purposes. The results obtained are comparable.     

Two-range scaling for tip vortex cavitation inception

The suggested scaling concept for cavitation inception index takes into account boundary layer separation caused by tip vortices, as well as difference in dependencies of boundary layer thickness on Reynolds number for laminar and turbulent boundary layer. Provided asymptotic analysis gives two scaling laws with different values of power for high and low Reynolds numbers. Deduced values correspond well to known experimental data.     

Analysis on Link Between the Macroscopic and Microscopic Air–Water Properties in Self-Aerated Flows

Self-aeration in high-speed free surface flows occurs commonly and is of interest to ocean engineering, hydraulic engineering, and environmental engineering. For two-phase air–water flows, macroscopic air–water flow properties develop gradually, accompanied by the change of microscopic air–water structures. In this article, representational experimental studies on macroscopic and microscopic characteristics of self-aerated open-channel flows are summarized and compared. The isolated effect of the flow Reynolds number and air quantity on the differences in air count rate and chord size are analyzed and discussed. The results show that the characterized flow depth y₅₀, affected by the turbulence transfer, is a specific criterion to distinguish the interior air–water structure development. Two distinct linear trends of self-aeration are found, depending on the y₅₀/y₉₀ variation with a breaking point at C_mean = 0.50. The air count rate and size scale in self-aerated flows are affected by the air quantity of self-aerated flows, even with identical flow Reynolds numbers. Thus, a specific parameter is proposed to assess the air–water structures and a series of self-similarity relationships in self-aeration properties are obtained. The link between macroscopic and microscopic air–water properties results in significant scale effect on air–water structures in self-aerated flows.     

内孤立波浅化破碎过程斜坡沉积物孔压响应特征实验分析

观测资料显示内孤立波沿斜坡浅化过程对海底沉积物的作用犹如一台水中吸尘器,在破碎转换阶段达到最强,甚至会触发一系列地质活动,引发地质灾害。为界定此过程中沉积物的动力响应特征和影响因素,在大型重力式分层流水槽中模拟不同振幅内孤立波和不同类型沉积物斜坡连续作用过程,利用孔隙水压力采集系统实时记录孔隙水压力变化,对比分析不同水动力、坡度、沉积物类型情况下沉积物中超孔压变化特征。分析结果表明,内孤立波破碎过程,破波位置海床表层波压力和不同深度超孔隙水压力都存在相似的"U"型负压力变化过程;破碎波经过位置沉积物表现为和表面波压力正相关的孔压响应特征。破碎点沉积物中超孔压幅值随深度减小,约在6%波长深度位置减少到坡面压力的50%。超孔压幅值和内孤立波振幅、沉积物类型和斜坡度密切相关,坡度由0.071变化到0.160时,波压力幅值可增大至1.6倍。内孤立波振幅变化不影响不同类型海床土动力响应规律,只与超孔隙水压力值大小有关,内孤立波对海床的动力作用可认为弹性作用。     

Numerical study on scale effect of form factor

The scale effect of form factor is investigated via a numerical approach in this paper, where the turbulent ship flow is computed by solving the steady and incompressible Reynolds-averaged Navier-Stokes and continuity equations. A wall function approach is employed to bridge the near-wall and outer turbulent flow region. The numerical scheme based on a finite-volume formulation is applied to discretize the coupled governing equation. For the sake of numerical stability, accuracy and economy, an identical grid is employed to compute ship flow at different Reynolds number, where the grid is optimized for the medium Reynolds number of the investigated range. Four surface ships and two sub-bodies with notably different geometrical characteristics are chosen as the investigated cases, where double-model flow without appendages is computed. The calculated total resistance coefficient shows a decreasing tendency against Reynolds number among all studied hulls. Similar to the calculated total resistance coefficient, the calculated friction resistance coefficient decreases with the Reynolds number and varies relatively little for a given Reynolds number among different hulls. The viscous pressure resistance coefficient is less insensitive to the Reynolds number but apparently depends on hull form. Compared with the form factor calculation based on empirical friction lines, the flat-plate friction prediction based on CFD approach clearly gives smaller Re-dependent form factor, which should more realistically reflect the scale effect of form factor. The form factor exhibits a near linear and increasing dependence on Reynolds number. The numerical results show that the dependence of r_P on Reynolds number mainly governs the scale effect of form factor.     

Experimental investigation on the wave-induced pore pressure around shallowly embedded pipelines

A series of regular wave experiments have been done in a large-scale wave flume to investigate the wave-induced pore pressure around the submarine shallowly embedded pipelines.The model pipelines are buried in three kinds of soils,including gravel,sand and silt with different burial depth.The input waves change with height and period.The results show that the amplitudes of wave-induced pore pressure increase as the wave period increase,and decay from the surface to the bottom of seabed.Higher pore pressures are recorded at the pipeline top and the lower pore pressures at the bottom,especially in the sand seabed.The normalized pressure around pipeline decreases as the relative water depth,burial depth or scattering parameters increase.For the silt seabed,the wavelet transform has been successfully used to analyze the signals of wave-induced pore pressure,and the oscillatory and residual pore pressure can be extracted by wavelet analysis.Higher oscillatory pressures are recorded at the bottom and the lower pressures at the top of the pipeline.However,higher residual pressures are recorded at the top and the lower pressures at the bottom of the pipeline.