• 대한기계학회논문집B
• /
• 제24권10호
• /
• pp.1317-1325
• /
• 2000

A new efficient numerical method for computing three-dimensional, unsteady, incompressible flows is presented. To eliminate the restriction of CFL condition, a fully-implicit time advancement in which the Crank-Nicolson method is used for both the diffusion and convection terms, is adopted. Based on an approximate block LU decomposition method, the velocity -pressure decoupling is achieved. The additional decoupling of the intermediate velocity components in the convection term is made for the fully -implicit time advancement scheme. Since the iterative procedures for the momentum equations are not required, the velocity components decouplings bring forth the reduction of computational cost. The second-order accuracy in time of the present numerical algorithm is ascertained by computing decaying vortices. The present decoupling method is applied to minimal channel flow unit with DNS (Direct Numerical Simulation).

• Journal of Mechanical Science and Technology
• /
• 제19권spc1호
• /
• pp.481-486
• /
• 2005

This paper introduces a numerical method for integration of the linear and nonlinear differential dynamic equation of motion. The variation of acceleration in two time steps is approximated as a combination of both trigonometric cosine and hyperbolic cosine functions with weighted coefficient. From which all necessary formulae are elaborated for the direct integration of the governing equation. A number of linear and nonlinear dynamic problems with various degrees of freedom are analysed using both the suggested method and Newmark method for the comparison. The numerical results show high advantages and effectiveness of the new method.

• Ocean Systems Engineering
• /
• 제7권4호
• /
• pp.399-412
• /
• 2017

The Gauss-Legendre integral method is applied to numerically evaluate the Green function and its derivatives in finite water depth. In this method, the singular point of the function in the traditional integral equation can be avoided. Moreover, based on the improved Gauss-Laguerre integral method proposed in the previous research, a new methodology is developed through the Gauss-Legendre integral. Using this new methodology, the Green function with the field and source points near the water surface can be obtained, which is less mentioned in the previous research. The accuracy and efficiency of this new method is investigated. The numerical results using a Gauss-Legendre integral method show good agreements with other numerical results of direct calculations and series form in the far field. Furthermore, the cases with the field and source points near the water surface are also considered. Considering the computational efficiency, the method using the Gauss-Legendre integral proposed in this paper could obtain the accurate numerical results of the Green function and its derivatives in finite water depth and can be adopted in the near field.

• 한국전산유체공학회지
• /
• 제6권4호
• /
• pp.26-34
• /
• 2001

The flow past a circular cylinder forced to vibrate transversely is numerically simulated by solving the two-dimensional Navier-Stokes equations modified by the vibration velocity of a circular cylinder at a Reynolds number of 164. The higher-order finite difference scheme is employed for the spatial discretization along with the second order Adams-Bashforth and the first order backward-Euler time integration. The calculated cylinder vibration frequency is between 0.60 and 1.30 times of the natural vortex-shedding frequency. The calculated oscillation amplitude extends to 25% of the cylinder diameter and in the case of the lock-in region it is 60%. It is made clear that the cylinder oscillation has influence on the wake pattern, the time histories of the drag and lift forces, power spectral density and phase diagrams, etc. It is found that these results include both the periodic (lock-in) and the quasi-periodic (non-lock-in) state. The vortex shedding frequency equals the driving frequency in the lock-in region but is independent in the non-lock-in region. The mean drag and the maximum lift coefficient increase with the increase of the forcing amplitude in the lock-in state. The lock-in boundaries are also established from the present direct numerical simulation.

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

The flow past a circular cylinder forced to vibrate transversely is numerically simulated by solving the two-dimensional Wavier-Stokes equations modified by the vibration velocity of a circular cylinder at a Reynolds number of 164. The higher-order finite difference scheme is employed for the spatial discretization along with the second order Adams-Bashforth and the first order backward-Euler time integration. The calculated cylinder vibration frequency is between 0.60 and 1.30 times of the natural vortex-shedding frequency. The calculated oscillation amplitude extends to $25\%$ of the cylinder diameter and in the case of the lock-in region it is $60\%$. It is made clear that the cylinder oscillation has influence on the wake pattern, the time histories of the drag and lift forces, power spectral density and phase diagrams, etc. It is found that these results include both the periodic (lock-in) and the quasi-periodic (non-lock-in) state. The vortex shedding frequency equals the driving frequency in the lock-in region but is independent in the non-lock-in region. The mean drag and the maximum lift coefficient increase with the increase of the forcing amplitude in the lock-in state. The lock-in boundaries are also established from the present direct numerical simulation.

