Advances in the three-dimensional hydroelasticity of ships

Author(s):  
Y-S Wu ◽  
W-C Cui

The study of hydroelasticity of ships first gained momentum in the late 1970s with the work of Bishop and Price, who established the two-dimensional (2D) hydroelasticity theory of ships. The concept and basic principle presented in their work to embody the structure and the surrounding fluid as a coupled entirety was further employed and extended in the creation of the general linear three-dimensional (3D) theory of hydroelasticity for an arbitrary shaped flexible marine structure travelling with a forward speed in a seaway in the middle of 1980s (Wu, 1984; Price and Wu, 1985a; Bishop et al., 1986). Since then, great progress has been achieved in the development and application of 3D hydroelasticity theories. These include the more rigorous methods of frequency-domain linear analysis accounting for the forward speed effect and the steady flow effect, the time-domain linear 3D theory, the non-linear 3D theory and the numerical methods for a floating structure travelling in rough seas with large motions, experimental techniques of 3D flexible ship models, the hydroelasticity-based design and safety assessment, etc. This paper presents an overview of these developments and achievements of linear and non-linear 3D hydroelasticity theories of ships, and the corresponding numerical and experimental techniques.

Author(s):  
L. Salles ◽  
M. Vahdati

The aim of this paper is to study the effects of mistuning on fan flutter and to compare the prediction of two numerical models of different fidelity. The high fidelity model used here is a three-dimensional, whole assembly, time-accurate, viscous, finite-volume compressible flow solver. The Code used for this purpose is AU3D, written in Imperial College and validated for flutter computations over many years. To the best knowledge of authors, this is the first time such computations have been attempted. This is due to the fact that, such non-linear aeroelastic computations with mistuning require large amount of CPU time and cannot be performed routinely and consequently, faster (low fidelity) models are required for this task. Therefore, the second model used here is the aeroelastic fundamental mistuning model (FMM) and it based on an eigenvalue analysis of the linearized modal aeroelastic system with the aerodynamic matrix calculated from the aerodynamic influence coefficients. The influence coefficients required for this algorithm are obtained from the time domain non-linear Code by shaking one blade in the datum (tuned) frequency and mode. Once the influence coefficients have been obtained, the computations of aero damping require minimal amount of CPU time and many different mistuning patterns can be studied. The objectives of this work are to: 1. Compare the results between the two models and establish the capabilities/limitations of aeroelastic FMM, 2. Check if the introduction of mistuning would bring the experimental and computed flutter boundaries closer, 3. Establish a relationship between mistuning and damping. A rig wide-chord fan blade, typical of modern civil designs, was used as the benchmark geometry for this study. All the flutter analyses carried out in this paper are with frequency mistuning, but the possible consequences of mistuned mode shapes are briefly discussed at the end of this paper. Only the first family of modes (1F, first flap) is considered in this work. For the frequency mistuning analysis, the 1F frequency is varied around the annulus but the 1F mode shapes remain the same for all the blades. For the mode shape mistuning computations, an FE analysis of the whole assembly different mass blades is performed. The results of this work clearly show the importance of mistuning on flutter. It also demonstrates that when using rig test data for aeroelastic validation of CFD codes, the amount mistuning present must be known. Finally, it should be noted that the aim of this paper is the study of mistuning and not steady/unsteady validation of a CFD code and therefore minimal aerodynamic data are presented.


Author(s):  
Yoshiyuki Inoue ◽  
Md. Kamruzzaman

The LNG-FPSO concept is receiving much attention in recent years, due to its active usage to exploit oil and gas resources. The FPSO offloads LNG to an LNG carrier that is located close to the FPSO, and during this transfer process two large vessels are in close proximity to each other for daylong periods of time. Due to the presence of neighboring vessel, the motion response of both the vessels will be affected significantly. Hydrodynamic interactions related to wave effects may result in unfavorable responses or the risk of collisions in a multi-body floating system. Not only the motion behavior but also the second order drift forces are influenced by the neighboring structures due to interactions of the waves among the structures. A study is made on the time domain analysis to assess the behavior and the operational capability of the FPSO system moored in the sea having an LNG carrier alongside under environmental conditions such as waves, wind and currents. This paper presents an analysis tool to predict the dynamic motion response and non-linear connecting and mooring forces on a parallel-connected LNG-FPSO system due to non-linear exciting forces of wave, wind and current. Simulation for the mooring performance is also investigated. The three-dimensional source-sink technique has been applied to obtain the radiation forces and the transfer function of wave exciting forces on floating multi-bodies. The hydrodynamic interaction effect between the FPSO and the LNG carrier is included to calculate the hydrodynamic forces. For the simulation of a random sea and also for the generation of time depended wind velocity, a fully probabilistic simulation technique has been applied. Wind and current loads are estimated according to OCIMF. The effects of variations in wave, wind and current loads and direction on the slowly varying oscillations of the LNG and FPSO are also investigated in this paper. Finally, some conclusions are drawn based on the numerical results obtained from the present time domain simulations.


Author(s):  
D. C. Hong ◽  
S. Y. Hong ◽  
H. G. Sung

