• Journal of the korean society of oceanography
• /
• v.31 no.4
• /
• pp.173-182
• /
• 1996

A median-density fluid was injected into the upper layer of a two-layered fluid in a rotating cylindrical container. Several sets of the top and bottom boundary configurations were employed and the flow pattern of each layer including the injected fluid was observed to determine the factors that affect the path of the injected intermediate fluid. The axisymmetric path of the intermediate fluid when the upper layer had a free surface, changed into the asymmetric path with bulged-shape radial spreading whenever either the upper layer or the lower layer had ${\beta}$-effect. The internal Fronds number that controls the shape of the interface turned out to be the most important parameter that determines the radial spreading in terms of location and strength. When the upper and lower layer had the ${\beta}$-effect, convective overturning produced anticyclonic vortices at the frontal edge of the intermediate fluid, and that could enhance the vertical mixing of different density fluids. The intermediate fluid did not produce any topographic effect on the upper-layer motion during its spreading over the interface, since its thickness was very small. However, its anticyclonic motion within the bulged-shape produced a cyclonic motion in the lower layer just beneath the bulge.

• Journal of Mechanical Science and Technology
• /
• v.19 no.spc1
• /
• pp.292-304
• /
• 2005

In this paper we present the linearized equations of motion for a bicycle as a benchmark. The results obtained by pencil-and-paper and two programs are compared. The bicycle model we consider here consists of four rigid bodies, viz. a rear frame, a front frame being the front fork and handlebar assembly, a rear wheel and a front wheel, which are connected by revolute joints. The contact between the knife-edge wheels and the flat level surface is modelled by holonomic constraints in the normal direction and by non-holonomic constraints in the longitudinal and lateral direction. The rider is rigidly attached to the rear frame with hands free from the handlebar. This system has three degrees of freedom, the roll, the steer, and the forward speed. For the benchmark we consider the linearized equations for small perturbations of the upright steady forward motion. The entries of the matrices of these equations form the basis for comparison. Three diffrent kinds of methods to obtain the results are compared : pencil-and-paper, the numeric multibody dynamics program SPACAR, and the symbolic software system Auto Sim. Because the results of the three methods are the same within the machine round-off error, we assume that the results are correct and can be used as a bicycle dynamics benchmark.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.2
• /
• pp.86-97
• /
• 1993

A numerical method is developed for the nonlinear motion of two-dimensional wedges and axisymmetric-forced-heaving motion using Semi-Largrangian scheme under assumption of potential flows. In two-dimensional-problem Cauchy's integral theorem is applied to calculate the complex potential and its time derivative along boundary. In three-dimensional-problem Rankine ring sources are used in a Green's theorem boundary integral formulation to salve the field equation. The solution is stepped forward numerically in time by integrating the exact kinematic and dynamic free-surface boundary condition. Numerical computations are made for the entry of a wedge with a constant velocity and for the forced harmonic heaving motion from rest. The problem of the entry of wedge compared with the calculated results of Champan[4] and Kim[11]. By Fourier transform of forces in time domain, added mass coefficient, damping coefficient, second harmonic forces are obtained and compared with Yamashita's experiment[5].

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.4 no.4
• /
• pp.187-200
• /
• 1992

A general scheme is developed to determine the Langrangian motions of water particles by the Eulerian velocity at their mean positions by using Taylor's theorem. Utilizing the Stokes finite-amplitude wave theory, the orbital motions and the mass transport velocity including the effects of higher-order wave components are determined. The fifth-order approximation of orbital motion gives very good predictions of actual water particle motion in Stokes fifth-order wave theory except near the free-surface. The fifth-order theory predicts the mass transport velocity less than that given by the existing second-order theory over the whole water depth.

• Journal of the Society of Naval Architects of Korea
• /
• v.45 no.2
• /
• pp.109-123
• /
• 2008

The present paper presents the CFD result for a beam wave test case. An ONR tumblehome ship model with bilge keels is used. The beam wave test is for zero forward speed and roll and heave 2DOF with wave slope $a_k=0.156$ and wavelength ${\lambda}=1.12L_{PP}$, with $L_{PP}$ the ship length. The problems is solved numerically with an unsteady Reynolds averaged Navier-Stokes approach. The free surface flow is computed using a single-phase level-set method and the motions in each time step are integrated using a predictor-corrector iteration approach which uses dynamic overset grids moving with relative ship motion. The predicted CFD results for motions and forces are compared with experimental data, showing a reasonable agreement.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.12 no.1
• /
• pp.667-679
• /
• 2020

