• Water for future
• /
• 제21권3호
• /
• pp.299-307
• /
• 1988

A model is developed to predict the long-term beach evolution near the long groin considering the combined effects of variation of sea level, wave refraction and diffraction. A numerical solution for this problem is solved by considering the equation as a system subject to the boundary condition for longshore transport rate. One possible method is the centered Crank-Nicolson type implicit scheme. The results which ard obtained by applying this numerical model at Songdo beach, Pohang are as follows. Owing to the approximation used in the calculation of the refraction and diffraction coefficients, the discrepancy between the predicted and actual shoreline occurs to the interior of long groin. However, the shape of shoreline at the exterier of long groins agrees well.

• Bulletin of the Society of Naval Architects of Korea
• /
• 제27권3호
• /
• pp.19-30
• /
• 1990

This note describes an application of Hamiton's principle to nonlinear free-surface flow problems. Two functionals are constructed based on classical Hamilton's principle with a modification due to the presence of a free surface. As an effort towards the development of an efficient numerical scheme for our problem, we present the following three test results: i) The bounding principles of the eigenvalues for the linear dispersion relation. ii) By assuming steady solitary waves, an approximate relation between the amplitudes and the speeds of solitary waves are derived from the two functionals constructed. Their numerical results are compared with those of Longuet-Higgins & Fenton(1974). iii) The shapes and charicteristics of solitary waves are computed from two sets of functionals by varying the number of total finite elements in the fluid domain.

• Journal of the Society of Naval Architects of Korea
• /
• 제28권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.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• 제19권6호
• /
• pp.607-613
• /
• 2009

A new formulation of NDIF method to the algebraic eigenvalue problem is introduced to efficiently extract natural frequencies of arbitrarily shaped plates with the simply supported boundary condition. NDIF method, which was developed by the authors for the free vibration analysis of arbitrarily shaped membranes and plates, has the feature that it yields highly accurate natural frequencies compared with other analytical methods or numerical methods(FEM and BEM). However, NDIF method has the weak point that it needs the inefficient procedure of searching natural frequencies by plotting the values of the determinant of a system matrix in the frequency range of interest. A new formulation of NDIF method developed in the paper doesn't require the above inefficient procedure and natural frequencies can be efficiently obtained by solving the typical algebraic eigenvalue problem. Finally, the validity of the proposed method is shown in several case studies, which indicate that natural frequencies by the proposed method are very accurate compared to other exact, analytical, or numerical methods.

• Journal of the Society of Naval Architects of Korea
• /
• 제30권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 the Computational Structural Engineering Institute of Korea
• /
• 제27권3호
• /
• pp.173-182
• /
• 2014

The mixed convolved action principle provides a new rigorous weak variational formalism for a broad range of initial boundary value problems in mathematical physics and mechanics in terms of mixed formulation, convolution, and fractional calculus. In this paper, its potential in the development of numerical methods for transient problems in various dynamical systems when adopting temporally second order approximation is investigated. For this, the classical single-degree-of-freedom linear elastic dynamical systems are primarily considered to investigate computational characteristics of the developed algorithms. For the undamped system, all the developed algorithms are symplectic with respect to the time step. For the damped system, they are shown to be accurate with good convergence characteristics.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• 제23권3호
• /
• pp.331-340
• /
• 2010

The basic theory and application of new adaptive finite element algorithm have been proposed in this study including the adaptive hp-refinement strategy, and the effective method for constructing hp-approximation. The hp-adaptive finite element concept needs the integrals of Legendre shape function, nonuniform p-distribution, and suitable constraint of continuity in conjunction with irregular node connection. The continuity of hp-adaptive mesh is an important problem at the common boundary of element interface. To solve this problem, the constraint of continuity has been enforced at the common boundary using the connectivity mapping matrix. The effective method for constructing of the proposed algorithm has been developed by using hierarchical nature of the integrals of Legendre shape function. To verify the proposed algorithm, the problem of simple cantilever beam has been solved by the conventional h-refinement and p-refinement as well as the proposed hp-refinement. The result obtained by hp-refinement approach shows more rapid convergence rate than those by h-refinement and p-refinement schemes. It it noted that the proposed algorithm may be implemented efficiently in practice.

• Journal of the Earthquake Engineering Society of Korea
• /
• 제1권4호
• /
• pp.59-68
• /
• 1997

In analysis of underground structures, the effects of artificial boundary conditions are considered as one of the major reasons for differences from experimental results. These phenomena can be overcome by using the boundary elements which satisfy the multi-layered half space conditions. The fundamental solutions of multi-layered half-space for boundary element method is formulated satisfying the transmission and reflection of waves at each layer interface and radiation conditions at bottom layer. The governing equations can be obtained from the displacements at each layer which are expressed in terms of harmonic functions. All types of waves can be included using the complete response from semi-infinite integrals with respect to horizontal wavenumbers using expansion of Fourier series and Hankel transformation. Two dimensional Green's functions are derived from cylindrical Navier equations and potentials performing infinite integration in y-direction. In this case, it is effective to transform into two dimensional problem using semi-analytical integration and sinusoidal Bessel function. Some verifications are given to show the accuracy and efficiency of the developed method, and numerical examples to demonstrate the dynamic behavior of underground with various properties.

• Bulletin of the Society of Naval Architects of Korea
• /
• 제7권1호
• /
• pp.27-36
• /
• 1970

선박유체역학(船舶流體力學)의 경계조건(境界條件)의 선형화(線型化)를 위한 퍼터베이션(Perturbation)방법(方法)을 따라 배를 가늘고 길게 보느냐, 엷은 수평판(水平板)으로 보느냐 등등에 따라 여러가지 다른 종류(種類)의 수학적(數學的) 정의(定義)에 따르는 이론(理論)을 추궁함으로서, 1) Neumann Ship, 2) Thin Ship, 3) Michell Ship, 4) Slender Ship 등을 이론적(理論的)인 면(面)에서와 수치계산(數値計算) 예(例)로서 비교검토(比較檢討) 하였다. 동시(同時)에 이에 따르는 "직접문제(直接問題)" "간접문제(間接問題)" 등(等)에 관(關)하여도 언급(言及)하고 그 차이점(差異點)을 지적하였다. 이에 의(依)하면 정속전진시(定速前進時)에 L/B가 10인 경우에는 최대광폭의 차(差)가 $5{\sim}8%$정도 되고 또 선수수선반각(船首水線半角)에 상당한 차(差)가 있음으로 같은 이론적(理論的)인 계산치(計算値)에 대(對)하여 저항(抵抗)은 적어도 $10{\sim}15%$정도 차(差)가 생길 것이 아닌가 보게끔 되었다.

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

극초음속 유동장의 정확한 해석을 위해 AUSMPW+ 수치기법과 충격파 포착시 생기는 수치오차를 제거하기 위해 충격파 정렬 기법(Shock-Aligned Grid Technique)을 개발하였다. AUSMPW+ 수치기법은 자체 수치점성이 적은 수치기법으로 점성 경계층 계산시 정확한 계산결과를 보여주며 기존의 AUSM 계열이 가지는 문제점인 물성치의 진동 현상을 제거한 수치기법이다. 원통형과 무딘 물체 주위의 극초음속 유동장 해석을 통해 공력이 진동현상 없이 정확하게 계산됨을 확인하였다. 그리고 충격파 정렬 기법의 특성을 파악하기 위해 충격파 반사문제와 충격파-충격파 상호작용 문제를 해석하여 수치오차 없이 충격파를 포착할 수 있음을 보였다. 또한 화학적 평형 비평형 유동 영역까지 충격파 정렬 격자 기법을 확장하였다.