Investigation of drag crisis phenomenon using CFD methods

When fluid flow passes a cylinder, the drag crisis phenomenon occurs between the sub-critical and the super-critical Reynolds numbers. The focus of the present studies was on the numerical prediction of the drag crisis based on CFD methods. In this work, block structured meshes with refined grids near the cylinder surface and in the downstream were employed. Both 2D and 3D simulations were performed using various turbulence models, including the SST k − ω model, the k − ϵ model, the SST with LCTM, the DES model, and the LES model. In the convergence studies, the effects of the grid size, the time step, the first grid size and the aspect ratio (for 3D simulations) on the solutions were examined. The errors due to spatial and time discretizations were quantified according to a V&V procedure. Validation studies were carried out for various Reynolds numbers between Re = 6.31 × 10⁴ and 7.57 × 10⁵. The averaged drag force, the RMS of lift force and the Strouhal number were compared with experimental data. The studies indicated that standard 2D and 3D RANS methods were inadequate to capture the drag crisis phenomenon. The LES method however has the potential to address the problem.     

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.     

Large-eddy simulation of the flow past a circular cylinder at sub- to super-critical Reynolds numbers

Large-eddy simulation of turbulent flow past a circular cylinder at sub- to super-critical Reynolds numbers is performed using a high-fidelity orthogonal curvilinear grid solver. Verification studies investigate the effects of grid resolution, aspect ratio and convection scheme. Monotonic convergence is achieved in grid convergence studies. Validation studies use all available experimental benchmark data. Although the grids are relatively large and fine enough for sufficiently resolved turbulence near the cylinder, the grid uncertainties are large indicating the need for even finer grids. Large aspect ratio is required for sub-critical Reynolds number cases, whereas small aspect ratio is sufficient for critical and super-critical Reynolds number cases. All the experimental trends were predicted with reasonable accuracy, in consideration the large facility bias, age of most of the data, and differences between experimental and computational setup in particular free stream turbulence and roughness. The largest errors were for under prediction of turbulence separation.     

Force and flow characteristics of a circular cylinder with uniform surface roughness at subcritical Reynolds numbers

The main purpose of this study is to establish a better understanding of the relationship between drag reduction and surface roughness. Experiments were conducted to measure the force and flow characteristics of a circular cylinder with different types of artificial surface roughness over the range 6 × 10³ < Re < 8 × 10⁴ (Re is based on the cylinder diameter D). The roughness cylinder was formed by covering the exterior surface of the cylinder with uniformly distributed (1) sandpaper, (2) netting, and (3) dimples. The roughness coefficient ranged from k/D = 0.0028 to 0.025 (k is the roughness height). A detailed quantitative measurement of the flow field around the cylinder using Particle Imaging Velocimetry (PIV) was carried out. The hydrodynamic force coefficients (drag and lift) of the rough cylinders are compared against those of a smooth cylinder measured under the same flow conditions. It is found that certain configuration of surface roughness significantly reduces the mean drag coefficient of the cylinder, particularly at large Reynolds numbers. In addition, the root-mean-square (r.m.s.) lift coefficient of the rough cylinders is considerably lower than that of a smooth cylinder.     

海洋立管试验支座设计及其对涡激振动影响研究

针对大长细比顶张力立管内外流共同作用时涡激振动的试验研究,考虑端部铰接和固接两种边界条件,设计了可施加张力和内流的铰接支座和固接支座,并对两种支座进行数值模拟及受力分析,最后用于深水立管的涡激振动试验。立管模型材料采用6.2 m长铜管,外径20 mm,壁厚1.5 mm,分别支撑在铰接支座和固接支座上,施加不同的顶部张力、外流流速和内流流速。通过在立管表面粘贴光纤光栅应变计获得动态应变数据,分别从动态响应幅值、频率等方面对不同边界条件立管的试验数据进行对比分析,得出在顶张力、内流速相同的情况下,外流速相等时两端铰接立管的振动幅值大于两端固接立管,随着外流速的增加,两端铰接立管比两端固接立管率先发生锁振现象。     

On the planing of a flat plate at high Froude numbers in a two-dimensional case

