Numerical modeling of flow and scour below a pipeline in currents: Part I. Flow simulation

The performance of the standard k–ε, Wilcox high-Reynolds-number k–ω, Wilcox low-Reynolds-number k–ω and Smagorinsky's subgrid scale (SGS) turbulence models is examined against the flow around a circular cylinder 0.37 diameter above a rigid wall. The governing equations are solved using finite difference method in a non-orthogonal boundary-fitted curvilinear coordinate system. A mesh dependence study for the four turbulence models is carried out on computational meshes with different densities. In addition, the performance of the k–ω models with either wall function or no-slip boundary condition on the cylinder surface is examined on the finest mesh. It is found that the SGS model over-predicts the shedding of vortices from the cylinder and is sensitive to the computational mesh and the model constant C_s used. The standard k–ε and the Wilcox k–ω models predict the mean velocity field quite well but generally under-predict the velocity and hydrodynamic force oscillations using wall functions on the cylinder surface. It is also found that the Wilcox k–ω models with the no-slip boundary condition on the cylinder surface give better predictions on the shedding of vortices than their counterparts using the wall function boundary condition.     

RANS v2–f simulation of a swash event: Detailed flow structure

The flow structure of a swash event over a uniform slope is studied using a RANS-VOF numerical model coupled with a v²–f turbulence closure. The model is compared with experimental data of recent laboratory experiments. The ability of the turbulence modelling for simulating swash flow and the evolution of the computed bed shear stress during run-up and run-down are investigated. The agreement between numerical results and measured data, such as water depth, depth-averaged velocity and bed shear stress is very good during run-up. Main discrepancies are found during run-down. The paper also examines the aeration of the water layer in the swash flow, taking advantage of the PLIC method for computing the air–water interfaces. Air is continuously entrapped in the swash front and released at its rear during run-up. A detailed analysis indicates that the flow reversal is initiated near the bottom at the outer boundary of the swash zone and progresses landward. The study highlights the asymmetry between run-up and run-down. During run-up, the swash front propagation determines the turbulence properties and the bed shear stress profile on the beach, whereas the flow properties are more homogeneously distributed in the swash area during run-down.     

Internal flow structure of short wind waves

The internal flow structure of wind waves in a wind-wave tunnel was investigated on the bases of the measured vorticity distributions, streamline patterns, internal pressure fields, and stress distributions at the water surface for some waves in the field. In part I the experimental method and the internal vorticity structure relative to the individual wave crests are described. The measured vorticity distributions of distinct waves (waves with waveheight comparable with or larger than that of significant wavesH _1/3) in the field indicate that the surface vorticity layer is extraordinarily thickened near the crest, and the vorticity near the water surface shows a particularly large value below the crest. The flow near the crest of distinct waves is found to be in excess of the phase speed in a very thin surface layer, and the tangential stress distribution has a dominant peak near the crest. It is argued that the occurrence of the region of high vorticity in distinct waves is associated with the local generation of vorticity near the crest by tangential stress which attains a peak, under the presence of excess flow.     

Numerical modeling of flow and scour below a pipeline in currents: Part II. Scour simulation

A vertical two-dimensional (2D) numerical model for time dependent local scour below offshore pipelines subject to unidirectional steady flow is developed. The governing equations for the flow and sediment transport are solved by using finite difference method in a general curvilinear coordinate system. The performance of two turbulence models, the standard k–ɛ model and Smagorinsky subgrid scale (SGS) model, on modeling time dependent scour processes is examined. Both suspended load and bed load are considered in the scour model. The suspended-load model is verified against two channel sediment transport cases. The change of bed level is calculated from the continuity equation of total sediment transport. A new time marching scheme and a sand slide scheme are proposed for the scour calculation. It is found that the proposed time marching scheme and sand slide model work well for both clear-water and live-bed scour situations and the standard k–ɛ turbulence closure is more preferable than the SGS model in the 2D scour model developed in this study.     

Numerical flow analysis of single-stage ducted marine propulsor