As the LNG carrier operates in ice covered waters, it is key to ensure the overall safety, which is related to the coupling effect of ice-breaking process and internal liquid sloshing. This paper focuses on the sloshing simulation of the ice-breaking LNG carrier, and the numerical method is proposed using Circumferential Crack Method (CCM) and Volume of Vluid (VOF) with two main key factors (velocity νx and force Fx). The ship motion analysis is carried out by CCM when the ship navigates in the ice-covered waters with a constant propulsion power. The velocity νx is gained, which is the initial excitation condition for the calculation of internal sloshing force Fx. Then, the ship motion is modified based on iterative computations under the union action of ice-breaking force and liquid sloshing load. The sloshing simulation under the LNG tank is studied with the modified ship motion. Moreover, an ice-breaking LNG ship with three-leaf type tank is used for case study. The internal LNG sloshing is simulated with three different liquid heights, including free surface shape and sloshing pressure distribution at a given moment, pressure curves at monitoring points on the bulkhead. This present method is effective to solve the sloshing simulation during ice-breaking process, which could be a good reference for the design of the polar ice-breaking LNG carrier.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.12 no.3
• /
• pp.257-264
• /
• 2010

This study proposed the aseismic design concept for underground space based on site response analysis and laboratory tests. The results of this study showed that the location of the control points of input motions such as design response spectra and time history of acceleration and the assumption of bedrock properties such as elasticity or rigidity play an important role in aseismic design of underground space. Therefore, the appropriate ground response model among models applying motions such as free surface motion, bedrock motion, or bedrock outcropping motion must be utilized to provide reasonable boundary conditions of underground space under earthquake loading and practical aseismic design.

• Journal of the Korean Mathematical Society
• /
• v.49 no.3
• /
• pp.585-604
• /
• 2012

In this paper we consider the evolution of the rolling stone with a rotationally symmetric nonconvex compact initial surface ${\Sigma}_0$ under the Gauss curvature flow. Let $X:S^n{\times}[0,\;{\infty}){\rightarrow}\mathbb{R}^{n+1}$ be the embeddings of the sphere in $\mathbb{R}^{n+1}$ such that $\Sigma(t)=X(S^n,t)$ is the surface at time t and ${\Sigma}(0)={\Sigma}_0$. As a consequence the parabolic equation describing the motion of the hypersurface becomes degenerate on the interface separating the nonconvex part from the strictly convex side, since one of the curvature will be zero on the interface. By expressing the strictly convex part of the surface near the interface as a graph of a function $z=f(r,t)$ and the non-convex part of the surface near the interface as a graph of a function $z={\varphi}(r)$, we show that if at time $t=0$, $g=\frac{1}{n}f^{n-1}_{r}$ vanishes linearly at the interface, the $g(r,t)$ will become smooth up to the interface for long time before focusing.

• Journal of Ocean Engineering and Technology
• /
• v.24 no.5
• /
• pp.16-21
• /
• 2010

During the loading/offloading operation of a liquefied natural gas carrier (LNGC) that is moored in a side-by-side configuration with an offshore plant, sloshing that occurs due to the partially filled LNG tank and the interactive effect between the two floating bodies are important factors that affect safety and operability. Therefore, a time-domain software program, called CHARM3D, was developed to consider the interactions between sloshing and the motion of a floating body, as well as the interactions between multiple bodies using the potential-viscous hybrid method. For the simulation of a floating body in the time domain, hydrodynamic coefficients and wave forces were calculated in the frequency domain using the 3D radiation/diffraction panel program based on potential theory. The calculated values were used for the simulation of a floating body in the time domain by convolution integrals. The liquid sloshing in the inner tanks is solved by the 3D-FDM Navier-Stokes solver that includes the consideration of free-surface non-linearity through the SURF scheme. The computed sloshing forces and moments were fed into the time integration of the ship's motion, and the updated motion was, in turn, used as the excitation force for liquid sloshing, which is repeated for the ensuing time steps. For comparison, a sloshing motion coupled analysis program based on linear potential theory in the frequency domain was developed. The computer programs that were developed were applied to the side-by-side offloading operation between the offshore plant and the LNGC. The frequency-domain results reproduced the coupling effects qualitatively, but, in general, the peaks were over-predicted compared to experimental and time-domain results. The interactive effects between the sloshing liquid and the motion of the vessel can be intensified further in the case of multiple floating bodies.

• Journal of the Korean Society of Safety
• /
• v.17 no.1
• /
• pp.99-104
• /
• 2002

DQM(differential quadrature method) is applied to computation of eigenvalues of the equations of motion governing the free in-plane vibration for circular curved beams including mid-surface extension and the effects of rotatory inertia. Fundamental frequencies are calculated for the members with various end conditions and opening angles. The results are compared with numerical solutions by other methods for cases in which they are available. The differential quadrature method gives good accuracy even when only a limited number of grid points is used.