We seek the solution of the planing of a flat plate at high Froude numbers by a perturbation procedure. The angle of attack of the plate is assumed to vary with the speed of the plate in the present study. A harmonic function K is introduced for the solution of the first-order disturbance potential which becomes the Green function in the limiting case when the Froude number tends to infinity. We get the solution of the first-order potential from Green's theorem applied to K and the first-order potential. Then we obtain the asymptotic solutions of the angle of attack α, lift L and drag D as follows:

where α1. Here W, L_W, and U are the weight of the plate per unit width, wetted length, and speed of the plate, respectively.     

Reynolds and Mass-Damping Effect on Prediction of the Peak Amplitude of A Freely Vibrating Cylinder

The Reynolds effect and mass-damping effect on the peak amplitude of a freely vibrating cylinder is studied by using forced oscillating data from Gopalkrishnan＇ s research in 1993, in which all experimental cases were carried out at a fixed Reynolds and the tested cylinder was recognized as a body that had no mass and damping. However, the Reynolds and roass-damping are the very important parameters for the peak amplitude of a freely vibrating cylinder. In the present study, a function F is introduced to connect the forced oscillation and free vibration. Firstly the peak amplitude AG^＊ can be obtained from the function F using forced oscillation data of Gopalkrishnan＇ s experimental at Re = 10^4, and then the Reynolds effect is taken into account in the function f（Re）, while the mass-damping effect is considered in the function K（ α ）, where a is the mass-damping ratio. So the peak amplitude of a freely vibrating cylinder can be predicted by the expression： A ^＊ = K（ α ）f（ Re ）AG^＊ . It is found that the peak transverse amplitudes predicted by the above equation agree very well with many recent experimental data under both high and low Reynolds conditions while roass-damping varies. Furthermore, it is seen that the Reynolds number does have a great effect on the peak amplitude of a freely vibrating cylinder. The present idea in this paper can be applied as an update in the empirical models that also use forced oscillation data to predict the vortex induced vibration （VIV） response of a long riser in the frequency domain.     

亚临界区雷诺数下圆柱绕流场电磁力控制数值研究

在亚临界区雷诺数下,采用脱体涡模拟方法对弱电解质中电磁力作用下圆柱绕流场及其升阻力特性进行了数值模拟与分析。结果表明,电磁力可以提高圆柱体边界层内的流体动能,延缓圆柱体近壁面流动分离,减弱绕流场中流向和展向大尺度涡的强度,减小圆柱体阻力及其升力脉动幅值;当电磁力作用参数大于某个临界值后,流动分离角消失,在圆柱体尾部产生射流现象,电磁力产生净推力作用,出现负阻力现象,而且升力脉动幅值显著减小且接近于零。     

Numerical modelling of short- and long-wave transformation on a barred beach

This work aims to demonstrate an advancement towards the integrated modelling of surf zone hydrodynamics by means of a VOF-type numerical model (COBRAS-UC) based on the Reynolds-Averaged Navier–Stokes equations. In this paper, the numerical model is adapted and validated for the study of nearshore processes on a mildly-sloping beach. The model prediction of wave energy transformation and higher order statistics (skewness and asymmetry) are in good agreement with detailed laboratory observations from a barred beach [Boers, M. (1996). “Simulation of a surf zone with a barred beach; Report 1: Wave heights and wave breaking”. Tech. Rep.96-5, Comm. on Hydrol. and Geol. Eng., Dept. of Civil Engineering, Delft University of Technology]. Moreover, the numerical model allows us to study the low-frequency motions inside the surf zone. It is found that in order to achieve a satisfactory simulation of both short- and long-wave transformation, the numerical model must achieve: (i) the simultaneous second-order wave generation and absorption, (ii) the energy transfer between triad of components, (iii) the short- and long-wave energy dissipation inside the surf zone, and (iv) the wave reflection at the shoreline. Comparisons between numerical and experimental results demonstrate the model capability to satisfactorily simulate all the aforementioned processes.     

Numerical study of effect of fluke inclination on behavior of drag anchor in clay with linearly increasing shear strength under unidirectional and combined loading