The present work has solved 3D incompressible RANS equations on a rotating, non-orthogonal multi-blocked grid to efficiently analyze a ducted marine propulsor with rotor–stator interaction. To handle the interface boundary between a rotor and a stator, the sliding multi-block technique using the cubic spline interpolation and the bilinear interpolation was applied. To validate the present code, the flow of a single stage turbine flow was simulated. Time averaged pressure coefficients were compared with experiments and good agreements were obtained. After the code validation, the flowfield around a single-stage ducted marine propulsor having a single stage of rotor and stator was successfully simulated and the hydrodynamic performance coefficients were computed.     

Savonius型转轮三角阵列减流性能研究

为了助力海洋牧场减流防护工程, 研究Savonius型转轮阵列减流性能。作者建立Savonius型转轮三角阵列尾流场数值模型, 并通过水池实验验证准确性, 基于可靠数值模型探究转轮阵列尾涡减流机理, 研究三角阵列结构参数L_X、L_Y, 以及动力参数TSR、初始流速、旋向对整体减流性能的影响规律。结果表明,下游转轮产生的涡流呈现非对称分布, 并且产生更多涡流的转轮拥有更好的减流效果。另外, L_X为3D和L_Y为2D时减流性能最佳。最后对比发现, 在叶尖速比为0.9~1.1减流效果更好; 初始流速大小不影响减流效果; 下游转子对称分布时, 随着上游转子改变旋转方向, 减流效果出现明显差异。     

Acoustic monitoring of gas emissions from the seafloor. Part I: quantifying the volumetric flow of bubbles

Three decades of continuous ocean exploration have led us to identify subsurface fluid related processes as a key phenomenon in marine earth science research. The number of seep areas located on the seafloor has been constantly increasing with the use of multi-scale imagery techniques. Due to recent advances in transducer technology and computer processing, multibeam echosounders are now commonly used to detect submarine gas seeps escaping from the seafloor into the water column. A growing number of en-route surveys shows that sites of gas emissions escaping from the seafloor are much more numerous than previously thought. Estimating the temporal variability of the gas flow rate and volumes escaping from the seafloor has thus become a challenge of relevant interest which could be addressed by sea-floor continuous acoustic monitoring. Here, we investigate the feasibility of estimating the volumetric flow rates of gas emissions from horizontal backscattered acoustic signals. Different models based on the acoustic backscattering theory of bubbles are presented. The forward volume backscattering strength and the inversion volumetric flow rate solutions were validated with acoustic measurements from artificial gas flow rates generated in controlled sea-water tank experiments. A sensitivity analysis was carried out to investigate the behavior of the 120-kHz forward solution with respect to model input parameters (horizontal distance between transducer and bubble stream, bubble size distribution and ascent rate). The most sensitive parameter was found to be the distance of the bubble stream which can affect the volume backscattering strength by 20 dB within the horizontal range of 0–200 m. Results were used to derive the detection probability of a bubble stream for a given volume backscattering strength threshold according to different bubble flow rates and horizontal distance.     

Characteristic friction coefficient and scale effects in oscillatory porous flow

This work presents a simple method to evaluate the performance of a porous breakwater when it is impinged with normal incidence by a non-breaking monochromatic wave train. It is based on: 1) a potential flow model for wave interaction with permeable structures and 2) a set of experimental tests on a rectangular porous structure with uniform granular distribution. A characteristic friction diagram is obtained considering wave energy balance in a control volume, minimising the error between the numerical model and the experimental results for the wave transmission coefficient. Results show that, for large breakwater widths, the reflection process reaches a saturation regime before the waves exit the structure at a distance from the seaside between the interval 0.2 < x/L < 0.45. For larger breakwater widths, the reflection coefficient is almost constant (except for “resonant” conditions) and wave transmission decreases exponentially. Under such conditions, the wave propagation through the porous medium depends on the relative diameter D/L and the porosity of the material; the dependence on the relative breakwater width B/L and the ratio diameter wave height D/H is weak. This diagram intends to be useful for preliminary engineering studies of breakwater's efficiency and performance and as an adequate selection criteria of the experimental stone diameter to minimize scale effects in laboratory studies.     

The inception of sheet flow in oscillatory flow

A simple model is developed to study the inception of sheet flow in oscillatory flow based on the available experimental data. The inception of sheet flow in oscillatory flow is well defined by the simple model: A/d=KA²ω/ν+B, where A is the semi-excursion of wave orbital motion near the bed, d is the grain size, ω is the angular frequency, ν is the kinematic viscosity of water, and K and B are the coefficients and dependent on sediment properties only. The inception velocity of sheet flow derived from the model is shown to be the function of grain size d, oscillatory period T and specific sediment density s. For a given sediment, the inception velocity is found to increase sharply initially with T and then approach a constant at T>6.0 s. The present model is quite simple and gives good agreement with the available experimental data.     

