Numerical Prediction of Impact-Related Wave Loads on Ships

2006 ◽  
Vol 129 (1) ◽  
pp. 39-47 ◽  
Author(s):  
Thomas E. Schellin ◽  
Ould el Moctar

We present a numerical procedure to predict impact-related wave-induced (slamming) loads on ships. The procedure was applied to predict slamming loads on two ships that feature a flared bow with a pronounced bulb, hull shapes typical of modern offshore supply vessels. The procedure used a chain of seakeeping codes. First, a linear Green function panel code computed ship responses in unit amplitude regular waves. Ship speed, wave frequency, and wave heading were systematically varied to cover all possible combinations likely to cause slamming. Regular design waves were selected on the basis of maximum magnitudes of relative normal velocity between ship critical areas and wave, averaged over the critical areas. Second, a nonlinear strip theory seakeeping code determined ship motions under design wave conditions, thereby accounting for the nonlinear pressure distribution up to the wave contour and the frequency dependence of the radiation forces (memory effect). Third, these nonlinearly computed ship motions constituted part of the input for a Reynolds-averaged Navier–Stokes equations code that was used to obtain slamming loads. Favorable comparison with available model test data validated the procedure and demonstrated its capability to predict slamming loads suitable for design of ship structures.

Author(s):  
Thomas E. Schellin ◽  
Ould El Moctar

We present a numerical procedure to predict impact-related wave-induced (slamming) loads on ships. The procedure was applied to predict slamming loads on two ships that feature a flared bow with a pronounced bulb, hull shapes typical of modern offshore supply vessels. The procedure used a chain of seakeeping codes. First, a linear Green function panel code computed ship responses in unit amplitude regular waves. Wave frequency and wave heading were systematically varied to cover all possible combinations likely to cause slamming. Regular design waves were selected on the basis of maximum magnitudes of relative normal velocity between ship critical areas and wave, averaged over the critical areas. Second, a nonlinear strip theory seakeeping code determined ship motions under design wave conditions, thereby accounting for the ship’s forward speed, the swell-up of water in finite amplitude waves, as well as the ship’s wake that influences the wave elevation around the ship. Third, these nonlinearly computed ship motions constituted part of the input for a Reynolds-averaged Navier-Stokes equations (RANSE) code that was used to obtain slamming loads. Favourable comparison with available model test data validated the procedure and demonstrated its capability to predict slamming loads suitable for design of ship structures.


2018 ◽  
Vol Vol 160 (A2) ◽  
Author(s):  
M Pawłowski

Application of strip theory for the prediction of ship motions in waves relies on sectional hydrodynamic coefficients; i.e. the added mass and damping coefficients. These coefficients apply to linearised problems and are normally computed for inviscid fluids. It is possible to account for viscosity but this cannot be done by the RANS equations, as in linear problems there is no room for turbulence. The hydrodynamic coefficients can include the effect of viscosity but this can be done rightly through the classic Navier–Stokes equations for laminar (non-turbulent) flows. For solving these equations commercial RANS software can be used, assuming no Reynolds stresses.


Author(s):  
M Pawłowski

Application of strip theory for the prediction of ship motions in waves relies on sectional hydrodynamic coefficients; i.e. the added mass and damping coefficients. These coefficients apply to linearised problems and are normally computed for inviscid fluids. It is possible to account for viscosity but this cannot be done by the RANS equations, as in linear problems there is no room for turbulence. The hydrodynamic coefficients can include the effect of viscosity but this can be done rightly through the classic Navier–Stokes equations for laminar (non-turbulent) flows. For solving these equations commercial RANS software can be used, assuming no Reynolds stresses.


Author(s):  
Bogdan Iwanowski ◽  
Henrik Grigorian ◽  
Ingar Scherf

Subsidence of the Ekofisk platforms creates several operational challenges. For safety of the platforms, it is of great importance to find the wave impact loads acting on the platforms’ decks. The paper describes how such loads can be computed. Three theoretical wave models are discussed in the paper: the Airy wave, Airy wave modified through Wheeler stretching and the 5th order non-linear Stokes wave. The wave loads for these wave models are computed by various methods. The method based on momentum displacement approach and Morison-type equation developed by Dr. Kaplan is used as a reference point. The loads are also computed through a solution of complete Navier-Stokes equations, with the Volume of Fluid (VOF) method used to trace motion of the fluid’s free surface. Results of different wave models and different computational methods are compared and discussed.


2003 ◽  
pp. 55-82
Author(s):  
M. Despotovic ◽  
Milun Babic ◽  
D. Milovanovic ◽  
Vanja Sustersic

