• Journal of computational fluids engineering
• /
• v.15 no.1
• /
• pp.71-80
• /
• 2010

The violent free-surface motions and the corresponding impact loads are numerically simulated by using the Moving Particle Semi-implicit (MPS) method, which was originally proposed by Koshizuka and Oka (1996) for incompressible flows. In the original MPS method, there were several shortcoming including non-optimal source term, gradient and collision models, and search of free-surface particles, which led to less-accurate fluid motions and non-physical pressure fluctuations. In the present study, how those defects can be remedied is illustrated by step-by-step improvements in respective processes of the revised MPS method. The improvement of each step is explained and numerically demonstrated. The numerical results are also compared with the experimental results of Martin and Moyce (1952) for dam-breaking problem. The current numerical results for violent free-surface motions and impact pressures are in good agreement with their experimental data.

• Tribology and Lubricants
• /
• v.33 no.6
• /
• pp.309-314
• /
• 2017

This paper presents a theoretical investigation of the dynamic characteristics of externally pressurized air pad bearings with closed loop grooves. These grooves are made on the surface of bearings to reduce the number of supply holes so that manufacturing costs can be reduced. The semi-implicit method is applied to calculate the time varying pressure profile on the air bearing surface owing to the advantages of numerical stability and fast time tracing characteristics. The static pressure of the groove bearings is much higher than that without grooves, so the groove bearings can provide high load carrying capacity. The equation of motion considering vertical motion and tilting motion are also solved using the Runge-Kutta 4th order method. By combining the semi-implicit method and the Runge-Kutta method, fast calculations of the dynamic behavior of the air bearing can be achieved. The variations of bearing reaction force, air film reaction moment, height, and tilting angle are investigated for the step force input, which is 20% higher than the bearing reaction, when the nominal clearance is 6 mm. The effect of the groove width and the groove depth are investigated by calculating the dynamic behavior. The possibility of the air hammering with the depth of the groove is found and discussed.

• Structural Engineering and Mechanics
• /
• v.3 no.2
• /
• pp.135-144
• /
• 1995

Since structural systems may fail in any one of several failure modes, computation of system reliability is always difficult. A method using numerical quadrature for computing structural system reliability with either one or more than one failure mode is presented in this paper. Statistically correlated safety margin equations are transformed into a group of uncorrelated variables and the joint density function of these uncorrelated variables can be generated by using the Maximum Entropy Method. Structural system reliability is then obtained by integrating the joint density function with the transformed safety domain enclosed within a set of linear equations. The Gaussian numerical integration method is introduced in order to improve computational accuracy. This method can be used to evaluate structural system reliability for Gaussian or non-Gaussian variables with either linear or nonlinear safety boundaries. It is also valid for implicit safety margins such as computer programs. Both the theory and the examples show that this method is simple in concept and easy to implement.

• Journal of the Society of Naval Architects of Korea
• /
• v.37 no.1
• /
• pp.67-81
• /
• 2000

An accurate and efficient numerical method for two-dimensional nonlinear radiation problem has been developed. The wave motion due to a moving body is described by the assumption of ideal fluid flow, and the governing Laplace equation can be effectively solved by the higher-order boundary element method with the help of the GMRES (Generalized Minimal RESidual) algorithm. The intersection or corner problem is resolved by utilizing the so-called discontinuous elements. The implicit trapezoidal rule is used in updating solutions at new time steps by considering stability and accuracy. Traveling waves caused by the oscillating body are absorbed downstream by the damping zone technique. It is demonstrated that the present method for time marching and radiation condition works efficiently for nonlinear radiation problem. To avoid the numerical instability enhanced by the local gathering of grid points, the regriding technique is employed so that all the grids on the free surface may be distributed with an equal distance. This makes it possible to reduce time interval and improve numerical stability. Special attention is paid to the local flow around the body during time integration. The nonlinear radiation force is calculated by the "acceleration potential technique". Present results show good agreement with other numerical computations and experiments.

• Journal of Korea Game Society
• /
• v.18 no.5
• /
• pp.123-132
• /
• 2018