Phase-lag effects in sheet flow transport

The inception of the sheet flow regime as well as the effects of the phase lag when the sheet flow regime is established were investigated for oscillatory flows and combined steady and oscillatory flows. A new criterion for the inception of sheet flow is proposed based on around 300 oscillatory flow cases from experiments. This criterion was introduced in the Camenen and Larson [Camenen, B., Larson, M., 2005. A bedload sediment transport formula for the nearshore. Estuarine, Coastal and Shelf Science 63, 249–260.] bed load formula in order to take into account phase-lag effects in the sheet flow regime. The modification of the Camenen and Larson formula significantly improves the overall agreement with data and yields a correct behavior in relation to some of the main governing parameters, which are the median grain size d₅₀, the orbital wave velocity U_w, and the wave period T_w. The calibration of the new formula was based on more than 200 experimental data values on the net sediment transport rate for a full wave cycle. A conceptual model was also proposed to estimate the ratio between sediment transport rate with and without phase lag, (r_pl = q_s,net / q_s,net,ϕ=0). This simple model provides accurate results and may be used together with any quasi-steady model for bed load transport.     

Effect of wave amplitude on turbulent flow in a wavy channel by direct numerical simulation

The present study numerically investigates the characteristics of three-dimensional turbulent flow in a wavy channel. For the purpose of a careful observation of the effect of the wave amplitude on the turbulent flow, numerical simulations are performed at a various range of the wave amplitude to wavelength ratio (0.01?α/λ?0.05), where the wavelength is fixed with the same value of the mean channel height (H). The immersed boundary method is used to handle the wavy surface in a rectangular grid system, using the finite volume method. The Reynolds number (Re=U_bH/ν) based on the bulk velocity (U_b) is fixed at 6760. The present computational results for a wavy surface are well compared with those of references. When α/λ=0.02, the small recirculating flow occurs near the trough at the instant, but the mean reverse flow is not observed. In the mean flow field, the reverse flow appears from α/λ=0.03 among the wave amplitude considered in this study. The domain of the mean reverse flow defined by the locations of separation and reattachment depends strongly on the wave amplitude. The pressure drag coefficient augments with increasing the wave amplitude. The friction drag coefficient shows the increase and decrease behavior according to the wave amplitude. The quantitative information about the flow variables such as the distribution of pressure and shear stress on the wavy surface is highlighted.     

Steady streaming around a circular cylinder in an oscillatory flow

Steady streaming induced by an oscillatory flow around a circular cylinder is investigated using a numerical method. Two-dimensional Reynolds-averaged Navier-Stokes equations are solved using a finite element method with a k-ω turbulent model closure. The range of the Keulegan-Carpenter (KC) number investigated is between 2 and 40, which is substantially higher than those reported in literature related to steady streaming to date. A constant value of Stokes number (β) of 196 is chosen in this study. The steady streaming structures and velocity distribution are analysed in detail. It is found that the characteristics of steady streaming are strongly related to the vortex shedding flow regimes.     

Modelling and measurement of sand transport processes over full-scale ripples in oscillatory flow

