An Investigation on Two-Dimensional Nonlinear Sloshing in Rectangular Tank

Author(s):  
Dongya Zhao ◽  
Zhiqiang Hu ◽  
Gang Chen

Two-dimensional liquid sloshing in rectangular tank of FLNG system is investigated both numerically and experimentally. In numerical simulation, a time-domain scheme has been developed based on potential flow theory in boundary element method. Tank movement is defined by wall boundary condition to produce a reciprocating oscillation. Nonlinear free surface condition is adopted to capture free surface elevation. Energy dissipation caused by viscous effects is considered by applying artificial damping term to the dynamic free surface condition, which is also vital to achieve a steady-state solution. For comparison, experiments of a rectangular tank filled with water subjected to specified oscillation are carried out. As coupling effects between sloshing and tank motion is not included in this research, the testing apparatus is required to produce consistent oscillation movement and not affected by the change of filling condition and sloshing load. Liquid surface elevations in several typical places of the tank were measured. Sloshing related parameters including oscillation amplitude, frequency and filling level are analyzed systematically. It’s found that numerical simulation results have good agreement with phenomenon observed under small amplitude excitation, and this nonlinear analysis method is proved to be effective in capturing liquid surface elevation. It is found that sloshing in tank is sensitive to filling level as well as excitation frequency, especially in the crucial combination cases of them. For given filling level, sloshing tends to be violent near corresponding natural frequencies, and viscous damping has limited contribution to sloshing amplitude when resonance occurs. This fundamental investigation also paves path for the study of more complicated sloshing problems.

Author(s):  
J. Lee ◽  
J.-M. Vanden-Broeck

AbstractThe motion of a two-dimensional bubble rising at a constant velocity U in an inclined tube of width H is considered. The bubble extends downwards without limit, and is bounded on the right by a wall of the tube, and on the left by a free surface. The same flow configuration describes also a jet emerging from a nozzle and falling down along an inclined wall. The acceleration of gravity g and the surface tension T are included in the free surface condition. The problem is characterized by the Froude number the angle β between the left wall and the horizontal, and the angle γ between the free surface and the right wall at the separation point. Numerical solutions are obtained via series truncation for all values of 0 < β < π. The results extend previous calculations of Vanden-Broeck [12–14] for β = π/2 and of Couët and Strumolo [3] for 0 < β < π/2. It is found that the behavior of the solutions depends on whether 0 < β 2π/3 or 2π/3 ≤ β < π. When T = 0, it is shown that there is a critical value F of Froude number for each 0 < β 2π/3 such that solutions with γ = 0, π/3 and π - β occur for F > Fc F = Fc and F < Fc respectively, and that all solutions are characterized by γ = 0 for 2π/3 ≤ β < π. When a small amount of surface tension T is included in the free surface condition, it is found that for each 0 < β < π there exists an infinite discrete set of values of F for which γ = π - β. A particular value F* of the Froude number for which T = 0 and γ = π - β is selected by taking the limit as T approaches zero. The numerical values of F* and the corresponding free surface profiles are found to be in good agreement with experimental data for bubbles rising in an inclined tube when 0 < β < π/2.


Author(s):  
D. C. Hong ◽  
S. Y. Hong ◽  
G. J. Lee ◽  
M. S. Shin

The radiation-diffraction potential of a ship advancing in waves is studied using the three-dimensional frequency-domain forward-speed free-surface Green function (Brard 1948) and the forward-speed Green integral equation (Hong 2000). Numerical solutions are obtained by making use of a second-order inner collocation boundary element method which makes it possible to take account of the line integral along the waterline in a rigorous manner (Hong et al. 2008). The present forward-speed Green integral equation includes not only the usual free surface condition for the potential but also the adjoint free surface condition for the forward-speed free-surface Green function as indicated by Brard (1972). Comparison of the present numerical results of the heave-heave wave damping coefficients and the experimental results for the Wigley ship models I, II and III (Journee 1992) has been presented. These coefficients are compared with those calculated without taking into account of the line integral along the waterline in order to show the forward speed effect represented by the waterline integral when it is properly included in the free-surface Green integral equation. Comparison of the present numerical results and the equivalent time-domain results (Hong et al. 2013) has also been presented.


2009 ◽  
Author(s):  
Jérémie Raymond ◽  
Jean-Marie Finot ◽  
Jean-Michel Kobus ◽  
Gérard Delhommeau ◽  
Patrick Queutey ◽  
...  

The discussion is based on results gathered during the first two years of a 3 years research program for the benefits of Groupe Finot-Conq, Naval Architects. The introduction presents the objectives of the program: Setting up a practical method using numerical and experimental available tools to design fast planing sailing yachts. The aim of this paper is to compare advantages and disadvantages of four different kinds of CFD codes which are linear and non-linear potential flow approach, RANSE solver using finite differences method and RANSE solver using volume of fluid method. The Fluid Mechanics Laboratory of the Ecole Centrale de Nantes (France) has developed those three approaches so those homemade codes will be used for this study. The first one is REVA, a potential flow code with a linearised free surface condition. ICARE is a RANSE solver using finite differences method with a non linear free surface condition. It is extensively used for industrial projects as for sailing yachts projects (ACC for example). ISIS-CFD is a RANSE solver using finite volume method to build the spatial discretization of the transport equations with unstructured mesh. The latter is able to compute sprays for fast planing ships but is also the slower in terms of CPU time. In addition, we had the opportunity to test FS-FLOW which is a potential flow code with a non linear free surface condition distributed by FRIENDSHIP CONSULTING. Numerical results for the four codes are compared with the other codes' results as with tank tests data. Those tank tests were made using captive model test technique on two Open60' models. Reasons of the choice of the captive model technique are explained and experimental procedures are briefly described. Comparisons between codes are mainly based on the easiness of use, the cost in CPU time and the confidence we can have in the results as a function of the boat speed. Flow visualizations, pressure maps, free surface deformation are shown and compared. Analysis of local quantities integrated or by zone is also presented. Results are analyzed focusing on the ability of each code to represent flow dynamics for every speed with a special attention to high speeds. The practical question raised is to know which kind of answers each code can bring in terms of tendencies evaluation or sensitivity to hull geometry modifications. The main goal is to be able to judge if those codes are able to make reliable and consistent comparisons of different designs. Conclusion is that none of the codes is perfect and gather all the advantages. It is still difficult to propose a definitive methodology to estimate hydrodynamic performances at every speed and at every stage of the design process. Knowing each code limitations, it appears more coherent to use each of them at different stages of the design process: the quickest and less reliable to understand the main tendencies and the longest and more precise to validate the final options.


