• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.709-714
• /
• 2000

A splitting method for the direct numerical simulation of solid-liquid mixtures is presented, where a symmetric pressure equation is newly proposed. Through numerical experiment, it is found that the newly proposed splitting method works well with a matrix-free formulation fer some bench mark problems avoiding an erroneous pressure field which appears when using the conventional pressure equation of a splitting method. When deriving a typical pressure equation of a splitting method, the motion of a solid particle has to be approximated by the 'intermediate velocity' instead of treating it as unknowns since it is necessary as a boundary condition. Therefore, the motion of a solid particle is treated in such an explicit way that a particle moves by the known form drag (pressure drag) that is calculated from the pressure equation in the previous step. From the numerical experiment, it was shown that this method gives an erroneous pressure field even for the very small time step size as a particle velocity increases. In this paper, coupling the unknowns of particle velocities in the pressure equation is proposed, where the resulting matrix is reduced to the symmetric one by applying the projector of the combined formulation. It has been tested over some bench mark problems and gives reasonable pressure fields.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.1876-1881
• /
• 2003

When a projectile travels at high speed underwater, supercavitating flow arises, in which a huge cavity is generated behind the projectile so that only the nose, i.e., the cavitator, of the projectile is wetted, while the rest of it should be surrounded by the cavity. In that case, the projectile can achieve very high speed due to the reduced drag. Furthermore if the nose of the body is shaped properly, the attendant pressure drag can be maintained at a very low value, so that the overall drag is also reduced dramatically. In this study, shape optimization technique is employed to determine the optimum cavitator shape for minimum drag, given certain operating conditions. Shape optimization technique is also used to solve the potential flow problem for any given cavitator, which is a free boundary value problem having the cavity shape as unknown a priori. Analytical sensitivities are derived for various shape parameters in order to implement a gradient-based optimization algorithm. Simultaneous optimization technique is proposed for efficient cavitator shape optimization, in which the cavity and cavitator shape are determined in a single optimization routine.

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.3
• /
• pp.14-25
• /
• 1998

This paper describes the numerical modelling for the steady and unsteady forces of 3-D wings flying near the free surface based on a potential based panel method. For the steady problem where a wing flies over the fixed float surface, steady lift and drag forces are calculated for wings with and without end-plates having different sections, angle of attacks, aspect patios and flying heights. These numerical results are compared with the wind tunnel test results. The unsteady problem is treated as a boundary value one where a wing flies over the described wavy surface. The unsteady lift force variations of a wing due to different wave lengths and heights are calculated at different flying heights.

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.3
• /
• pp.33-44
• /
• 1992

In the present paper, three types of analytic and numerical method are applied to the analysis of sloshing problem. Analytic solution with linear free-surface boundary condition is introduced and numerical methods are used to analyze flued flow trapped in two-and three-dimensional tanks. Source-distribution method is applied to two- and three-dimensional rectangular tanks and sphere tank. Finite difference method is utilized to compute fluid motion and pressure evolution in two dimensional tanks with girders or slopes. Calculated results are compared with those of experiment or other numerical techniques.

• Journal of Advanced Research in Ocean Engineering
• /
• v.4 no.4
• /
• pp.146-162
• /
• 2018

The unsteady problem of a rigid body impact onto a floating plate is studied. Both the plate and the water are at rest before impact. The plate motion is caused by the impact force transmitted to the plate through an elastic layer with viscous damping on the top of the plate. The hydrodynamic force is calculated by using the second-order model of plate impact by Iafrati and Korobkin (2011). The present study is concerned with the deceleration experienced by a rigid body during its collision with a floating object. The problem is studied also by a fully-nonlinear computational-fluid-dynamics method. The elastic layer is treated with a moving body-fitted grid, the impacting body with an immersed boundary method, and a discrete-element method is used for the contact-force model. The presence of the elastic layer between the impacting bod- ies may lead to multiple bouncing of them, if the bodies are relatively light, before their interaction is settled and they continue to penetrate together into the water. The present study is motivated by ship slamming in icy waters, and by the effect of ice conditions on conventional free-fall lifeboats.

• Coupled systems mechanics
• /
• v.11 no.4
• /
• pp.335-355
• /
• 2022

It's important to note that the number of studies on the lateral vibration of steel liquid storage tanks has been quite modest in the past. The aim of this research has to look at the variables that affect vibration of storage tanks and to highlight the characteristics of a construction that hasn't received much attention in the literature. The storage tank has pre-sized in the study, and aluminum and steel have chosen as components. The specified material qualities and the factors utilized in the investigation has used to calculate vibration frequency values. The resulting calculations are backed up by tables and graphs, and it's an important to look into the parameters that affect the vibration frequencies that will occur on the designed storage tank vary. In the literature, water tanks are usually modelled as lumped masses. The horizontal stiffness of the column on which it is placed is assumed to be constant throughout. This is an approximation method of solving this problem. The column is handled in this study with a more realistic approach that fits the continuum mechanics in the analysis. The reservoir part is incorporated directly into the problem as the boundary condition.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.14 no.1
• /
• pp.8-18
• /
• 2002

