Drag force on a circular cylinder in unseparated laminar high Reynolds number nonlinear random oscillatory flow

The stochastic properties of the drag force maxima on a circular cylinder subjected to nonlinear random waves are investigated. Unseparated laminar high Reynolds number flow is considered. A simplified approach based on second order Stokes waves is presented, including the sum-frequency effect only. It is demonstrated how a drag force formula valid for regular linear waves can be used to find the cumulative distribution function of individual drag force maxima for nonlinear irregular waves. Here the [Wang, 1968] drag force coefficient is used.     

Flow around a partly buried cylinder in a steady current

The flow around a cylinder, surface mounted or partly buried has been examined at Reynolds number of 1.3×10⁴ and 2.6×10⁴ by using the flow visualization method in the case of steady current. The lengths of the upstream and downstream separation regions were determined and the characteristics of the flow area were discussed for the cylinder having different burial depth to the diameter ratio (G/D). For comparsion, some numerical results obtained from the application of the FLUENT, computational fluid dynamics (CFD) software package have also been presented. Both experimental and numerical results indicated that the lengths of the separation regions near the upstream and downstream of the cylinder decreased with the increasing burial ratio (G/D=0, 0.20, 0.50). The results have significant implications both for the illustration of the potential for using CFD software and highlighting the need for data of two dimensional flow around the buried cylinder.     

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.     

Propulsive efficiency and structural response of a sandwich composite propeller

Studying the sandwich composite propeller (SCMP) is of great significance since the sandwich structure is lightweight and possesses high strength. This study proposes and verifies a fluid–structure interaction (FSI) method for a 3D underwater sandwich composite structure to calculate the performance of the propeller. The Reynolds-averaged Navier–Stokes formula-based computational fluid dynamics is adopted to solve for propeller loads, whereas the finite element method (FEM) is applied to solve for propeller deformations. ANSYS Workbench’s system coupling is utilized to deliver the loads and deformations in the FSI. The paper also compares the propulsive performance and structural response of the SCMP and conventional composite propeller (CMP). The impact of the structural form and core material on the SCMP is explored. The results show that the weight reduction effect of the SCMP is better than that of the CMP, the propulsive efficiency of the SCMP is higher at low advance coefficients and lower at high advance coefficients, and the maximum pitch angles of the SCMP decrease at all conditions, unlike the case for the CMP. Moreover, the thinner the facing of the SCMP, the greater the influence of the higher twist–deformation ratio of the resulting structural form on the intrinsic frequency.     

Experimental study of forces on a cylinder in oscillatory flow with a cross current

Forces on a circular cylinder have been measured with nominally two-dimensional current and oscillatory flow at right angles. Previous results for purely oscillatory flow defined by a Keulegan-Carpenter number, Kc, have been extended for reduced velocities, Vr, in the range 3–10. For Kc<7 modification of the Karman street by oscillation is complex and locking-on has a dominant influence. For Kc > 7 simply adding forces due to the current and oscillation as though in isolation generally gives conservative results. The ‘current’ drag shows considerable variation and can even be negative. The Morison fit to the in-line force is generally less satisfactory when there is a current and can be wholly inadequate.     

Flow around circular and noncircular cylinders by a low-order panel method

Time-dependent cross-flow was studied around cylinders with circular and noncircular cross-sections. The numerical approach for the analysis was a low-order panel method based on constant source and dipole values along each panel. The method was previously used successfully for several applications, such as calculation of the added mass and damping coefficients. In simulating the viscous time-dependent flow around the cylinder, the time-dependent wake feature of the code was used. For the circular and D-cylinders, the results agreed well with the experiments. Suggestions for improving the results for T-cylinders with angle of attack are included.     

Flow induced vibrations of a sharp edge square cylinder in the wake of a circular cylinder