Various methods have been proposed to efficiently animate the motion of soft objects in realtime. In order to maintain the topology between the elements of the objects, it is required to employ constraint forces, which limit the size of the time steps for the numerical integration and reduce the efficiency. To tackle this, an implicit method with larger steps was proposed. However, the method is, in essence, a linear system with a large matrix, of which solution requires heavy computations. Several approximate methods have been proposed, but the approximation is obtained with an increased damping and the loss of accuracy. In this paper, new integration method based on harmonic oscillation with better stability was proposed, and it was further stabilized with the hybridization with approximate implicit method. GPU parallelism can be easily implemented for the method, and large-scale soft objects can be simulated in realtime.

• Journal of applied mathematics & informatics
• /
• v.29 no.3_4
• /
• pp.937-948
• /
• 2011

The method of Asymptotic Inner Boundary Condition for Singularly Perturbed Two-Point Boundary value Problems is presented. By using a terminal point, the original second order problem is divided in to two problems namely inner region and outer region problems. The original problem is replaced by an asymptotically equivalent first order problem and using the stretching transformation, the asymptotic inner condition in implicit form at the terminal point is determined from the reduced equation of the original second order problem. The modified inner region problem, using the transformation with implicit boundary conditions is solved and produces a condition for the outer region problem. We used Chawla's fourth order method to solve both the inner and outer region problems. The proposed method is iterative on the terminal point. Some numerical examples are solved to demonstrate the applicability of the method.

• Journal of Ocean Engineering and Technology
• /
• v.21 no.6
• /
• pp.42-46
• /
• 2007

The impact of a single wave generated by a dam break with a tall structure is modeled with a three-dimensional version of the Moving particle semi-implicit (MPS) method. The particle method is more feasible and effective than methods based on grid connection problems involving violent free surface motions. In the present study, the Tsunami impact load and the change of longitudinal velocity component around the structure, which are obtained from the numerical simulation, are compared to those from experiments.

• Journal of the Korean Society of Surveying, Geodesy, Photogrammetry and Cartography
• /
• v.8 no.2
• /
• pp.9-16
• /
• 1990

This is a geodetic study on the the shoreline changes occurred by the facilities constructed in the beach. There are emperical, hydrological and numerical methods in predicting of the shorline changes. Numerical method is the most suitable method in the field of geodesy. There are many predicting models. This study adopted one-line model because it has a few hydrological factor and simplify the natural phenomena. This study established the ideal seawall boundary condition, applied the explicit model and the implicit model in the Dongsan harbour in East Sea, and could predict the optimum seawall position for protection of shore. The results are following ; Seawall protect shore of which input angle of wave is below 20$^\circ$, a ratio of wave height bleak/line does not effect in shoreline changes. The implicit model is accuracy but can not predict longtime change. But the explicit model is the opposite of the implicit model.

• International Journal of Ocean System Engineering
• /
• v.1 no.1
• /
• pp.37-45
• /
• 2011

The nonlinear free-surface motions interacting with a floating body were investigated using the Moving Particle Semi-implicit (MPS) method proposed by Koshizuka and Oka [6] for incompressible flow. In the numerical method, more realistic Lagrangian moving particles were used for solving the flow field instead of the Eulerian approach with a grid system. Therefore, the convection terms and time derivatives in the Navier-Stokes equation can be calculated more directly, without any numerical diffusion, instabilities, or topological failure. The MPS method was applied to a numerical simulation of predicting the efficiency of floating-type breakwater interacting with waves.

• Structural Engineering and Mechanics
• /
• v.48 no.5
• /
• pp.711-736
• /
• 2013

In this paper new implicit time integration called N-IHOA is presented for dynamic analysis of high damping systems. Here, current displacement and velocity are assumed to be functions of the velocities and accelerations of several previous time steps, respectively. This definition causes that only one set of weighted factors is calculated from the Taylor series expansion which leads to a simple approach and reduce the computational efforts. Moreover a comprehensive study on stability of the proposed method i.e., N-IHOA compared with IHOA integration which is performed based on amplification matrices proves the ability of the N-IHOA in high damping vibrations such as control systems. Also, wide range of numerical examples which contains single/multi degrees of freedom, damped/un-damped, free/forced vibrations from finite element/finite difference demonstrate that the accuracy and efficiency of the proposed time integration is more than the common approaches such as the IHOA, the Wilson-${\theta}$ and the Newmark-${\beta}$.