• International Journal of Naval Architecture and Ocean Engineering
• /
• v.8 no.2
• /
• pp.153-168
• /
• 2016

Sloshing phenomenon is a complicated free surface flow problem that increases the dynamic pressure on the sidewalls and the bottom of the storage tanks. When the storage tanks are partially filled, it is essential to be able to evaluate the fluid dynamic loads on the tank's perimeter. In this paper, a numerical code was developed to determine the pressure distribution on the rectangular and trapezoidal storage tanks' perimeters due to liquid sloshing phenomenon. Assuming the fluid to be inviscid, the Laplace equation and the nonlinear free surface boundary conditions were solved using coupled boundary element - finite element method. The code performance for sloshing modeling was validated using Nakayama and Washizu's results. Finally, this code was used for partially filled rectangular and trapezoidal storage tanks and free surface displacement, pressure distribution and horizontal and vertical forces exerted on the tanks' perimeters due to liquid sloshing phenomenon were estimated and discussed.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.5 no.3
• /
• pp.204-211
• /
• 1993

Boundary element method is applied to simulate nonlinear water waves using Green's identity formula in a numerical wave flume. A system of linear equations is formulated from the governing equation and free surface boundary conditions in order to calculate velocity potential and water surface elevation at each nodal point. The velocity square terms are included in the dynamic free surface boundary condition. The free surface is treated as a moving boundary. the vertical variation of velocity potential being considered in calculating the time derivative of the velocity potential at the free surface. The present method is applied to simulate solitary wave and Stokes 2nd order wave, and shows excellent agreements with their theoretical values.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.119-125
• /
• 2004

Unsteady nonlinear wave motions on the free surface over a plane beach of constant slope are numerically simulated using a finite difference method in rectangular grid system. Two-dimensional Navier-Stokes equations and the continuity equation are used for the computations. Irregular leg lengths and stars are employed near the boundaries of body and free surface to satisfy the boundary conditions. Also, the free surface which consists of markers or segments is determined every time step with the satisfaction of kinematic and dynamic free surface conditions. Moreover, marker-density method is also adopted to allow plunging jets impinging on the free surface. The second-order Stokes wave theory and solitary wave theory are employed for the generation of waves on the inflow boundary. For the simulation of wave breaking phenomena, the computations are carried out with the plane beach of constant slope in surf zone. The results are compared with each other. The marker-density method is better then the hybrid method. Also they are compared with other existing experimental results. The Agreement between the experimental data and the computation results is good.

• Journal of Ship and Ocean Technology
• /
• v.8 no.2
• /
• pp.61-69
• /
• 2004

Unsteady nonlinear wave motions on the free surface over a plane beach of constant slope are numerically simulated using a finite difference method in rectangular grid system. Two-dimensional Navier-Stokes equations and the continuity equation are used for the computations. Irregular leg lengths and stars are employed near the boundaries of body and free surface to satisfy the boundary conditions. Also, the free surface which consists of markers or segments is determined every time step with the satisfaction of kinematic and dynamic free surface conditions. Moreover, marker-density method is also adopted to allow plunging jets impinging on the free surface. The second-order Stokes wave theory is employed for the generation of waves on the inflow boundary. For the simulation of wave breaking phenomena, the computations are carried out with the plane beach of constant slope in surf zone. The results are compared with other existing experimental results. Agreement between the experimental data and the computation results is good.

• Steel and Composite Structures
• /
• v.35 no.3
• /
• pp.449-462
• /
• 2020

This work aims to study effects of the crack and the surface energy on the free longitudinal vibration of axially functionally graded nanorods. The surface energy parameters considered are the surface stress, the surface density, and the surface Lamé constants. The cracked nanorod is modelled by dividing it into two parts connected by a linear spring in which its stiffness is related to the crack severity. The surface and bulk material properties are considered to vary in the length direction according to the power law distribution. Hamilton's principle is implemented to derive the governing equation of motion and boundary conditions. Considering the surface stress causes that the derived governing equation of motion becomes non-homogeneous while this was not the case in works that only the surface density and the surface Lamé constants were considered. To extract the frequencies of nanorod, firstly the non-homogeneous governing equation is converted to a homogeneous one using an appropriate change of variable, and then for clamped-clamped and clamped-free boundary conditions the governing equation is solved using the harmonic differential quadrature method. Since the present work considers effects of all the surface energy parameters, it can be claimed that this is a comprehensive work in this regard.

• Journal of KSNVE
• /
• v.4 no.2
• /
• pp.137-146
• /
• 1994

An analytical method for linear free vibration of fully liquid-filled circular cylindrical shell with various boundary conditions is developed by the Fourier series expansion based on the Stokes' transformation. A set of modal displacement functions and their derivatives of a circular cylindrical shell is substituted into the Sanders' shell equations in order to explicitily represent the Fourier coefficients as functions of the end point displacements, forces, and moments. For the vibration relevant to the liquid motion, the velocity potential of liquid is assumed as a sum of linear combination of suitable harmonic functions in the axial directions. The unknown parameter of the velocity potential is selected to satisfy the boundary condition along the wetted shell surface. An explicit expression of the natural frequency equation can be obtained for any kind of classical boundary conditions. The natural frequencies of the liquid-filled cylindrical shells with the clamped-free, the clamped-clamped, and the simply supported-simply supported boundary conditions examined in the previous works, are obtained by the analytical method. The results are compared with the previous works, and excellent agreement is found for the natural frequencies of the shells.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.156-161
• /
• 1998