The response of an oscillating circular cylinder at the wake of an upstream fixed circular cylinder was classified by different researchers as galloping, wake induced galloping or wake induced vibration. Furthermore it is already known that a sharp edge square cylinder would undergo galloping if it is subjected to uniform flow. In this study the influence of the wake of a fixed circular cylinder on the response of a downstream square cylinder at different spacing ratios (S/D = 4, 8, 11) is experimentally investigated. The subject appears not to have received previous attention. The lateral displacements, lift forces and the pressure data from gauges mounted in the wake of the oscillating cylinder are recorded and analyzed. The single degree of freedom vibrating system has a low mass-damping parameter and the Reynolds number ranges from 7.7 × 10² to 3.7 × 10⁴.In contrast to that for two circular cylinders in tandem arrangement, the freely mounted downstream square cylinder displays a VIV type of response at all spacing ratios tested. There is no sign of galloping or wake induced galloping with the square cylinder. With increase at the spacing ratio the cross-flow oscillations decrease. It is shown that the vortices arriving from the upstream fixed circular cylinder play a major role on the shedding mechanism behind the downstream square cylinder and cause the square cylinder to shed vortices with frequencies above Strouhal frequency of the fixed square cylinder (St = 0.13). The VIV type of oscillations in the downstream square cylinder is most probably caused by the vortices newly generated behind the square cylinder.     

Experimental and numerical investigation of local scour around a submerged vertical circular cylinder in steady currents

Local scour around a submerged vertical circular cylinder in steady currents was studied both experimentally and numerically. The physical experiments were conducted for two different cylinder diameters with a range of cylinder height-to-diameter ratios. Transient scour depth at the stagnation point (upstream edge) of the cylinder was measured using the so-called conductivity scour probes. Three-dimensional (3D) seabed topography around each model cylinder was measured using a laser profiler. The effect of the height-to-diameter ratio on the scour depth was investigated. The experimental results show that the scour depth at the stagnation point is independent on cylinder height-to-diameter ratio when the later is smaller than 2. The increase rate of equilibrium scour depth with cylinder height increases with an increase in Shields parameter.     

Viscosity and nonlinearity effects on the forces and waves generated by a floating twin hull under heave oscillation

Nonlinear hydrodynamics of a twin rectangular hull under heave oscillation is analyzed using numerical methods. Two-dimensional nonlinear time-domain solutions to both inviscid and viscous problems are obtained and the results are compared with linear, inviscid frequency-domain results obtained in [26] to quantify nonlinear and viscous effects. Finite-difference methods based on boundary-fitted coordinates are used for solving the governing equations in the time domain [2]. A primitive-variables based projection method [6] is used for the viscous analysis and a mixed Eulerian–Lagrangian formulation [11] for inviscid analysis. The algorithms are validated and the order of accuracy determined by comparing the results obtained from the present algorithm with the experimental results of Vugt [22] for a heaving rectangle in the free surface. The present study on the twin-hull hydrodynamics shows that at large and non-resonant regular frequencies, and small amplitude of body oscillation, the fluid viscosity does not significantly affect the wave motion and the radiation forces. At low frequencies however the viscosity effect is found to be significant even for small amplitude of body oscillation. In particular, the hydrodynamic force obtained from the nonlinear viscous analysis is found to be closer to the linear inviscid force than the nonlinear inviscid force to the linear inviscid force, the reason for which is attributed to the wave dampening effect of viscosity. Since the wave lengths generated at smaller frequencies of oscillation are longer and therefore the waves could have a more significant effect on the dynamic pressure on the bottom of the hulls which contribute to the heave force, the correlation between the heave force and the wave elevation is found to be larger at smaller frequencies. Because of nonlinearity, the wave radiation and wave damping force remained nonzero even at and around the resonant frequencies – with the resonant frequencies as determined in [26] using linear potential flow theory. As to be expected, the nonlinear effect on the wave force is found to be significant at all frequencies for large amplitude of oscillation compared to the hull draft. The effect of viscosity on the force, by flow separation, is also found to be significant for large amplitude of body oscillation.     

Behavior of vortex-induced vibration of a circular cylinder near a deformable wall with two degrees of freedom in steady flow

The behavior of vortex-induced vibration of a two-degree-of-freedom cylinder near a deformable wall in steady flow is investigated experimentally.The typical phenomenon of the two-degree-of-freedom cylinder’s VIV is discussed.The influences of initial gap between the cylinder and the wall on the dynamic responses of the cylinder are analyzed.The comparison is made about dynamic responses of the cylinder with one and two degrees of freedom.Experimental results show that the vibration of the cylinder near a deformable wall with a small value of initial gap-to-diameter ratios can generally be divided into two phases.The initial gap-to-diameter ratios have a noticeable influence on the occurrence of transverse vibration.The transverse maximum amplitude of the cylinder with two degrees of freedom is larger than that of the cylinder with one degree of freedom under the condition with the same values of other parameters.However,the vibration frequency of the cylinder for the two degrees of freedom case is smaller than that for the one degree of freedom case at the same value of Vr number.     