Drag anchor is a widely used anchor type in offshore engineering for the mooring system. The prediction of the anchor trajectory installation and the final position is important for anchor selection in design. The existing method using yield envelope method for trajectory prediction ignored the shallow anchor behavior but applied the deep yield envelope from a deeply embedded horizontal fluke in uniform clay for the whole drag-in installation process. However, the anchor fluke embedment depth and inclination angle change continually during installation in clay with linearly increasing shear strength soil profile in practice. Studies on the effect of fluke inclination angle on the anchor behavior in clay with such non-uniform soil profile under unidirectional and combined loading are important and necessary for the improvement of the yield envelope method to ensure a reasonable prediction. The current 2D finite element studies investigate the anchor behavior for inclined fluke in clay with linearly increasing shear strength under unidirectional vertical, horizontal and rotational loadings first. Then the effects of the fluke inclination angle, soil non-homogeneity and embedment depth ratio on the shallow yield envelopes are investigated. It is found that the effect of fluke inclination angle on the vertical capacity factors for anchor in clay with non-uniform and uniform soil profile is largely different. The resultant large impact on the yield envelopes shown here illustrates the importance of considering the fluke inclination angle and soil non-homogeneity in the prediction of anchor trajectory using yield envelope method.     

Practical evaluation of sinkage and trim effects on the drag of a common generic freely floating monohull ship

A practical method to account for the influence of sinkage and trim on the drag of a freely floating (free to sink and trim) common monohull ship at a Froude number F ≤ 0.45 is considered. The sinkage and the trim are estimated via two alternative simple methods, considered previously. The drag is also estimated in a simple way, based on the classical Froude decomposition into viscous and wave components. Specifically, well-known semiempirical expressions for the friction drag, the viscous pressure drag and the drag due to hull roughness are used, and the wave drag is evaluated via a practical linear potential flow method. This simple approach can be used for ship models as well as full-scale ships with smooth or rough hull surfaces, and is well suited for early ship design and optimization. The method considered here to determine the sinkage and the trim, and their influence on the drag, yields theoretical predictions of the drag of the Wigley, S60 and DTMB5415 hulls that are much closer to experimental measurements than the corresponding predictions for the hull surfaces of the ships in equilibrium position at rest. These numerical results suggest that sinkage and trim effects, significant at Froude numbers 0.25 < F, on the drag of a typical freely floating monohull ship can be realistically accounted for in a practical manner that only requires simple potential flow computations without iterative computations for a sequence of hull positions.     

Impact of a scallop drag on the sediment chemistry, microbiota, and faunal assemblages of a shallow subtidal marine benthic community

Scallops are usually obtained by means of a heavy metal dredge that is pulled over the sea bottom. Most studies of the impact of this gear type have shown that larger invertebrates, in particular, are severely disturbed. These studies, however, have been conducted on coarse sediments, ranging from sands to cobble, and have dealt only with faunal changes. In this paper the impact of a New England type scallop dredge on the fauna and sedimentary nutritional characteristics of a silty sand community is detailed. The site, in the Damariscotta River, Maine, USA, was sampled during the fall and winter prior to, then again immediately following, the dragging event, and twice more over the ensuing six months. Loss of surficial sediment, lowered food quality of the sediment (as measured by microbial populations, enzyme hydrolysable amino acids, and chlorophyll a), and changes in the faunal composition of the dragged site were observed. While some taxa returned to the drag track relatively quickly, others such as the cumaceans, phoxocephalid and photid amphipods, and nephtyid polychaetes, were not seen in abundances comparable to those of the adjacent undragged site until the food quality also recovered.     

Effect of peak perforation on flow past a conic cylinder at <Emphasis Type="Italic">Re</Emphasis>=100: drag,lift and Strouhal number