This paper describes a three-dimensional compressible Navier-Stokes code, which has been developed for analysis of turbocompressor blade rows and other internal flows. Despite numerous numerical techniques and statement that Computational Fluid Dynamics has reached state of the art, issues related to successful simulations represent valuable database of how particular tech?nique behave for a specifie problem. This paper deals with rapid numerical method accurate enough to be used as a design tool. The mathematical model is based on System of Favre averaged Navier-Stokes equations that are written in relative frame of reference, which rotates with constant angular velocity around axis of rotation. The governing equations are solved using finite vol?ume method applied on structured grids. The numerical procedure is based on the explicit multistage Runge-Kutta scheme that is coupled with modem numerical procedures for convergence acceleration. To demonstrate the accuracy of the described numer?ical method developed software is applied to numerical analysis of flow through impeller of axial turbocompressor, and obtained results are compared with available experimental data.


1987 ◽  
Vol 109 (4) ◽  
pp. 345-352 ◽  
Author(s):  
M. Reggio ◽  
R. Camarero

A numerical procedure to solve three-dimensional incompressible flows in arbitrary shapes is presented. The conservative form of the primitive-variable formulation of the time-dependent Navier-Stokes equations written for a general curvilinear coordiante system is adopted. The numerical scheme is based on an overlapping grid combined with opposed differencing for mass and pressure gradients. The pressure and the velocity components are stored at the same location: the center of the computational cell which is used for both mass and the momentum balance. The resulting scheme is stable and no oscillations in the velocity or pressure fields are detected. The method is applied to test cases of ducting and the results are compared with experimental and numerical data.


Author(s):  
Vladimir Shigunov ◽  
Ould el Moctar ◽  
Thomas E. Schellin ◽  
Jan Kaufmann ◽  
Rasmus Stute

The dynamic stability was investigated of a typical offshore service vessel operating under stability critical operating conditions. Excessive roll motions and relative motions at the stern were studied for two loading conditions for ship speeds ranging from zero to the design speed. A linear frequency-domain seakeeping analysis was followed by nonlinear time-domain simulations of ship motions in waves. Based on results from these methods, critical scenarios were selected and simulated using finite-volume solvers of the Reynolds-averaged Navier-Stokes equations to understand the phenomena related to dynamically unstable ship motions as well as to confirm the results of the simpler analysis methods. Results revealed the possibility of excessive roll motions and water run-up on deck; counter measures such as a ship-specific operational guidance are discussed.


1970 ◽  
Vol 12 (6) ◽  
pp. 404-420 ◽  
Author(s):  
S. C. Kacker ◽  
J. H. Whitelaw

An existing numerical procedure for solving the steady, two-dimensional, constant property form of the Navier–Stokes equations, has been used to obtain predictions of mean and fluctuating properties downstream of a two-dimensional wall jet. The Prandtl–Kolmogorov model of turbulence, with a simple empirical expression for the length scale, is shown to permit satisfactory predictions over a wide range of flow situations. These flow situations are relevant to the design of film-cooling slots.


Author(s):  
Pablo Sánchez-Puga ◽  
Javier Tajuelo Rodríguez ◽  
Juan Manuel Pastor ◽  
Miguel Ángel Rubio

Flow field based methods are becoming increasingly popular for the analysis of interfacial shear rheology data. Such methods take properly into account the subphase drag by solving the Navier-Stokes equations for the bulk phases flows, together with the Boussinesq-Scriven boundary condition at the fluid-fluid interface, and the probe equation of motion. Such methods have been successfully implemented at the double wall-ring (DWR), the magnetic rod (MR), and the bicone interfacial shear rheometers. However, a study of the errors introduced directly by the numerical processing is still lacking. Here we report on a study of the errors introduced exclusively by the numerical procedure corresponding to the bicone geometry at an air-water interface. In our study we directly input a preset the value of the complex interfacial viscosity and we numerically obtain the corresponding flow field and the complex amplitude ratio for the probe motion. Then we use the standard iterative procedure to obtain the calculated complex viscosity value. A detailed comparison of the set and calculated complex viscosity values is made upon changing different parameters such as real and imaginary parts of the complex interfacial viscosity and frequency. The observed discrepancies yield a detailed landscape of the numerically introduced errors.


Author(s):  
Kuki Junichi ◽  
Kazuyuki Toda ◽  
Makoto Yamamoto

This paper presents a numerical procedure to predict a three-dimensional sand erosion phenomenon and the interaction between the flow field and the eroded surface. To simulate this phenomenon, the turbulent flow field, the particle trajectory and the amount of erosion on the eroded wall are calculated repeatedly. In computations of the flow field, compressible Navier-Stokes equations and low-Reynolds-number type k–ε turbulence model are adopted. Assuming that the concentration of suspended particle is dilute, particle-particle collision and the influence of particle motions on the flow field are neglected. The Neilson-Gilchrist erosion model is used to estimate the weight loss due to erosion. To verify the developed code, two types of 90-degree bends are computed. The results show that the present procedure can reasonably reproduce the sand erosion process and the temporal change of both the flow field and the wall surface qualitatively.


Sign in / Sign up

Export Citation Format

Share Document