The radiation and diffraction potentials of a ship advancing in waves are calculated in the time-domain using the three-dimensional time-domain forward-speed free-surface Green function and the Green integral equation on the basis of the Neumann-Kelvin linear wave hypothesis. The Green function approximated by Newman for large time is used together with the Green function by Lamb for small time. The time-domain diffraction problem is solved for the time derivative of the potential by using the time derivative of the impulsive incident wave potential represented by using the complementary complex error function. The integral equation for the potential is discretized according to a second-order boundary element method where the collocation points are located inside the panel. It makes it possible to take account of the line integral along the waterline in a rigorous manner. The six-degree-of-freedom motion and memory functions as well as the diffraction impulse response functions of a hemisphere and the Wigley seakeeping model are presented for various Froude numbers. Comparisons of the wave damping and exciting force and moment coefficients for zero forward speed, calculated by using the Fourier transforms of the time-domain results and the frequency-domain coefficients calculated by using the improved Green integral equation which is free of the irregular frequencies, have been shown to be satisfactory. The wave damping coefficients for non-zero forward speed, calculated by using Fourier transforming of the present time-domain results have also been compared to the experimental results and agreement between them has been shown to be good. A simulation of coupled heave-pitch motion of the Wigley seakeeping model advancing in regular head waves of unit amplitude has been carried out.


2003 ◽  
Vol 9 (7) ◽  
pp. 839-862 ◽  
Author(s):  
Khaled S. Youssef ◽  
Dean T. Mook ◽  
Ali H. Nayfeh ◽  
Saad A. Ragab

Roll motion is an undesirable feature of the behavior of a ship in rough seas, and so it is natural to consider ways of reducing it. The most common devices for increasing roll damping are bilge keels. However, the effectiveness of keels is limited, and anti-roll tanks and fins are used when more control is required. Moreover, unlike keels, anti-roll tanks can be used when the ship is not underway. Our objective is to develop design procedures for passive tanks for roll reduction in rough seas. To this end, we develop an improved model of the passive tank-liquid motion in this paper. This tank consists of U-shaped tubes placed side by side along the length of the ship. The equations of six-degrees-of-motion (6DOF) that govern the tank-liquid are coupled with those that govern the 6DOF motion of the ship, and all of the equations are integrated simultaneously in the time domain using the Large Amplitude Motion Program (“LAMP”). LAMP is a three-dimensional time-domain simulation of the motion of ships in waves. The unstabilized and stabilized roll motions of a S60-70 ship with forward speed and beam waves have been analyzed. For high-amplitude waves, the variation of the roll angle with the encounter-wave frequency exhibits typical nonlinear phenomena: a shift in the resonance frequency, multi-valued responses, and jumps. The performance of passive tanks on a S60-70 ship with forward speed is investigated in an irregular sea with different encounter-wave directions. It is found that passive anti-roll tanks tuned in the nonlinear range are very effective in reducing the roll motion. The effect of the tank mass and distribution of tank tubes on the performance of the tank system is studied. Also, it is found that passive anti-roll tanks are very effective in reducing the roll motion in sea state five of a ship whose pitch frequency is nearly twice its roll frequency.


Author(s):  
G. Jacobs ◽  
F. Theunissen

In order to understand how the algorithms underlying neural computation are implemented within any neural system, it is necessary to understand details of the anatomy, physiology and global organization of the neurons from which the system is constructed. Information is represented in neural systems by patterns of activity that vary in both their spatial extent and in the time domain. One of the great challenges to microscopists is to devise methods for imaging these patterns of activity and to correlate them with the underlying neuroanatomy and physiology. We have addressed this problem by using a combination of three dimensional reconstruction techniques, quantitative analysis and computer visualization techniques to build a probabilistic atlas of a neural map in an insect sensory system. The principal goal of this study was to derive a quantitative representation of the map, based on a uniform sample of afferents that was of sufficient size to allow statistically meaningful analyses of the relationships between structure and function.


Algorithms ◽  
2021 ◽  
Vol 14 (5) ◽  
pp. 129
Author(s):  
Yuan Li ◽  
Ni Zhang ◽  
Yuejiao Gong ◽  
Wentao Mao ◽  
Shiguang Zhang

Compared with continuous elements, discontinuous elements advance in processing the discontinuity of physical variables at corner points and discretized models with complex boundaries. However, the computational accuracy of discontinuous elements is sensitive to the positions of element nodes. To reduce the side effect of the node position on the results, this paper proposes employing partially discontinuous elements to compute the time-domain boundary integral equation of 3D elastodynamics. Using the partially discontinuous element, the nodes located at the corner points will be shrunk into the element, whereas the nodes at the non-corner points remain unchanged. As such, a discrete model that is continuous on surfaces and discontinuous between adjacent surfaces can be generated. First, we present a numerical integration scheme of the partially discontinuous element. For the singular integral, an improved element subdivision method is proposed to reduce the side effect of the time step on the integral accuracy. Then, the effectiveness of the proposed method is verified by two numerical examples. Meanwhile, we study the influence of the positions of the nodes on the stability and accuracy of the computation results by cases. Finally, the recommended value range of the inward shrink ratio of the element nodes is provided.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Souma Jinno ◽  
Shuji Kitora ◽  
Hiroshi Toki ◽  
Masayuki Abe

AbstractWe formulate a numerical method on the transmission and radiation theory of three-dimensional conductors starting from the Maxwell equations in the time domain. We include the delay effect in the integral equations for the scalar and vector potentials rigorously, which is vital to obtain numerically stable solutions for transmission and radiation phenomena in conductors. We provide a formalism to connect the conductors to any passive lumped-parameter circuits. We show one example of numerical calculations, demonstrating that the new formalism provides stable solutions to the transmission and radiation phenomena.


Sign in / Sign up

Export Citation Format

Share Document