• 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 the Computational Structural Engineering Institute of Korea
• /
• v.23 no.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
• /
• v.1 no.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.

• The Proceeding of the Korean Institute of Electromagnetic Engineering and Science
• /
• v.5 no.3
• /
• pp.14-21
• /
• 1994

In this paper, a numerical method to be obtain the radar cross section(RCS) of three- dimensional bodies with arbitrary geometry and material compositions on the electromagnetic field with arbitrary incident angle is described. The RCS is obtained by solving the individual surface integral equation about multilayers scatterer using the three-dimensional bonudary element method(BEM). To show propriety and usefulness as to the three-dimensional BEM in this paper, the choice of a geometry is a multi-regular hexahedron and multi-right-angled hexahedron out of oblique incident electric field, and is considered to apply to every condition with loss sufficiently.

• Annual Conference of KIPS
• /
• 2005.05a
• /
• pp.785-788
• /
• 2005

일정한 형태의 물체를 분석하기 위해 사용되는 각-거리 그래프를 이용하여 임의의 물체의 경계선 내부 영역의 면적을 측정하는 영상처리 알고리듬을 제안한다. 물체의 경계선 내부의 한 점을 중심으로 1차 각-거리 그래프를 생성하고 이 그래프로부터 거리 값이 급격히 변화하는 위치를 추출하여 1 차 그래프에서 접근하지 못한 영역을 인식하여 새 영역에서의 한 점을 중심으로 2 차 각-거리 그래프를 생성한다. 이와 같이 계층적으로 구성된 각-거리 그래프 그룹에 대해 거리의 제곱을 각도 방향으로 적분하여 물체의 경계선 내부 영역의 면적을 측정한다.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1997.04a
• /
• pp.608-614
• /
• 1997

2차원 압축성 오일러 방정식을 지배 방정식으로 하여 음향학적 문제를 직접 수치 계산하였다. Roe의 개략적인 리만해(Approximate Riemann Solver)에 의거하여 대류항을 구성하였고 높은 정확도를 유지하기 위해 Vanleer의 MUSCL 방법을 사용하여 공간을 이산화하였다. 시간 적분으로는 2차의 Runge-Kutta 방법을 사용하였다. Resonance tube 내의 음향학적 문제 해석을 위해 최근까지 제시된 음향 경계 조건들을 비교하였고 각 경계 조건들이 갖는 물리적 의미를 고찰하였다. 이들을 통하여 압력파의 물리적인 행태를 모사하였고 실험과 비교하였다. 계산결과들을 통하여 사용된 방법은 전형적인 압축성 유동 해석 문제 뿐만 아니라 선형, 비선형의 음향학적 문제들을 적절한 경계 조건과 결합하여 해석할 수 있음을 알 수 있었고 resonance tube와 같은 공학적으로 실제적인 문제들에 대한 응용이 기대된다.

• Communications of Mathematical Education
• /
• v.18 no.1 s.18
• /
• pp.239-248
• /
• 2004

본 논문에서는 맨델브로트(Mandelbrot) 집합의 개념을 3차의 복소 다항식 z^3$+c 에 확장시켜 3차 분기집합을 정의하고, 이 집합의 2-주기 성분의 경계선 방정식과 관련 기하학적 성질을 고등학교 및 대학에서 다루는 미적분학 관점에서 분석하고자 한다. 복소수, 삼각함수, 매개함수, 함수의 극값, 미분 및 적분 등의 기초 이론을 활용하여 2-주기 성분의 경계선 방정식을 매개함수로 표시하고, 경계선의 내부 면적, 둘레 길이, 무게중심 등을 이론적으로 기술한다. 수학 소프트웨어인 매스매티카(Mathematica)를 활용하여 2-주기성분의 작도 및 기하학적 성질에 관한 수치 해석적 결과를 제시하고자 한다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2006.05a
• /
• pp.1910-1913
• /
• 2006

To solve the interaction of incident monochromatic waves with a bottom-fixed vertical circular cylinder, a numerical analysis by boundary element method is developed using three-dimensional linear potential theory. A numerical analysis by boundary element method is based on Green's theorem and introduce to an integral equation for the fluid velocity potential around the vertical circular cylinder. These numerical results are compared with those of ManCamy and Fuchs(1954) and Williams and Mansour(2002), and it has shown good relationship with their results. This numerical analysis developed by boundary element method will be applied for various offshore structures to be constructed in coastal zones in the future.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.2
• /
• pp.146-158
• /
• 1991

The flow around a two-dimensional foil section Is measured by a LDV(Laser Doppler Velocimetry) system which is capable of measuring the datailed flow field without interfering the original flow field. A 2-color 3-beam LDV system, which is capable of mea,;tiring 2 velocity components simultaneously and uses 2W Ar-Ion laser source, is used to measure the flow field around an NACA0012 foil section. The measured flow velocities are analysed iii order to study the boundary layer characteristics, flow separation and the detail structure of the flow near the trailing edge of the foil. The boundary layer characteristics are compared with the results by the head's momentum integral method. For the case of small angle of attack at relatively higher Reynolds number, both results show good agreements. The measured data of the velocity field around an NACA0012 foil section would be valuable data to validate the CFD(Computational Fluid Dynamic) calculation results. The developed experimental technique to evaluate the characteristics of two-dimensional foil sections is essential tool to develope new blade sections which have good lift characteristics and better cavitation performances.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.22 no.1
• /
• pp.36-44
• /
• 2018

The thermal response of carbon/phenolic used in a solid rocket nozzle liner was analyzed. In this paper, the numerical analysis of the thermal response of carbon/phenolic consists of (1) the integration equation of the boundary layer to obtain the convective heat transfer coefficient of the combustion gas on the rocket nozzle wall and (2) 1-D finite difference method for heat conduction of carbon/phenolic to calculate the ablation, char, and temperature. The calculated result was compared with the result of a blast-tube-type test motor. It is found that the calculated result shows good agreement with the thermal response of the test motor, except at the vicinity of the throat insert.