水流作用下渔网养殖空间变化的计算方法

渔网的动力分析是养殖网箱设计和优化的基础。通过把一系列杆相互铰接在一起,建立网的模型,模拟网在来流作用下的变形。作用在杆上的水动力可以通过Morison公式来计算,采用有限单元的离散方法,并利用网目合并的方法来减少总的自由度数目。模型的计算结果同试验成果做了对比,计算与试验在柔性网的变形方面,结果吻合良好。     

Hydroelastic analysis of floating structures with liquid tanks and comparison with experimental tests

In this study, we develop a numerical method for a 3D linear hydroelastic analysis of floating structures with liquid tanks subjected to surface regular water waves and compare the numerical results with experimental tests. Considering direct couplings among structural motion, sloshing, and water waves, a mathematical formulation and a numerical method extended from a recent work [1] are developed. The finite element method is employed for the floating structure and internal fluid in tanks, and the boundary element method is used for the external fluid. The resulting formulation completely incorporates all the interaction terms including hydrostatic stiffness and the irregular frequency effect is removed by introducing the extended boundary integral equations. Through various numerical tests, we verify the proposed numerical method. We also performed 3D hydroelastic experimental tests of a floating production unit (FPU) model in an ocean basin. The measured dynamic motions are compared with the numerical results obtained using the proposed method.     

Numerical study of nonlinear freak wave impact underneath a fixed horizontal deck in 2-D space

Freak waves are extreme and unexpected surface waves with huge wave heights that may lead to severe damage to ships and offshore structures. However, few researches have been conducted to investigate the impact underneath fixed horizontal decks caused by freak waves. To study these phenomena, a 2-D numerical wave tank is built in which nonlinear freak waves based on the Peregrine breather solution are generated. As a validation, a regular-wave-induced underneath impact is simulated and compared to the existing experimental measurements. Then the nonlinear freak-wave-induced impact is investigate with different values of deck clearance above the mean free surface. In addition, a comparative simulation of a “large” regular wave based on the 2nd-order Stokes wave theory with the same crest height and wave length of the nonlinear freak wave is carried out to reveal the unique features of the nonlinear freak-wave-induced impact. By applying a fluid–structure interaction (FSI) algorithm in which the bottom deck and front side wall are simplified as Euler beams in 2-D and discretized by the finite element method (FEM), the hydroelastic effects are considered during the impact event. The vertical force acting underneath the bottom deck, the transversal force acting on the front side wall, the structural displacements of the elastic deck and wall are analyzed and discussed respectively, from which meaningful conclusions are drawn.     

Wave-induced mean flows in vertical rubble mound structures

Waves impinging on rubble mound breakwaters and seawalls induce a mean flow within the breakwater, analogous to the so-called undertow within the surf zone. Here, using a plane wave approximation (kh<1.5), a second-order problem is solved for an idealized breakwater with a rectangular cross-section to show the origin and the nature of the mean flow within the porous structure. The mean flow is expressed in terms of a mean stream function analytically derived, obtained based on the mass flux balance between the incident, reflected and transmitted waves. Furthermore, the evolution of other second-order magnitudes such as mean water level and mass flux is analyzed under different incident wave conditions, structure geometry and porous material characteristics. Results show that the evolution of the different mean quantities is controlled mainly by reflection and consequently depends highly on structure geometry and porous material characteristics. Furthermore, it is shown that the return flow is stronger with increasing mass flux decay. Some qualitative experiments to show the described mechanism are also presented.     

Multi-resolution MPS for incompressible fluid-elastic structure interactions in ocean engineering