In this paper, wave energy absorption of OWC(oscillating water column) device is analyzed. The analytic model consists of a partially immersed circular vertical cylinder open at its end and an air turbine connected with the air chamber. The boundary value problem is decomposed into scattering problem related to scattering by an incident wave in the absence of a pressure variation and radiation problem describing the flow due to an oscillating pressure in the absence of an incident wave. By invoking the continuity of an air flow inside the chamber, the oscillating pressure in a chamber is derived. With oscillating pressure, the mean power absorbed by OWC device and the capture width are obtained. In numerical calculation, the induced volume flux across the internal free surface of the chamber in the scattering and radiation problem and the maximum capture width are compared with various design parameters such as radius and submergence depth of chamber and wave conditions. The maximum capture width obtained by choosing the optimal value of turbine constant occurs at the first resonant mode (Helmholtz mode) among the natural frequencies of a circular cylinder chamber.

• Journal of Astronomy and Space Sciences
• /
• v.15 no.1
• /
• pp.209-220
• /
• 1998

The stability problem of steady motion of a rigid body with multi-elastic appendages and multi-liquid-filled cavities, in the presence of no external forces or torque, is considered in this paper. The flexible appendages are modeled as the clamped -free-free-free rectangular plates, or/and as the discrete mass- spring sub-system. The motion of liquid in every single ellipsoidal cavity is modeled as the uniform vortex motion with a finite number of degrees of freedom. Assuming that stationary holonomic constraints imposed on the body allow its rotation about a spatially fixed axis, the equation of motion for such a systematic configuration can be very complex. It consists of a set of ordinary differential equations for the motion of the rigid body, the uniform rotation of the contained liquids, the motion of discrete elastic parts, and a set of partial differential equations for the elastic appendages supplemented by appropriate initial and boundary conditions. In addition, for such a hybrid system, under suitable assumptions, their equations of motion have four types of first integrals, i.e., energy and area, Helmholtz' constancy of liquid - vortexes, and the constant of the Poisson equation of motion. Chetaev's effective method for constructing Liapunov functions in the form of a set of first integrals of the equations of the perturbed motion is employed to investigate the sufficient stability conditions of steady motions of the complete system in the sense of Liapunov, i.e., with respect to the variables determining the motion of the solid body and to some quantities which define integrally the motion of flexible appendages. These sufficient conditions take into account the vortexes of the contained liquids, the vibration of the flexible components, and coupling among the liquid-elasticity solid.

• Steel and Composite Structures
• /
• v.41 no.6
• /
• pp.899-910
• /
• 2021

A 4-unknown shear deformation theory is applied to investigate the vibration of functionally graded plates under thermal environment. The plate is fabricated from a functionally graded material mixed of ceramic and metal with continuously varying material properties through the plate thickness. Three types of thermal loadings, uniform, linear and nonlinear temperature rises along the plate thickness are taken into account. The present theory contains four unknown functions as against five or more in other higher order shear deformation theories. The through-the-thickness distributions of transverse shear stresses of the plate are considered to vary parabolically and vanish at upper and lower surfaces. The present model does not require any problem dependent shear correction factor. Analytical solutions for the free vibration analysis are derived based on Fourier series that satisfy the boundary conditions (Navier's method). Benchmark solutions are firstly considered to evaluate the accuracy of the proposed model. Comparisons with the solutions available in literature revealed the good capabilities of the present model for the simulations of vibration responses of FG plates. Some parametric studies are carried out for the frequency analysis by varying the volume fraction profile and the temperature distribution across the plate thickness.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.1
• /
• pp.12-18
• /
• 1991

Three-dimensional nonlinear sloshing effects due to tank motions are simulated by solving boundary value problem using the panel method based on boundary integral technique. While Shinkai used boundary elements on which source strengths vary linearly between nodes, the source of constant strength is distributed on each triangular panel in the present study. The source strength at each time step is determined by solving the Fredholm integral equation of the second kind obtained from Green's theorem. To avoid cumulative numerical errors as time elapses, Adam-Bashforth-Moulton method is employed. Numerical examples for the case of partially filled spherical tank on board oscillating in harmonic sway mode or pitch mode are solved. The elevation of the free surface is compared with the result by Shinkai and confirmed in good agreement during early time. The input and the output energy are comparatively evaluated to check the overall accuracy of the present numerical scheme. Although some leakage of energy are found as time marches, it is plausible when we take into account nonlinearities of the problem and the number of panels of the model.