Three-Dimensional Large Amplitude Body Motions in Waves

Author(s):  
Xinshu Zhang ◽  
Robert F. Beck

Three-dimensional, time-domain, wave-body interactions are studied in this paper for cases with and without forward speed. In the present approach, an exact body boundary condition and linearized free surface boundary conditions are used. By distributing desingularized sources above the calm water surface and using constant-strength panels on the exact submerged body surface, the boundary integral equations are solved numerically at each time step. Once the fluid velocities on the free surface are computed, the free surface elevation and potential are updated by integrating the free surface boundary conditions. After each time step, the body surface and free surface are regrided due to the instantaneous changing submerged body geometry. The desingularized method applied on the free surface produces non-singular kernels in the integral equations by moving the fundamental singularities a small distance outside of the fluid domain. Constant strength panels are used for bodies with any arbitrary shape. Extensive results are presented to validate the efficiency of the present method. These results include the added mass and damping computations for a hemisphere. The calm water wave resistance for a submerged spheroid and a Wigley hull are also presented. All the computations with forward speed are started from rest and proceed until a steady state is reached. Finally, the time-domain forced motion results for a modified Wigley hull with forward speed are shown and compared with the experiments for both linear computations and body-exact computations.

Author(s):  
Xinshu Zhang ◽  
Robert F. Beck

Three-dimensional, time-domain, wave-body interactions are studied in this paper for cases with and without forward speed. In the present approach, an exact body boundary condition and linearized free surface boundary conditions are used. By distributing desingularized sources above the calm water surface and using constant-strength flat panels on the exact wetted body surface, the boundary integral equations are numerically solved at each time step. Once the fluid velocities on the free surface are computed, the free surface elevation and potential are updated by integrating the free surface boundary conditions. After each time step, the body surface and free surface are regrided due to the instantaneous wetted body geometry. The desingularized method applied on the free surface produces nonsingular kernels in the integral equations by moving the fundamental singularities a small distance outside of the fluid domain. Constant-strength flat panels are used for bodies with any arbitrary shape. Extensive results are presented to validate the efficiency of the present method. These results include the added mass and damping computations for a hemisphere. The calm water wave resistance for a submerged spheroid and a Wigley hull are also presented. All the computations with forward speed are started from rest and proceeded until a steady state is reached. Finally, the time-domain forced motion results for a modified Wigley hull with forward speed are shown and compared to the experiments for both linear computations and body-exact computations.


2010 ◽  
Vol 54 (02) ◽  
pp. 79-94 ◽  
Author(s):  
Xinshu Zhang ◽  
Piotr Bandyk ◽  
Robert F. Beck

Large-amplitude, time-domain, wave-body interactions are studied in this paper for problems with forward speed. Both two-dimensional strip theory and three-dimensional computation methods are shown and compared by a number of numerical simulations. In the present approach, an exact body boundary condition and linearized free surface boundary conditions are used. By distributing desingularized sources above the calm water surface and using constant-strength flat panels on the exact body surface, the boundary integral equations are solved numerically at each time step. The strip theory method implements Radial Basis Functions to approximate the longitudinal derivatives of the velocity potential on the body. Once the fluid velocities on the free surface are computed, the free surface elevation and potential are updated by integrating the free surface boundary conditions. After each time step, the body surface and free surface are regrided due to the instantaneous changing wetted body geometry. Extensive results are presented to validate the efficiency of the present methods. These results include the added mass and damping computations for a Wigley III hull and an S-175 hull with forward speed using both two-dimensional and three-dimensional approaches. Exciting forces acting on a Wigley III hull due to regular head seas are obtained and compared using both the fully three-dimensional method and the two-dimensional strip theory. All the computational results are compared with experiments or other numerical solutions.


1999 ◽  
Vol 396 ◽  
pp. 37-71 ◽  
Author(s):  
LEONID BREVDO ◽  
PATRICE LAURE ◽  
FREDERIC DIAS ◽  
THOMAS J. BRIDGES