The paper presents a multi-resolution MPS (Moving Particle Semi-implicit)-based FSI (Fluid-Structure Interaction) solver for efficient and accurate simulations of incompressible fluid flows interacting with elastic structures. The fluid model is founded on the projection-based MPS solution of continuity and Navier-Stokes equations. The structure model is set based on MPS-based discretization of linear and angular momenta corresponding to an isotropic elastic solid. Fluid-structure coupling is conducted in a mathematically-physically consistent manner along with implementation of a multi-resolution scheme comprising of common radius of influence, revised weight function, revised number density and potential number density concepts to enhance i) consistency of particle-based discretizations, ii) imposition of boundary conditions and iii) volume conservation at fluid-structure interface. A set of previously developed enhanced schemes are also adopted for the fluid model. The robustness and efficiency of proposed Enhanced Multi-resolution MPS-MPS FSI solver are investigated through a set of ocean engineering-related benchmark tests. To the best knowledge of authors, this study presents the first multi-resolution particle method for FSI corresponding to incompressible fluid and elastic structures.     

Comparison of flow patterns in the near wake of a circular cylinder and a square cylinder placed near a plane wall

The flow patterns in the near wake of a cylinder (either circular or square in shape, D=25 mm) placed in the proximity of a fully developed turbulent boundary layer (thickness δ=0.4D) are investigated experimentally using particle image velocimetry (PIV). The effects of changing the gap height (S) between the cylinder bottom and the wall surface, over the gap ratio range S/D=0.1–1.0, have been investigated. The results show that both the ensemble-averaged and instantaneous flow fields are strongly dependent on S/D. The flow patterns for the two types of cylinders share many similarities with respect to the change in S/D, such as the reduced recirculation length and increased velocity fluctuation in the near wake with increasing S/D, as well as the trend of suppression of vortex shedding at small S/D and onset of vortex shedding at large S/D. However, developments of the shear layers, in terms of wake width, flow curvature, etc., are considerably different for these two types of cylinders. In general, the wake development and momentum exchange for the square cylinder are slower those for the circular cylinder at the same gap ratio. Correspondingly, it is shown that the periodic vortex shedding is delayed and weakened in the case of square cylinder, as compared to that of the circular cylinder at the same S/D.     

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.     

Hybrid three-dimensional finite-difference and finite-element analysis of seismic wave induced fluid–structure interaction of a vertical cylinder

The two-dimensional finite-difference scheme has been extended to three dimensions to solve nonlinear hydrodynamic pressures and structural responses of a deformable, vertical and circular surface-piercing offshore cylinder during earthquakes. A complete three-dimensional analysis has been made with both the three-dimensional equations of motion and the simultaneous action of three components of ground acceleration included in the analysis. Not only the magnitude but also the direction of the acting ground motion can be varied with time. The dynamic response of a cylinder is approximated by the displacements in the fundamental modes of vibration. A comparison of the dynamic displacement of the cylinder with and without surrounding sea water has been made. The flexibility of the offshore cylinder can significantly increase the hydrodynamic pressures acting on cylinder faces, that is, the fluid-structure interaction is necessary in offshore cylinder analysis. Although the hydrodynamic pressure induced by the vertical ground acceleration of the El Centro 1979 earthquake is significant, the calculated structural dynamic response of a cylinder is very small and the corresponding resultant hydrodynamic force is almost nil. The hydrodynamic force induced by two-horizontal ground acceleration is about the same as that by three simultaneous components of ground acceleration. For a solid and stubbier circular cylinder, the vertical component of ground acceleration may be neglected.     

Interaction mechanisms among waves,currents and a submerged plate

Wave-induced loads on a submerged plate, representative of submerged breakwater, coastal-bridge deck and a certain type of wave energy converter, in a uniform current are investigated in this study using fully nonlinear numerical wave tanks (NWTs) based on potential flow theory. The coupling effect of wave and current is explored, and the underlying interaction mechanisms of the hydrodynamic forces are described. The presence of a background current modifies the frequency dispersion. It produces changes of the water-surface elevation, and also has an effect on wave-induced loads. Depending on the nonlinearity, higher harmonic wave components are generated above the submerged plate. These contribute to the wave forces. It is found that the horizontal and the vertical force, hence the moment, are affected in the opposite way by the currents. The Doppler shifted effect dominates the vertical force and the moment on the plate. Whereas, the Doppler shifted effect and the generation of higher wave harmonics play opposite roles on the horizontal forces. The contribution of 2^nd order harmonics is found to be up to 30% of the linear component. The current-induced drag force, represented by the advection term ρU∂φ/∂x in the pressure equation, is found to lead to a decrease in the moment for the most range of wavelengths considered, and an increase in the moment for a small range of longer waves.