• International Journal of Naval Architecture and Ocean Engineering
• /
• v.10 no.5
• /
• pp.545-565
• /
• 2018

A numerical code is developed based on potential flow theory to investigate nonlinear sloshing in rectangular Liquefied Natural Gas (LNG) tanks under forced excitation. Using this code, internal free-surface elevation and sloshing loads on liquid tanks can be obtained both in time domain and frequency domain. In the mathematical model, acceleration potential is solved in the calculation of pressure on tanks and the artificial damping model is adopted to account for energy dissipation during sloshing. The Boundary Element Method (BEM) is used to solve boundary value problems of both velocity potential and acceleration potential. Numerical calculation results are compared with published results to determine the efficiency and accuracy of the numerical code. Sloshing properties in partially filled rectangular and membrane tank under translational and rotational excitations are investigated. It is found that sloshing under horizontal and rotational excitations share similar properties. The first resonant mode and excitation frequency are the dominant response frequencies. Resonant sloshing will be excited when vertical excitation lies in the instability region. For liquid tank under rotational excitation, sloshing responses including amplitude and phase are sensitive to the location of the center of rotation. Moreover, experimental tests were conducted to analyze viscous effects on sloshing and to validate the feasibility of artificial damping models. The results show that the artificial damping model with modifying wall boundary conditions has better applicability in simulating sloshing under different fill levels and excitations.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2001.09a
• /
• pp.75-82
• /
• 2001

A dynamic interaction analysis of an adjacent surface fecundation on a layered half-space is performed in the frequency domain. A semi-analytical approach is employed to reduce the integration range of the wavenumber in the surface fundamental solution for a layered half-space in boundary element (BE) formulations. The present study then adopts a combined boundary and finite element method to analyze the dynamic behavior of a system of flexible surface foundations on an elastic homogeneous and layered half-space. Numerical examples are presented to demonstrate the accuracy of the developed method. The examples show the feasibility of an extended application fur the complicated dynamic interaction of foundations on layered media.

• Journal of Ocean Engineering and Technology
• /
• v.18 no.1
• /
• pp.1-6
• /
• 2004

In this paper, the performance of a wave absorbing system using an inclined punching plate, was investigated. Using the multi-domain boundary element method (BEM), the reflection coefficients of an inclined punching plate were tested with various design parameters, such as inclined angle, porosity, and wave frequencies. To confirm the numerical solutions, the systematic model test was conducted at 2-D tank and square tank. The numerical results were in good agreement with the experimental results within the entire frequency range. It was found that an inclined punching plate had an excellent wave absorbing efficiency, compared to a horizontal one. Also, the optimal range of an inclined angle had an excellent wave absorbing efficiency, compared to a horizontal one the optimal range of the inclined angle is 10°<β<20°. The developed wave absorber was installed at KRISO's square basin, and is working effectively for various model tests.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.22 no.3
• /
• pp.9-18
• /
• 1985

The numerical calculation for solving boundary-value problem related to potential flows with a free surface is carried out by application of the localized finite element method. Only forced motion of 2-D body in infinitely deep fluid is considered, although this schemes is equally applicable to any first order time-harmonic problems of similar nature. The infinite domain of the fluid is separated into the inner flow field and the outer flow field with common inter-surface boundary. The finite element method is applied to obtain the solution in the inner flow field and the Green functions are utilized to represent the solution in the outer flow field. At the inter-surface boundary, the continuity of the value of potential and the normal derivative of the potential(i.e. matching condition) is conserved. The present method has better computational efficiency than the previous LFEM and the integral equation method of Frank. This enhanced computational efficiency is presumably due to the fact that the present method gives a symmetric coefficient matrix and requires less computational time in calculating the influence coefficient matrix of Green function than the integral equation method. And the irregular frequency desen't exist because the uniqueness of the solution is assured by the such that the exact free surface condition is satisfied on the boundary of the localized finite element region(i.e. inner region). As an example of the above method, the hydrodynamic forces for the circular cylinder and the rectangular cylinders are calculated. In the computed results, the small number of singularity distribution segments($3{\sim}6$) give good result relative to Ursell's and Vugts'.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.19 no.5
• /
• pp.333-339
• /
• 1999

It is well recognized that ultrasonic technique is one of the most common and reliable nondestructive evaluation techniques for quantitative estimation of defects in structures. For the quantitative and accurate estimation of internal defects. the characteristics of scattered ultrasonic wavefields must be understood. In this study. the scattered near-field and far-field due to a circular cavity embedded in infinite media subjected to incident SH-waves were calculated by the boundary element method. The frequency response of the scattered ultrasonic far-field was transformed into the time-domain signal by obtaining its inverse Fourier transform. It was found that the amplitude of time-domain signal decreases and its time delay increases as the distance between the detecting point of ultrasonic scattered field and the center of internal cavity increases.

