• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2004년도 춘계 학술대회논문집
• /
• 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
• /
• 제8권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.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제6권2호
• /
• pp.206-218
• /
• 2014

The wave attenuation by floating breakwaters in high amplitude waves, which can lead to wave overtopping and breaking, is examined by numerical simulations. The governing equations, the Navier-Stokes equations and the continuity equation, are calculated in a fixed Cartesian grid system. The body boundaries are defined by the line segment connecting the points where the grid line and body surface meet. No-slip and divergence free conditions are satisfied at the body boundary cell. The nonlinear waves near the moving body is defined using the modified marker-density method. To verify the present numerical method, vortex induced vibration on an elastically mounted cylinder and free roll decay are numerically simulated and the results are compared with those reported in the literature. Using the present numerical method, the wave attenuations by three kinds of floating breakwaters are simulated numerically in a regular wave to compare the performance.

• Wind and Structures
• /
• 제31권1호
• /
• pp.75-84
• /
• 2020

The aerostatic stability analysis of a long-span suspension bridge by the Element-free Galerkin (EFG) method is presented in this paper. Nonlinear effects due to wind structure interactions should be taken into account in determining the aerostatic behavior of long-span suspension bridges. The EFG method is applied to investigate torsional divergence of suspension bridges, based on both the three components of wind loads and nonlinearities of structural geometric. Since EFG methods, which are based on moving least-square (MLS) interpolation, require only nodal data, the description of the geometry of bridge structure and boundaries consist of defining a set of nodes. A numerical example involving the three-dimensional EFG model of a suspension bridge with a span length of 888m is presented to illustrate the performance and potential of this method. The results indicate that presented method can effectively be applied for modeling suspension bridge structure and the computed results obtained using present modeling strategy for nonlinear suspension bridge structure under wind flow are encouragingly acceptable.

• 한국지하수토양환경학회:학술대회논문집
• /
• 한국지하수토양환경학회 2001년도 총회 및 춘계학술발표회
• /
• pp.77-80
• /
• 2001

The element-free Galerkin (EFG) method is one of meshless methods, which is an efficient method of modeling problems of fluid or solid mechanics with complex boundary shapes and large changes in boundary conditions. This paper discusses the theory of the EFG method and its applications to modeling of groundwater flow. In the EFG method, shape functions are constructed based on the moving least square (MLS) approximation, which requires only set of nodes. The EFG method can eliminate time-consuming mesh generation procedure with irregular shaped boundaries because it does not require any elements. The coupled EFG-FEM technique was introduced to treat Dirichlet boundary conditions. A computer code EFGG was developed and tested for the problems of steady-state and transient groundwater flow in homogeneous or heterogeneous aquifers. The accuracy of solutions by the EFG method was similar to that by the FEM. The EFG method has the advantages in convenient node generation and flexible boundary condition implementation.

• 한국정밀공학회:학술대회논문집
• /
• 한국정밀공학회 1996년도 추계학술대회 논문집
• /
• pp.743-747
• /
• 1996

A Potential Field Method is applied to the proposed algorithm for the planning of collision-free paths of redundant manipulators. The planning is carried out on the base of kinematic configuration. To make repulsive potentials, sources are distributed on the boundaries of obstacles. To escape from local minimum of the main potential and to attack other difficulties of the planning, various potentials are defined simultaneously, Inverse Kinematics Problems of the redundant manipulators are solved by unconstrained optimization method. Computer simulation result of the path planning is presented.

• 한국지진공학회:학술대회논문집
• /
• 한국지진공학회 1997년도 춘계 학술발표회 논문집 Proceedings of EESK Conference-Fall 1997
• /
• pp.128-135
• /
• 1997

The indirect boundary integral equation is formulated to analyze the behavior of a cavity in a multilayered half-plane subjected to earthquake waves. This formulation uses the fundamental solutions that are numerically calculated by the generalized transmission and reflection coefficient method. The free surface of the cavity without external excitation influences the behavior of the half-plane. Consequently this analysis adds the consideration of scattering-field into the analysis and the total motion field of the cavity is decomposed into the free-field and scattering-field motions. The free-field motion is obtained from the modification of the transmission and reflection coefficient method. The scattering-field motion is calculated is calculated by the indirect boundary value problem which has the ficticious boundaries and sources. In this study, P wave, SV wave, SH wave, and Rayleigh wave are analyzed respectively.

• 한국재료학회지
• /
• 제24권11호
• /
• pp.585-590
• /
• 2014

Microstructural evolution in the thickness direction of an oxygen free copper processed by accumulative rollbonding (ARB) is investigated by electron back scatter diffraction (EBSD) measurement. For the ARB, two copper alloy sheets 1 mm thick, 30 mm wide and 300 mm long are first degreased and wire-brushed for sound bonding. The sheets are then stacked and roll-bonded by about 50% reduction rolling without lubrication at an ambient temperature. The bonded sheet is then cut to the two pieces of the same dimensions and the same procedure was repeated on the sheets up to eight cycles. The specimen after 1 cycle showed inhomogeneous microstructure in the thickness direction so that the grains near the surface were finer than those near the center. This inhomogeneity decreased with an increasing number of ARB cycles, and the grain sizes of the specimens after 3 cycles were almost identical. In addition, the aspect ratio of the grains decreased with an increasing number of ARB cycles due to the subdivision of the grains by shear deformation. The fraction of grains with high angle grain boundaries also increased with continuing process of the ARB so that it was higher than that of the low angle grain boundaries in specimens after 3 cycles. A discontinuous dynamic recrystallization occurred partially in specimens after 5 cycles.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제19권1호
• /
• pp.1-21
• /
• 2015

This paper presents two algorithms based on the Jamshidian equation which is from the Black-Scholes partial differential equation. The first algorithm is for American call options and the second one is for American put options. They compute numerically free boundary and then option price, iteratively, because the free boundary and the option price are coupled implicitly. By the upwind finite-difference scheme, we discretize the Jamshidian equation with respect to asset variable s and set up a linear system whose solution is an approximation to the option value. Using the property that the coefficient matrix of this linear system is an M-matrix, we prove several theorems in order to formulate a bisection method, which generates a sequence of intervals converging to the fixed interval containing the free boundary value with error bound h. These algorithms have the accuracy of O(k + h), where k and h are step sizes of variables t and s, respectively. We prove that they are unconditionally stable. We applied our algorithms for a series of numerical experiments and compared them with other algorithms. Our algorithms are efficient and applicable to options with such constraints as r > d, $r{\leq}d$, long-time or short-time maturity T.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 춘계 학술대회논문집
• /
• pp.55-60
• /
• 1998

An elliptic grid generation scheme using Laplace's equations guarantees the resulting grids to be crossing-free as a result of maximum principle in its analytic form. Numerical results, however, often show the grid lines overlapping each other or crossing the boundaries, especially for very sharp convex corners. The cause of this problem is investigated, and it is found that this problem can be handled by properly modifying the coefficients of transformed Laplace's equations in the computational domain.