• 한국정밀공학회:학술대회논문집
• /
• 한국정밀공학회 1997년도 춘계학술대회 논문집
• /
• pp.552-556
• /
• 1997

The Direct Metal Prototyping(DMP), one of the rapid prototyping technologies, allows the manufacturing of three-dimensional metallic parts using metal powders directly from the CAD data. Laser power and scanning speed are the most important variables of the process. The objective of this study is to obtain the design data for laser power and scanning speed to bond metal powders effectively using the finite element method. To obtain the design values, a numerical analysis considering two-dimensional heat transfer during the sintering of metal powder layers of the process was performed. The laser beam has been modeled to have directionality in its heat flux distribution, i. e., in the scanning direction a Gaussian beam mode distribution has been assumed and in the thickness direction a square beam mode distribution. The three-dimensional irregular distribution of metal powders of the powder layer is idealized as two-dimensional distribution in which metal powders are located regularly and periodically on the plate. In this study the design values of laser power vs scanning speed have been obtained. Temperature distribution and temperature variation of the powder layers with respect to time have been predicted. The commputed dsign data will be useful in determining the initial conditions of the process.

• Journal of applied mathematics & informatics
• /
• 제32권3_4호
• /
• pp.389-403
• /
• 2014

In this paper, we first provide comparison results of preconditioned AOR methods with direct preconditioners $I+{\beta}L$, $I+{\beta}U$ and $I+{\beta}(L+U)$ for solving a linear system whose coefficient matrix is a large sparse irreducible L-matrix, where ${\beta}$ > 0. Next we propose how to find a near optimal parameter ${\beta}$ for which Krylov subspace method with these direct preconditioners performs nearly best. Lastly numerical experiments are provided to compare the performance of preconditioned iterative methods and to illustrate the theoretical results.

• 대한기계학회논문집B
• /
• 제25권3호
• /
• pp.389-396
• /
• 2001

Direct numerical simulation of bubble growth and merger process on a single nucleation site during partial nucleate boiling is performed. The equations governing conservation of mass, momentum and energy are solved using a finite difference method combined with a level set method for capturing the vapor-liquid interface. The level set method is modified to include the effects of phase change at the interface and contact angle at the wall. Also, a simplified formulation for predicting the evaporative heat flux in a thin liquid micro-layer is developed and incorporated into the level set formulation. Based on the numerical results, the bubble growth and merger pattern and its effect on the heat transfer are discussed.

• 한국지반공학회논문집
• /
• 제33권12호
• /
• pp.93-105
• /
• 2017

본 연구에서는 하수관 손상에 의한 지반함몰 발생 과정에서 지반의 포화도 상승에 따른 흙의 불포화 강도 저하의 영향을 파악하기 위하여, 직접 전단 실험, 모형 실험, 그리고 수치해석을 수행하였다. 직접 전단 시험 결과, 흙의 마찰각은 포화도의 영향을 크게 받지 않으나, 점착력은 포화도의 영향을 크게 받음을 알 수 있다. 포화도 상승에 따른 강도저하의 영향만을 고려하기 위하여, 물의 침투효과를 배제한 모형 실험을 실시하여 지반 함몰 현상을 재현하였다. 지반 함몰은 대변위를 동반하며, 기존 유한요소법의 적용이 어렵다. 본 연구에서는 대변위 해석 기법인 일반 보간 재료점법을 사용하여 수행한 모형실험을 수치적으로 모사하였다. 비록 경계 조건 차이, 불완전 포화, 손상부 연결관의 폐색등에 의해 함몰 시간 등에는 차이가 있었지만, 유사한 토체의 변형 거동이 모형 실험과 수치해석에서 발견되었다.

• 한국통신학회논문지
• /
• 제11권6호
• /
• pp.396-401
• /
• 1986

Steepest Descent Method를 정전계 문제의 해석에 적용하였다. 반복법과 moment method와 같은 direct method와의 차이점을 기술하였다. 반복법을 이용하여 정확한 해에 단조적으로 수렴하며 대형 system의 해석에 적합함을 보였다. 정전계 문제에 대해서 해석한 수치결과를 제시하였으며 moment solution과 잘 일치하고 있음을 보였다.