A flow past a circular-section cylinder with a perforated conic shroud, in which the perforation is located at the peak of the conic disturbance as the shroud installed on the cylinder and uniformly distributed with several circular holes, is numerically simulated at a Reynolds number of 100. Two factors in the perforation are taken into account, i.e. the attack angle relative to the direction of incoming flow and diameter of holes. The effect of such perforation on the drag, lift and vortex-shedding frequency is mainly investigated. Results have shown that variation of the attack angle has a little effect, especially on the drag and vortex-shedding frequency, except in certain cases due to the varied vortex-shedding patterns in the near wake. The increasing hole diameter still exhibits a little effect on the drag and frequency of vortex shedding, but really reduces the lift, in particular at larger wavelength, such as the lift reduction reaching almost 66%–68% after introducing the perforation.     

高效几丁质降解菌CB101的分离、鉴定及其几丁质酶系的研究

几丁质是海洋中居首位的多糖物质,它是虾、蟹壳的主要成分,同时广泛存在于昆虫外骨骼及真菌的细胞壁中.许多生物包括植物、细菌、真菌和一些无脊椎动物与脊椎动物都能够产生几丁质酶降解几丁质,其中细菌利用几丁质酶降解几丁质作为生活所需碳源和氮源.从生态角度来看,几丁质酶在自然界几丁质循环过程中扮演了重要角色,作为对环境的适应,多种微生物都带有完整的几丁质酶系.     

Spartina anglica eradication experiment and in situ monitoring assess structuring strength of habitat complexity on marine macrofauna at high tidal level

The cordgrass Spartina anglica is an introduced species that tends to invade sheltered sand and mudflats, at the upper low marsh level. In Arcachon Bay, a lagoon in South West of France, the cordgrass can also replace Zostera noltii beds. The consequence of cordgrass presence on macrobenthic fauna was estimated and compared to adjacent habitats (bare sands, Z. noltii sea grass beds) during one year. The communities of the three habitats were characterised by low species richness, low abundance and biomass (when Hydrobia ulvae, 90% of abundance, is not considered) and high seasonal stability. The infaunal assemblages were particularly homogeneous between habitats without any characteristic species.Cordgrass eradication experiments were performed and zoobenthic recolonisation was observed the following year. Modifications in benthic fauna were observed on epifauna only.These results highlight the limited structuring effect of habitat heterogeneity at high tidal levels and in soft-bottom sediments where desiccation becomes the dominant factor determining infauna community structure.     

Effects of a high rise building on wind flow and beach characteristics at Atlantic City, NJ, USA

This study assesses the effects of a high rise building on the beach and dune at an intensively developed shoreline through field investigation of geomorphological and sedimentological characteristics and wind tunnel investigation of wind speed and direction. Results indicate that wind speed is increased close to the building, where wind is often deflected to a more onshore direction. Obliquely offshore northwesterly winds are the most dramatically altered and cause onshore flows (with pronounced upward flows) in the landward portion of the beach near the building. Onshore deflections in wind direction cause besanding of the boardwalk and create highly localized scour and deposition zones. A dune on the backbeach breaks up zones of sediment transport and scour and diminishes the likelihood that scour areas will cause increased hazard to landward buildings as a result of flooding. Generic large-scale studies are required to develop criteria for communities to use to specify optimum dimensions, locations and shapes of landforms that are modified or created by high rise structures.     

Analytical and experimental investigation on transverse vibrations of solid, circular and annular plates carrying a concentrated mass at an arbitrary position with marine applications

This investigation arose from the practical necessity of placing a centrifugal pump rigidly attached to a thin, circular cover plate of a water tank in a medium size ocean vessel. Due to lack of space, it was necessary to locate the system off—center of the circular configuration. It was considered necessary to calculate the fundamental frequency of the coupled system.The first part of the present study is concerned with the determination of the fundamental frequency of vibration of a circular plate carrying a concentrated mass at an arbitrary position, using a variational approach.Numerical results are obtained for the stated problem for several combinations of the intervening geometric and mechanical parameters.An experimental investigation is also performed in the case of clamped plates.Based on the results for solid circular plates, the fundamental frequency of annular plates with a free inner edge and a concentrated mass is also obtained.Circular plates are fundamental structural elements in ocean engineering applications: from off-shore platforms to underwater acoustic transducers. In a great variety of circumstances, they must carry operational systems in an eccentric fashion. Since the dynamic performance is always of interest, one must know at least some of the basic dynamic parameters.