• Structural Engineering and Mechanics
• /
• v.26 no.4
• /
• pp.363-376
• /
• 2007

A FE-BE procedure is presented for dynamic analysis of concrete arch dams. In this technique, dam body is discretized by finite elements, while foundation rock is handled by three dimensional boundary element formulation. This would allow a rigorous inclusion of dam-foundation rock interaction, with no limitations imposed on geometry of canyon shape. Based on this method, a previously developed program is modified, and the response of Morrow Point arch dam is studied for various ratios of foundation rock to dam concrete elastic moduli under an empty reservoir condition. Furthermore, the effects of canyon shape on response of dam, is also discussed.

• Structural Engineering and Mechanics
• /
• v.33 no.1
• /
• pp.47-66
• /
• 2009

Radiation damping due to wave propagation in unbounded domains may cause a significant reduction of structural vibrations when excited near resonance. Here a novel matrix-valued algebraic Pad$\acute{e}$-like stiffness formulation in the frequency-domain and a corresponding state equation in the time domain are elaborated for a soil-structure interaction problem with a layered soil excited in a transient manner by a flexible rotor during startup and shutdown. The contribution of radiation damping caused by a soil-layer upon a rigid bedrock is characterized by the corresponding amount of critical damping as it is used in structural dynamics.

• Journal of Ocean Engineering and Technology
• /
• v.27 no.1
• /
• pp.80-84
• /
• 2013

This study is to develop the simplified formulae for added mass coefficient of a 2-D floating body with a semi-circle section in a finite water depth. The semi-circle floating body may represent a simplified midship section transformed by Lewis form, which can be used for the ship motion analysis by strip theory. Since the added mass coefficient varies with motion frequencies and sea bottom effect, the correction factor representing the effect of water depth and frequencies is developed for accurate prediction of added mass. Using a two-dimensional numerical wave tank (NWT) technique based on the boundary element method (BEM) including sea bottom boundary the reference values of added mass are calculated to develop the correction factor. For verification and effectiveness of the formulae, the predicted added mass coefficients for various frequencies and water depth ratios are compared with the calculated values from NWT technique.

• Journal of Ocean Engineering and Technology
• /
• v.15 no.2
• /
• pp.22-32
• /
• 2001

The numerical investigation of obliquely incident wave interactions with fully submerged dual buoy/porous-membrane floating breakwaters placed in parallel with spacing is studied based on linear potential theory and Darcy's law. The numerical solutions are obtained by using a discrete-membrane dynamic model and second-kind modified Bessel function distribution over the entire boundaries of fluid regions. First, numerical solutions for an idealized dual submerged system without buoys are obtained. Second, a more practical dual submerged system with membrane tension provided by buoys at its tops is investigated by the multi-domain boundary element method particularly devised for dual buoy/porous-membrane problems with gaps. The velocity potentials of wave motion are coupled with porous-membrane deformation, and solved simultaneously since the boundary condition on porous-membrane is not known in advance. The effects of varying permeability on membranes and wave characteristics are discussed for the optimum design parameters of systems previously studied. The inclusion of permeability on membrane eliminates the resonances that aggravate the breakwater performance. The system is highly efficient when waves generated by the buoys and membranes were mutually canceled and its energy at resonance frequency dissipates through fine pores on membranes.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.18 no.11
• /
• pp.1157-1169
• /
• 2008

This paper presents spectral element formulation which approximates Lamb wave propagation by PZT transducers bonded on a thin plate. A two layer beam model under 2-D plane strain condition is introduced to simulate high-frequency dynamic responses induced by a piezoelectric (PZT) layer rigidly bonded on a base plate. Mindlin-Herrmann and Timoshenko beam theories are employed to represent the first symmetric and anti-symmetric Lamb wave modes on a base plate, respectively. The Euler-Bernoulli beam theory and 1-D linear piezoelectricity are used to model the electro-mechanical behavior of a PZT layer. The equations of motions of a two layer beam model are derived through Hamilton's principle. The necessary boundary conditions associated with the electro-mechanical properties of a PZT layer are formulated in the context of dual functions of a PZT layer as an actuator and a sensor. General spectral shape functions of response field and the associated boundary conditions are obtained through equations of motions converted into frequency domain. Detailed spectrum element formulation for composing the dynamic stiffness matrix of a two layer beam model is presented as well. The validity of the proposed spectral element is demonstrated through numerical examples.