Author(s):  
Takashi Ikeda

Abstract The nonlinear coupled vibrations of an elastic structure and liquid sloshing in a rectangular tank, partially filled with liquid, are investigated. The structure containing the tank is vertically subjected to a sinusoidal excitation. In the theoretical analysis, the resonance curves for the responses of the structure and liquid surface are presented by the harmonic balance method, when the natural frequency of the structure is equal to twice the natural frequency of one of the sloshing modes. From the theoretical analysis, the following predictions have been obtained: (a) Due to the nonlinearity of the fluid force, harmonic oscillations appear in the structure, while subharmonic oscillations occur on the liquid surface, (b) the shapes of the resonance curves markedly change depending on the liquid depth, and (c) when the detuning condition is slightly deviated, almost periodic oscillations and chaotic oscillations appear at certain intervals of the excitation frequency. These were qualitatively in good agreement with the experimental results.


1994 ◽  
Vol 116 (1) ◽  
pp. 91-94 ◽  
Author(s):  
E. P. Rood

An understanding of the process by which vorticity interacts with a free surface is sought by analytical examination of the free-surface condition for the vorticity flux. A novel mechanism is suggested that permits closed vortex loops to evolve into open loops terminating at the free surface. It is hypothesized that abrupt vortex “disconnection,” observed in physical experiments, arises from a smooth diffusion of vorticity through the interface, with a necessary coincident tangential acceleration of the interface attributed to viscous forces.


Author(s):  
D. C. Hong ◽  
Y. Y. Kim ◽  
S. H. Han

The hydrodynamic interaction of two bodies floating in waves is studied. The two-body hydrodynamic coefficients of added mass, wave damping and exciting forces and moments are calculated using the irregular frequency free radiation-diffraction potential solution of the improved Green integral equation associated with the free surface Green function (Hong 1987) according to the conventional two-body analysis. It is well known that the conventional two-body potential solution with usual grid fineness largely overestimates the hydrodynamic coefficients at and near the resonance frequency of the free surface in the gap between two floating bodies moored side-by-side in close proximity (Huijsmans et al. 2001, Hong et al. 2005). The two-body diffraction problem has been solved by both the conventional two-body analysis without damped free surface condition and a boundary matching method with and without damped free surface condition. Numerical results of the wave exciting force coefficients of two identical caissons floating side by side obtained by the two methods have been presented and the discrepancies between them have been discussed. Particular attention is paid to the wave elevation in the gap at the resonance frequencies. Amplitudes and phases of the scattering wave elevations in the gap at the first three free surface resonance frequencies computed by the boundary matching method without damped free surface condition have been presented. It has also been shown that the unrealistic wave elevation due to the resonance of the free surface in the gap can be reduced by imposing the damped free surface condition upon the flow in the gap as used in the oscillating water column hydrodynamics (Hong et al. 2004).


Uniqueness and solvability theorems are proved for a well-posed formulation of the two-dimensional Neumann-Kelvin problem (the modified Neumann-Kelvin problem) in the case, when a body is partly immersed in a supercritical stream. Uniqueness is provided by two supplementary conditions which prescribe (i) additional flux at infinity downstream due to presence of body and (ii) a linear relation between the free-surface elevation at stern point and the velocity circulation along wetted contour. Two versions of source method are developed to find a solution. The first version is simpler, but it fails for some irregular values of the body’s velocity. In the second ver­sion complex sources’ strengths are used, avoiding irregular values and establishing the unique existence of solutions.


1975 ◽  
Vol 69 (2) ◽  
pp. 405-416 ◽  
Author(s):  
G. Dagan

The linearized theory of free-surface gravity flow past submerged or floating bodies is based on a perturbation expansion of the velocity potential in the slenderness parameter e with the Froude number F kept fixed. It is shown that, although the free-wave amplitude and the associated wave resistance tend to zero as F → 0, the linearized solution is not uniform in this limit: the ratio between the second- and first-order terms becomes unbounded as F → 0 with ε fixed. This non-uniformity (called ‘the second Froude number paradox’ in previous work) is related to the nonlinearity of the free-surface condition. Criteria for uniformity of the thin-body expansion, combining ε and F, are derived for two-dimensional flows. These criteria depend on the shape of the leading (and trailing) edge: as the shape becomes finer the linearized solution becomes valid for smaller F.Uniform first-order approximations for two-dimensional flow past submerged bodies are derived with the aid of the method of co-ordinate straining. The straining leads to an apparent displacement of the most singular points of the body contour (the leading and trailing edges for a smooth shape) and, therefore, to an apparent change in the effective Froude number.


2006 ◽  
Vol 3 (6) ◽  
pp. 13189 ◽  
Author(s):  
K Ichinose ◽  
S Yoshida ◽  
K Gomi ◽  
K Taniuchi ◽  
K Fukuda ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document