URANS investigation of the interaction between the free surface and a shallowly submerged underwater vehicle at steady drift

An axisymmetric underwater vehicle (UV) at a steady drift angle experiences the complex three-dimensional crossflow separation. This separation arises from the unfavorable circumferential pressure gradient developed from the windward side toward the leeward side. As is well known, the separated flow in the leeward side gives rise to the formation of a pair of vortices, which affects considerably the forces and moments acting on the UV. In this regard, the main purpose of the present study is to evaluate the role of the leeward vortical flow structure in the hydrodynamic behavior of a shallowly submerged UV at a moderate drift angle traveling beneath the free surface. Accordingly, the static drift tests are performed on the SUBOFF UV model using URANS equations coupled with a Reynolds stress turbulence model. The simulations are carried out in the commercial code STARCCM+ at a constant advance velocity based on Froude number equal to F_n = 0.512 over submergence depths and drift angles ranging from h = 1.1D to h = ∞ and from β = 0 to β = 18.11°, respectively. The validation of the numerical model is partially conducted by using the existing experimental data of the forces and moment acting on the totally submerged bare hull model. Significant interaction between the low-pressure region created by the leeward vortical flow structure and the free surface is observed. As a result of this interaction, the leeward vortical flow structure appears to be largely responsible for the behavior of the forces and moments exerted on a shallowly submerged UV at steady drift.     

Numerical investigation of the seakeeping behavior of a catamaran advancing in regular head waves

A numerical study was undertaken in order to assess the capability of an unsteady RANS code to predict the seakeeping characteristics of a high-speed multi-hull vessel in high sea states. Numerical analysis includes evaluation of ship motions, effects of wave steepness on ship response, catamaran natural frequency and added resistance in waves. Computations were performed for the DELFT 372 catamaran by the URANS solver CFDSHIP-Iowa V.4. The code was validated with encouraging results for high ship speeds (0.3≤Fn≤0.75) and high wave amplitudes (0.025≤Ak≤0.1). Comparison with strip theory solutions shows that the RANS method predicts ship motions with higher accuracy and allows the detection of nonlinear effects. Current computations evidence that heave peaks occur at resonance for all Fn, and reach the absolute maximum at Fn=0.75. Maximum pitch occurs at frequencies lower than resonance, for each speed, and absolute maximum occurs at medium Fn=0.6. Maximum added resistance, R_aw, was computed at Fn=0.45, which, interestingly, is near the catamaran Fn_coincidence. Overall, we found similar results as Simonsen et al. (2008) for KCS containership, though, herein, a multi-hull geometry and higher speeds were tested. Also, our results are useful to further evaluate the exciting forces and their correlation with f_e and λ/L_pp.     

Hydrodynamic assessment of planing hulls using overset grids

In conjunction with high performance computers, recent developments in computational science paved the path to more accurate representation of body motions inside fluids. Small motions inside the flow can be computationally approximated by using rigid body motion but it is incapable of accurately predicting the large motions of a planing vessel. The implementation of overset grid has made it possible to better approximate the complex fluid-structure interaction problem of the planing regime. The focus of this study was to evaluate the opportunity of using an overset grid system to numerically solve the flow around a planing hull and to understand the planing regime with this invaluable tool. It was shown in this study that the overset grid better captures the large motions of the planing hull at high Froude numbers. Then, the results obtained by overset grid were used to calculate the resistance components of a planing hull in a wide Froude number range. The resistance components were discussed with respect to values generated by Savitsky approach. Using the benefits that the computational science brings, the flow was visualized to explain some underlying physics relevant to the planing regime.     

Numerical and experimental study on seakeeping performance of a SWATH vehicle in head waves

The seakeeping characteristics of a Small Waterplane Area Twin Hull (SWATH) vehicle equipped with fixed stabilizing fins was investigated by experimental and numerical methods The calculation methods range from viscous CFD simulation based on an unsteady RANS approach to Boundary Element Method (BEM) based on Three Dimensional Translating-pulsating Source Green Function (3DTP). Responses of ship motions in head regular waves and nonlinear effects on motion responses with increasing wave amplitude were analyzed. Numerical simulations have been validated by comparisons with experimental tests. The results indicate that the heave and pitch transfer functions depict two peaks with the increase of wave length. Comparisons amongst experimental data and different numerical calculations illustrate that the RANS method predicts ship motions with higher accuracy and allows the detection of nonlinear effects. The heave and pitch transfer functions see a downward trend with the increasing wave amplitude in the resonant zone at low speed.     

Development of hybrid URANS-LES methods for flow simulation in the ship stern area

The paper presents the results of the application of a new hybrid URANS-LES method for the investigations of the ship wake behind the tanker KVLCC2. The switching between URANS and LES models is based on the ratio between the turbulence scale and the cell size of the mesh. Ship resistance, fields of the axial velocity and turbulent kinetic energy in the propeller plane are calculated and compared with measurements. Much attention is paid to the analysis of the unsteady velocities, their PDF distributions and spectra. Numerical analysis shows that the instantaneous velocities deviate substantially from their mean values which are usually used as the estimated velocities in modern engineering methodologies. The thrust variation in the unsteady wake is more than twice as large as that in the time averaged (frozen) wake. The results of the present study point out that the unsteadiness in the wake behind full ships can be very large and should be taken into account when propulsion and unsteady loadings are determined.