The film flow down an inclined plane has several features that make it an interesting prototype for studying transition in a shear flow: the basic parallel state is an exact explicit solution of the Navier–Stokes equations; the experimentally observed transition of this flow shows many properties in common with boundary-layer transition; and it has a free surface, leading to more than one class of modes. In this paper, unstable wavepackets – associated with the full Navier–Stokes equations with viscous free-surface boundary conditions – are analysed by using the formalism of absolute and convective instabilities based on the exact Briggs collision criterion for multiple k-roots of D(k, ω) = 0; where k is a wavenumber, ω is a frequency and D(k, ω) is the dispersion relation function.The main results of this paper are threefold. First, we work with the full Navier–Stokes equations with viscous free-surface boundary conditions, rather than a model partial differential equation, and, guided by experiments, explore a large region of the parameter space to see if absolute instability – as predicted by some model equations – is possible. Secondly, our numerical results find only convective instability, in complete agreement with experiments. Thirdly, we find a curious saddle-point bifurcation which affects dramatically the interpretation of the convective instability. This is the first finding of this type of bifurcation in a fluids problem and it may have implications for the analysis of wavepackets in other flows, in particular for three-dimensional instabilities. The numerical results of the wavepacket analysis compare well with the available experimental data, confirming the importance of convective instability for this problem.The numerical results on the position of a dominant saddle point obtained by using the exact collision criterion are also compared to the results based on a steepest-descent method coupled with a continuation procedure for tracking convective instability that until now was considered as reliable. While for two-dimensional instabilities a numerical implementation of the collision criterion is readily available, the only existing numerical procedure for studying three-dimensional wavepackets is based on the tracking technique. For the present flow, the comparison shows a failure of the tracking treatment to recover a subinterval of the interval of unstable ray velocities V whose length constitutes 29% of the length of the entire unstable interval of V. The failure occurs due to a bifurcation of the saddle point, where V is a bifurcation parameter. We argue that this bifurcation of unstable ray velocities should be observable in experiments because of the abrupt increase by a factor of about 5.3 of the wavelength across the wavepacket associated with the appearance of the bifurcating branch. Further implications for experiments including the effect on spatial amplification rate are also discussed.


Author(s):  
D. C. Hong ◽  
T. B. Ha ◽  
K. H. Song

The added resistance of a ship was calculated using Maruo’s formula [1] involving the three-dimensional Kochin function obtained using the source and normal doublet distribution over the wetted surface of the ship. The density of the doublet distribution was obtained as the solution of the three-dimensional frequency-domain forward-speed Green integral equation containing the exact line integral along the waterline. Numerical results of the Wigley ship models II and III in head seas, obtained by making use of the inner-collocation 9-node second-order boundary element method have been compared with the experimental results reported by Journée [2]. The forward-speed hydrodynamic coefficients of the Wigley models have shown no irregular-frequencylike behavior. The steady disturbance potential due to the constant forward speed of the ship has also been calculated using the Green integral equation associated with the steady forward-speed free-surface Green function since the so-called mj-terms [3] appearing in the body boundary conditions contain the first and second derivatives of the steady potential over the wetted surface of the ship. However, the free-surface boundary condition was kept linear in the present study. The added resistances of the Wigley II and III models in head seas obtained using Maruo’s formula showing acceptable comparison with experimental results, have been presented. The added resistances in following seas obtained using Maruo’s formula have also been presented.


Author(s):  
Dakui Feng ◽  
Xianzhou Wang ◽  
Zhiguo Zhang ◽  
Yanming Guan

The catamaran is composed of two monohulls, the flow fields between the inner and outer side of each monohull are different, the bodies must be considered as lifting bodies. So it is very important to know the lifting effect on hydrodynamic characteristics of catamaran hull at the preliminary design stage of its hull form. The pressure Kutta condition is imposed on the trailing-surface of the lifting body by determining the dipole distribution, which generates required circulation on the lifting part. The method is based on Green’s second theorem. Rankine Sources and dipoles are placed on boundary surfaces. Time-stepping scheme is adopted to simulate the wave generated by the catamaran with a uniform speed in deep water. The values of the potential and position of the free surface are updated by integrating the nonlinear Lagrangian free surface boundary conditions for every time. A moving computational window is used in the computations by truncating the fluid domain (the free surface) into a computational domain. The grid regeneration scheme is developed to determine the approximate position of the free surface for the next time step. An implicit implement of far field condition is enforced automatically at the truncation boundary of the computational window, Radiation condition is satisfied automatically. The influences on the wave making resistance of the distance between the twin hulls of the Wigley catamaran on the hydrodynamic characteristics are discussed. The numerical results are presented compared with the existing simulation result. The method can be used to simulate the flow fields around the foil near free surface.