A new series of laboratory experiments was performed in the Aberdeen Oscillatory Flow Tunnel (AOFT) and the Large Oscillating Water Tunnel (LOWT) to investigate time-averaged suspended sand concentrations and transport rates over rippled beds in regular and irregular oscillatory flow. The wave-induced oscillatory near-bed flows were simulated at full-scale. Five series of experiments were carried out. During the two AOFT experimental series, ripple dimensions, ripple migration rates and net sand transport rates were measured under regular and irregular asymmetric flow for two different sand types. The three LOWT experimental series focussed on measurements of the ripple dimensions, ripple migration rates, time-averaged suspended sand concentrations and net sand transport rates under regular asymmetric and irregular weakly asymmetric flow for two different sand types. From analysis of new and other full-scale data, it is concluded that the lower part of the time- and bed-averaged concentration profile (up to two times the ripple height above the ripple crest level) has an exponential profile. A new reference concentration formula is proposed based on the formula of Bosman and Steetzel [Bosman, J.J., Steetzel, H.J., 1986. Time- and bed-averaged concentration under waves. Proc. 20th ICCE Taipei, ASCE, pp. 986–1000], which includes the grain-size influence. Furthermore, it is shown that the concentration decay length is strongly related to the ripple height and that the simple formula R_c = 1.27η gives good agreement with the data. A new transport model is proposed for the wave-related net transport over full-scale ripples based on a modified half wave cycle concept of Dibajnia and Watanabe [Dibajnia, M., Watanabe, A., 1992. Sheet flow under nonlinear waves and currents. Proc. 23rd ICCE Venice, ASCE, pp. 2015–2028; Dibajnia, M., Watanabe, A., 1996. A transport rate formula for mixed sands. Proc. 25th ICCE Orlando, ASCE, pp. 3791–3804]. The magnitudes of the half wave cycle transport contributions are related to the grain-related Shields parameter, the degree of wave asymmetry and a newly defined vortex suspension parameter P, which is the ratio between the ripple height and the median grain-size. The new model has been calibrated using transport data from the new regular flow experiments and has subsequently been validated using other data, including measurements from irregular flow experiments. The new model is seen to perform better overall than existing practical models for ripple regime net sand transport.     

An experimental investigation of the characteristics of a cylinder interacting with waves and flow in the dominant-load regime

We describe in this paper the experimental investigations of the interaction of a bottom-pivoted vertical cylinder with water waves and flow, to determine the dominant-load-regime map by application of response step functions and response RAO. A rigid circular cylindrical mass-damper-spring oscillator system is investigated in regular waves and uniform flow to determine the response characteristics in the frequency domain. Interaction with waves dominates in the high frequency range f* = f_osc/ω_v = 0.862–1.547, with magnitude in the range of 0.1 rad. On the other hand, interaction with flow dominates at lower frequency range, f* = 0.442–0.862, with magnitude in the range of 0.01 rad. These are caused by the non-overlap peak positions of the magnitude response in waves and flow due to the change in added mass of the cylinder moving in different types of fluid loads. The frequency f* = 0.862 is the point where the dominant factors are transferred. The location of separation points determines the pressure distribution to induce the added mass changed. Separation positions determine the magnitude response, but do not determine the configuration of response RAO. That allows to enhance or reduce the magnitude response of the cylinder by taking advantage of the dominant-load-regime map in the frequency domain.     

Wind tunnel experiments on a fully appended laminar flow submersible for oceanographic survey

This paper reports experimental measurements of the hydrodynamics forces acting on a low drag submersible. The model was manufactured using standard, cost-effective methods without undue attention paid to the roughness or waviness of the model's surface. It was seen, using flow visualisation techniques, that the laminar boundary layer extended over 70% of the vehicle length; furthermore, no separation of the flow was seen over the remaining 30% of the body. Hydrodynamic forces were measured, and a significant reduction in drag was found compared with standard turbulent submersible shapes. The low drag performance was quite resilient to changes in incidence, showing a large drag “bucket” for angles of incidence between ± 4°. With the addition of a set of cruciform fins at the rear of the body, incorporating control surfaces, it was found that for a centre of gravity position forward of 41.3% of the body length, the vehicle would be statically stable both stick-free and stick-fixed. Finally, using a pitot tube wake rake the velocity profile of the boundary layer at the propeller location was found to be that of an attached, turbulent boundary layer closely matching a power law. The importance of full interference measurement tests to allow for the support strut effects is clearly demonstrated.     

Distribution of mega fauna on sulfide edifices on the Eastern Lau Spreading Center and Valu Fa Ridge