This describes a numerical method for predicting the incompressible unsteady laminar three-dimensional flows of fluid behaviour with free-surface. The elliptic differential equations governing the flows have been linearized by means of finite-difference approximations, and the resulting equations have been solved via a fully-implicit iterative method. The free-surface is defined by the motion of a set of marker particles and interface behaviour was investigated by way of a 'Lagrangian' technique. Using the GALA concept of Spalding, the conventional mass continuity equation is modified to form a volumetric or bulk-continuity equation. The use of this bulk-continuity relation allows the hydrodynamic variables to be computed over the entire flow domain including both liquid and gas regions. Thus, the free-surface boundary conditions are imposed implicitly and the problem formulation is greatly simplified. The numerical procedure is validated by comparing the predicted results of a periodic standing waves problems with analytic solutions or experimental results from the literature. The results show that this numerical method produces accurate and physically realistic predictions of three-dimensional free-surface flows.

• International Journal of Ocean Engineering and Technology Speciallssue:Selected Papers
• /
• v.3 no.1
• /
• pp.14-22
• /
• 2000

The improved Green integral equation for the calculation of time-harmonic potentials in the radiation diffraction problem about a freely floating body in the presence of moderate or weak current is presented. The forward-speed Green function presented by Brard is used. The correct free surface boundary conditions on the physical free surface are employed as well as an appropriate boundary conditions on the non-physical inner free surface. The default in the existing Green integral equation as well as in the source integral equation is discussed in detail.

• Ocean Systems Engineering
• /
• v.6 no.3
• /
• pp.245-264
• /
• 2016

The iterative boundary element method (IBEM) developed originally before for cavitating two-dimensional (2-D) and three-dimensional (3-D) hydrofoils moving under free surface is modified and applied to the case of 2-D (two-dimensional) airfoils and 3-D (three-dimensional) wings over water. The calculation of the steady-state flow characteristics of an inviscid, incompressible fluid past 2-D airfoils and 3-D wings above free water surface is of practical importance for air-assisted marine vehicles such as some racing boats including catamarans with hydrofoils and WIG (Wing-In-Ground) effect crafts. In the present paper, the effects of free surface both on 2-D airfoils and 3-D wings moving steadily over free water surface are investigated in detail. The iterative numerical method (IBEM) based on the Green's theorem allows separating the airfoil or wing problems and the free surface problem. Both the 2-D airfoil surface (or 3-D wing surface) and the free surface are modeled with constant strength dipole and constant strength source panels. While the kinematic boundary condition is applied on the airfoil surface or on the wing surface, the linearized kinematic-dynamic combined condition is applied on the free surface. The source strengths on the free surface are expressed in terms of perturbation potential by applying the linearized free surface conditions. No radiation condition is enforced for downstream boundary in 2-D airfoil and 3-D wing cases and transverse boundaries in only 3-D wing case. The method is first applied to 2-D NACA0004 airfoil with angle of attack of four degrees to validate the method. The effects of height of 2-D airfoil from free surface and Froude number on lift and drag coefficients are investigated. The method is also applied to NACA0015 airfoil for another validation with experiments in case of ground effect. The lift coefficient with different clearance values are compared with those of experiments. The numerical method is then applied to NACA0012 airfoil with the angle of attack of five degrees and the effects of Froude number and clearance on the lift and drag coefficients are discussed. The method is lastly applied to a rectangular 3-D wing and the effects of Froude number on wing performance have been investigated. The numerical results for wing moving under free surface have also been compared with those of the same wing moving above free surface. It has been found that the free surface can affect the wing performance significantly.

• Structural Engineering and Mechanics
• /
• v.40 no.5
• /
• pp.665-688
• /
• 2011

This paper presents a theoretical analysis for the free vibration of rectangular tanks partially filled with an ideal liquid. Wet dynamic displacements of the tanks are approximated by combining the orthogonal polynomials satisfying the boundary conditions, since the rectangular tanks are composed of four rectangular plates. The classical boundary conditions of the tanks at the top and bottom ends are considered, such as clamped, simply supported, and clamped-free boundary conditions. As the facing rectangular plates are assumed to be geometrically and structurally identical, the vibration modes of the facing plates of the tanks can be divided into two categories: symmetric and antisymmetric modes with respect to the planes passing through the center of the tanks and perpendicular to the free liquid surface. The liquid displacement potentials satisfying the Laplace equation and liquid boundary conditions are derived, and the wet dynamic modal functions of a quarter of the tanks can be expanded by the finite Fourier transform for compatibility requirements along the contacting surfaces between the tanks and liquid. An eigenvalue problem is derived using the Rayleigh-Ritz method. Consequently, the wet natural frequencies of the rectangular tanks can be extracted. The proposed analytical method is verified by observing an excellent agreement with three-dimensional finite element analysis results. The effects of the liquid level and boundary condition at the top and bottom edges are investigated.