• Proceedings of the Korea Water Resources Association Conference
• /
• 2004.05b
• /
• pp.1376-1380
• /
• 2004

본 연구에서는 홍수범람지도 제작을 위해 사용되는 수치모형의 검증을 수행하고자 한다. 검증할 수치모형은 스위스의 Beffa에 의해 개발된 FLUMEN(FLUvial Modeling ENgine)으로서 수심 적분된 2차원 비선형 천수방정식에 불규칙 심각망을 이용한 유한체적법(finite volume method)이 적용된 수치모형이며 스위스, 독일, 오스트리아 등에서 홍수범람해석에 사용된 바 있는 모델이다. FLUMEN 모형의 검증을 위하여 범람해석시 가장 중요한 문제인 이동경계조건(moving boundary condition)을 포함하고 있는 원형섬에서의 고립파에 처오름높이를 계산하여 수리모형실험 결과와 비교한다. 수리모형실험은 미국 육군 공병단 산하의 해안공학연구소(CERC, Coastal Engineering Research Center)에서 수행되었으며(Liu 등, 1995) 수조의 중앙에 높이 0.625m, 해저지름 7.2m, 경사각 $14.04^{\circ}$의 원형섬이 위치한다. 본 연구 결과, 수치해석으로 계산된 섬에서의 실제 처오름높이와 입사파의 파고의 비(R/A)는 수리모형실험의 결과와 어느 정도 일치하였다.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.11
• /
• pp.1075-1080
• /
• 2010

For the aeroelastic analysis of a wing with friction damping, coupled time integration method was used to obtain time responses in the subsonic and transonic regions. To take into account aerodynamic nonlinearity induced by shock wave on the lifting surface, transonic small disturbance equation with in-phase periodic boundary condition was used for unsteady aerodynamic calculation. For 2-DOF airfoil system with displace-dependent friction dampers, the effects of normal load slope and Mach number on flutter boundary were investigated.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2001.11b
• /
• pp.1317-1322
• /
• 2001

What changes in the eigen values and eigen functions are produced if the boundary surface S is no longer rigid but has a specific acoustic admittance which may vary from point to point on S. In this paper, changes in eigen values and eigen functions are derived by using Kirchhoff-Helmholtz integral equation. And acoustic potential energy, which is representative measure describing the physical quantity in cavity, is defined. Acoustic potential energy can be divided into primary one and secondary one. Primary one is the acoustic potential energy through unchanged eigen functions, and secondary one is through changed eigen functions. Using these two term, we can find the eigenvalue problem, which gives the control performance when the boundary condition is changed.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.1
• /
• pp.61-68
• /
• 1991

The system such as railway bridge can be modelled as the restrained beam with intermediate supports. This kind of structures are subject to the moving load, which has a great effect on dynamic stresses and can cause sever motions, especially at high velocities. Therefore, to analyze the dynamic characteristics of the system due to the moving load is very important. In this paper, the governing equation of motion of a restrained beam subjected to the moving load is derived by using the Hamilton's principle. The orthogonal polynomial functions, which are trial functions and satisfying the geometric and dynamic boundary conditions, are obtained through simple procedure. The dynamic response of the system subjected to the moving loads is obtained by using the Galerkin's method and the numerical time integration technique. The numerical tests for various constraint, velocity and boundary conditions were preformed. Furthermore, the effects of mass of the moving load are studied in detail.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.6
• /
• pp.1478-1485
• /
• 1993

A shape optimization problem is formulated to determine the optimal position of heating lines in a compression molding die. The objective of the problem is that the cavity surface would be maintained by a prescribed uniform temperature. A boundary integral equation for the sensitivity of the temperature in terms of hole position is derived using the method of shape design sensitivity analysis. The boundary element method is employed to analyze the temperature and sensitivity field of the die. The sensitivity calculation algorithm is incorporated in an optimization routine. To demonstrate a numerical implementation, an example problem arising in thermal design of a compression molding die is dealt with, showing that the number of heating lines chosen for the die strongly affects the ultimate uniformity of the cavity surface temperature.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.5
• /
• pp.1084-1094
• /
• 1993

An efficient method for determining forced responses of a general linear structural system in time domain using subtructure modes and Lagrange multipliers is presented. Compared with the conventional mode synthesis methods, the suggested method does not construct the equations of motion of the combined whole structure and thus the modal parameters of the whole structure are not required. Only modal parameters of each substructure and geometric compatibility conditions are needed. Both the loaded interface free-free modes and free interface modes can be employed as the modal bases of each substructure. Recurrence discrete-time state equations based upon state transition matrix are formulated for the transient analysis of a parameter-changing system. It is shown form numerical examples that the suggested method is very accurate and efficient to calculate transient responses compares with the direct numerical integration method.

• 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 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.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.3
• /
• pp.725-733
• /
• 1990

An analysis of laminar mixed convection flow adjacent to the inclined flat surface which is subjected to a uniform temperature in a uniform free stream is performed. Nonsimilar boundary layed equation are derived by using the mixed convection parameters such that smooth transition from the purely forced convection limit to the purely free convection limit is possible. The governing equations are solved by a finite difference method using the coupled box scheme of sixth order. Numerical results are presented for prandtl numbers of 0.7 and 7 with the angle of inclination ranging from 0 to 90 degree from the vertical. The velocity distributions for the buoyancy assisting flow exhibit a significant overshoot above the free stream value in the region of intense mixed convection and the velocity field is found to be more sensitive to the buoyancy effect than the temperature field. The separation point near the wall was obtained for the buoyancy opposing flow. The local Nusselt number increases for buoyancy assisting flow and decreases for opposing flow with increasing value of the local Grashoff number in the mixed convection parameter. For large Prandtl number, the Nusselt number and the friction factor decrease significantly near the separation point. Present numerical predictions are in good agreement with recent experimental results by Ramachandran.

• Journal of the Society of Naval Architects of Korea
• /
• v.33 no.3
• /
• pp.35-47
• /
• 1996

A surface panel method treating a boundary-value problem of the Dirichlet type is presented to design a three dimensional body with free surface corresponding to a prescribed pressure distribution. An integral equation is derived from Green's theorem, giving a relation between total potential of known strength and the unknown local flux. Upon discretization, a system of linear simultaneous equations is formed including free surface boundary condition and is solved for an assumed geometry. The pseudo local flux, present due to the incorrect positioning of the assumed geometry, plays a role f the geometry corrector, with which the new geometry is computed for the next iteration. Sample designs for submerged spheroids and Wigley hull and carried out to demonstrate the stable convergence, the effectiveness and the robustness of the method. For the calculation of the wave resistance, normal dipoles and Rankine sources are distributed on the body surface and Rankine sources on the free surface. The free surface boundary condition is linearized with respect to the oncoming flow. Four-points upwind finite difference scheme is used to compute the free surface boundary condition. A hyperboloidal panel is adopted to represent the hull surface, which can compensate the defects of the low-order panel method. The design of a 5500TEU container carrier is performed with respect to reduction of the wave resistance. To reduce the wave resistance, calculated pressure on the hull surface is modified to have the lower fluctuation, and is applied as a Dirichlet type dynamic boundary condition on the hull surface. The designed hull form is verified to have the lower wave resistance than the initial one not only by computation but by experiment.