Hydrothermal vent sulfide edifices contain some of the most extreme thermal and chemical conditions in which animals are able to live. As a result, sulfide edifices in the East Pacific Rise, Juan de Fuca Ridge, and Mid Atlantic Ridge vent systems often contain distinct faunal assemblages. In this study, we used high-resolution imagery and in-situ physico-chemical measurements within the context of a Geographic Information System (GIS) to examine community structure and niche differentiation of dominant fauna on sulfide edifices in the Eastern Lau Spreading Center (ELSC) and Valu Fa Ridge (VFR) in the Western Pacific Ocean. Our results show that ELSC and VFR sulfide edifices host two distinct types of communities. One type, that covers the majority of sulfide edifice faces, is overall very similar to nearby lava communities and biomass is dominated by the same chemoautotrophic symbiont-containing molluscs that dominate lava communities, namely the provannid gastropods Alviniconcha spp. and Ifremeria nautilei and the mytilid bivalve Bathymodiolus brevior. The spatial distribution of the dominant molluscs is often a variation of the pattern of concentric rings observed on lavas, with Alviniconcha spp. at the tops of edifices where exposure to vent flow is the highest, and I. nautilei and B. brevior below. Our physico-chemical measurements indicate that because of rapid dispersion of vent fluid, habitable area for symbiont-containing fauna is quite limited on sulfide edifices, and the realized niches of the mollusc groups are narrower on sulfide edifices than on lavas. We suggest that competition plays an important role in determining the realized distributions of the mollusc groups on edifices. The other habitat, present in small patches of presumably hot, new anhydrite, is avoided by the dominant symbiont-containing molluscs and inhabited by crabs, shrimp and polynoids that are likely more heat tolerant. The ratio of sulfide concentration to temperature anomaly of vent fluids was significantly different between sulfide edifice sites and lava sites in the southern vent fields but not in the northern vent fields. We suggest that this is due to increased sulfide consumption by a large microbial consortium associated with the more friable andesitic lava substrates in the south.     

Experimental Investigation into the performance of the Axial-Flow-Type Waterjet according to the Variation of Impeller Tip Clearance

The tip clearance inside the duct from the tip of the impeller is very important to the performance of waterjet systems, which fact has been proven in the pump field. The tip clearance is especially important on the model scale because it is very difficult in manufacture to keep the tip clearance constant and minimally small along the inside of the duct. In the present study, a flush-type waterjet propulsion unit (duct, impeller, stator, and nozzle) was designed for an amphibious tracked vehicle. Two impellers of different inner diameter were designed and manufactured in order to investigate the gap effect. Resistance and self-propulsion tests with a 1/5-scale model were conducted in PNU towing tank. The flow rate at the nozzle exit, the static pressure at the various sections along the duct and also the nozzle, the revolution of the impeller, and the torque, thrust, and towing forces at various advanced speeds were measured. Based on these measurements, the performance was analyzed according to the ITTC 96 standard analysis method. Based on this analysis method, the full-scale effective and delivered power of the tracked vehicle was estimated according to the variation of tip clearance.     

单体喷水推进船舶动力定位控制方法研究

面向一型喷水推进器为动力定位执行机构的单体喷水推进船舶,开展了船舶仿真模型建立、推力分配及进一步动力定位相关控制算法设计与仿真模拟。针对传统的定步长Kalman滤波提出了改进方案,同时考虑到单体喷水推进船舶的特殊性,对推力分配算法进行了优化,并基于荷载前馈的PD控制算法开展了不同工况下的仿真试验。试验结果表明,所提出的喷水推进推力分配算法能够有效控制一型单体喷水推进船舶实现定点、定艏向、移位、转艏等动力定位典型功能,并具备一定抗外部干扰能力。     

Numerical simulation of a partially buried pipeline in a permeable seabed subject to combined oscillatory flow and steady current

Hydrodynamic forces exerting on a pipeline partially buried in a permeable seabed subjected to combined oscillatory flow and steady current are investigated numerically. Two-dimensional Reynolds-Averaged Navier-Stokes equations with a k−ω turbulent model closure are solved to simulate the flow around the pipeline. The Laplace equation is solved to calculate the pore pressure below the seabed with the simulated seabed hydrodynamic pressure as boundary conditions. The numerical model is validated against the experimental data of a fully exposed pipeline resting on a plane boundary under various flow conditions. Then the flow with different embedment depths, steady current ratios and KC numbers is simulated. The amplitude of seepage velocity is much smaller than the amplitude of free stream velocity as expected. The normalized Morison inertia, drag and lift coefficients based on the corresponding force coefficients of a fully exposed pipeline are investigated. The normalized Morison force coefficients reduce almost linearly with the increase of embedment depth and that the KC only has minor effect on the normalized Morison coefficients. It is also found that the permeable seabed condition causes a slight increase on the inline force and has a little effect on the lift force, compared with corresponding conditions in an impermeable bed.