Author(s):  
Xinshu Zhang ◽  
Wei Li ◽  
Yunxiang You

A three-dimensional time-domain approach has been developed to compute large-amplitude motion response and the second-order added wave resistance for ships traveling in waves. The proposed method is an extension of a well established linear approach developed in a previous paper [1]. The numerical model is developed based on boundary integral equation, which is solved at each time step by distributing desingularized sources above the calm water surface and employing constant-strength panels on body surface. The nonlinear Froude-Krylov and wave diffraction forces are computed. Equations of motion are solved with including the effects of Euler angles. A broad range of different hull forms, including two Wigley hulls, a Series 60 hull, and a S-175 hull, are employed to validate the present computational model. By comparing the obtained numerical results to experiments, it is demonstrated that the present model using double-body basis flow can well predict added wave resistance.


2016 ◽  
Author(s):  
Limin Chen ◽  
Guanghua He ◽  
Dazheng Wang ◽  
Zihao Zhang

A time-domain seakeeping numerical model based on a computational fluid dynamics (CFD) software FINE/Marine has been developed for nonlinear steady and unsteady viscous flows. Simulation of multi-phase flows around a Wigley hull with forward speed is performed by solving the Reynolds-average Navier-Stokes (RANS) and continuity equations with k-ω (SST-Menter) turbulence model. The water free surface is captured by Blend Reconstruction Interface Capturing Scheme (BRICS). Both steady and unsteady problems including wave-making, radiation and diffraction problems are simulated. Ship waves generated by the Wigley model advancing at a constant forward speed in calm water or incident waves are computed. The numerical results including the wave-making resistance and wave patterns for steady problem, hydrodynamic coefficients and forces for unsteady problems are illustrated and compared with experimental measurements in good agreement. It is confirmed that the present numerical model has the capability of evaluating the seakeeping performance of ships.


2013 ◽  
Vol 57 (01) ◽  
pp. 13-23
Author(s):  
Wei Qiu ◽  
Hongxuan (Heather) Peng

Motions of a floating body in waves are computed in the time domain by solving the body-exact problem with the panel-free method and exact geometry. In the present study, the body boundary condition is imposed on the instantaneous wetted surface exactly at each time step. The free surface boundary is assumed linear so that the time-domain Green function can be applied. The body geometry is represented by NonUniform Rational B-Spline surfaces. At each time step, the instantaneous wetted surface is obtained by trimming the entire body surface. With the panel-free method, the body-exact problems are solved without involving repanelization of the wetted hull surface at each time step. Validation studies have been carried out for a submerged sphere, a flared body, and a Wigley hull. The hydrodynamic forces on the submerged sphere undergoing large-amplitude motion were computed and compared with analytical solutions. For the flared body oscillating in a free surface and the Wigley hull in waves, numerical results were compared with experimental data and solutions by other numerical methods.


1991 ◽  
Vol 35 (04) ◽  
pp. 288-294
Author(s):  
F. Noblesse ◽  
D. M. Hendrix ◽  
L. Kahn

A nonlinear local analysis of the steady potential flow at a ship bow and stern, and more generally at any point along a ship waterline, is presented. The hull boundary condition and the nonlinear kinematic and dynamic free-surface boundary conditions are satisfied exactly, at the actual position of the free surface, in this analysis. The bow-flow analysis shows that the free surface at a ship bow is tangent to the stem. This theoretical result appears to agree with existing experimental measurements of steady bow waves of the Wigley hull. Simple analytical expressions defining the fluid velocity at the bow and the stern, and more generally at any point along the wave profile, in terms of the elevation of the free surface at the corresponding point are also given. These analytical expressions and the available experimental measurements of wave profiles along the Wigley hull show that the velocity of the flow disturbance due to this hull is fairly small compared to the hull speed everywhere along the wave profile except in very small regions around the bow and the stern, where the total fluid velocity is nearly equal to the hull speed in magnitude but directed vertically. Nonlinearities therefore appear to be quite important, although only in very small regions surrounding a ship bow and stern. A genuine nonlinear method of calculation must then be able to represent the very rapid variation in the direction of the fluid velocity occurring within small regions around a ship bow and stern. In particular, a sufficiently fine discretization is required in these regions.


Sign in / Sign up

Export Citation Format

Share Document