• Magazine of the Korean Society of Agricultural Engineers
• /
• v.37 no.2
• /
• pp.53-63
• /
• 1995

본 논문의 다경간 연속아치의 자유진동에 관한 연구이다. 다경간 연속아치의 고유진 동수 및 진ㄷㅇ형을 산출하기 위하여 내부지점의 지점조건에 다른 경계조건식을 유도하였다. 아치의 선형은 포물선을 택하였으며, 회전-로울러-회전, 고정-회전-고정의 지점 조건을 갖는 2경간 연속아치에 대한 수치해석 결과를 제시하였다. Runge-Kutta maethod을 이용 하였다. 실제 수치해석예에서는 회전관성이 고유진동수에 미치는 영향을 고찰 하였으며, 무차원 고유진동수와 아치높이 지간길이비 및 세장비 사이의 관계를 분석하였다. 또한 실험을 토아여 이론적인 해석결과를 검증하였다.

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.1
• /
• pp.29-42
• /
• 1992

A potential based panel method is used to predict the open water characteristics of spade-type rudders. The inflow velocity is assumed to be constant in lime and uniform in space. Source and dipole are distributed on the rudder surface. It is assumed that the wake surface is streaming from trailing edge and it is represented by dipole distribution. In this paper, wake geometry is assumed by imposing appropriate conditions at the trailing edge and far from the body. The effects of wake geometry are studied. The pressure Kutta condition is applied at the trailing edge, the effects of which are compared with those of two-dimensional Kutta condition. The results of calculations for a spade-type rudder are compared with published results. It is concluded that this approach shows fairly good agreement with experimental results and can be used in the initial design stage of a rudder.

• The Korean Journal of Rheology
• /
• v.8 no.3_4
• /
• pp.187-198
• /
• 1996

최근에 본 연구자에 의해서 단축 스크류 공정에서 카오스 스크류라고 명명되어진 카오스 혼합장치가 성공적으로 개발되었다. 기하학적 조건이나 공정조건에 대한 설계변수로 카오스 스크류를 설계하기 위하여 체류시간, 포인카레 단면 그리고 혼합패턴등에 대한 계산 과 해석이 이루어져야 하는데 이를 단지 Runge-Kutta 방법에 의해 속도장을 적분한다면 상당한 계산시간이 소비된다. 이러한 수치문제를 극복하기 위하여 본논문에서는 새로운 사 상법을 제안한다. 이 방법으　 사용하면 벽면 근처의 특이점 영역에서도 수치문제가 해결된 다. 본 논문에서 제안하는 수치사상법은 Runge-Kutta 방법에 비하여 수치계산의 효율성과 정확도 면에서, 특히 유안요소법으로 얻은 속도장에 대하여 우수함이 밝혀졌다. 이러한 사상 법은 공간주기 유동장뿐만 아니라 시간주기 유동장에서 적용할수 있다.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.40 no.2
• /
• pp.101-107
• /
• 2012

In this study, Constraint Force Equation Methodology (CFEM) is used to develop a multi-body dynamic analysis program with constraints. Seven constraint models are implemented to analyze constraint motions of multiple bodies. The augmented equations with the constraints are solved with the 4th order Runge-Kutta method for higher degree of accuracy. The analysis code is verified by comparing the analysis results of the motion of bodies with various constraints to published results.

• Journal of the Society of Naval Architects of Korea
• /
• v.34 no.4
• /
• pp.12-20
• /
• 1997

A potential-based panel method has been developed to investigate the resistance characteristics of a high speed catamaran advancing on the free surface. Normal dipoles and sources are distributed on the body surface while sources are distributed on the free surface. Linearised free surface conditions are used in the present analysis. To avoid the instabilities due to the velocity difference between inner and outer flow of a high speed catamaran, Kutta condition has been applied at the stern. Model test has been carried out not only to validate the numerical results but to confirm the capabilities of a CWC(Circulating Water Channel). It is believed that we can obtain the qualitatively reasonable results in the CWC. Computed results are compared with those of experiments and Insel's experimental values. Since the Kutta condition is applied at the stern, stable solutions are obtained at the high speed range. The present method, using linearised free surface conditions at the high speed range, seems to be a useful tool in the hull form design of a high speed catamaran.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.181-186
• /
• 1999

The effects of the frequency of upstream gust on the unsteady boundary characteristics on a downstream blade was simulated by using a Navier-Stokes code. The Navier-Stokes code is based on an unstructured finite volume method and uses a low Reynolds k-e turbulence model to close the momentum equations. The MIT flapping foil experiment set-up is used to simulate the interaction between the upstream wake and a blade. The frequency of the upstream wake is simulated by varying rate of pitching motion of the flapping airfoils. Three reduced frequencies. 3.62. 7.24. and 10.86. are simulated. As the frequency increases, the unsteady fluctuation on the surfaces of the downstream hydrofoil is shown to decrease while the upstream flapper wake has larger first harmonics of y-velocity component. The unsteady vortices are shown to interact with each other and. as a result. the upstream wake becomes undiscernible inside the inner layer. The turbulence kinetic energy shows a similar behavior. Limiting streamlines around the trailing edge of the flapper are shown to conform with the unsteady Kutta condition for a round trailing edge. while limiting streamlines around the trailing edge of the hydrofoil conforms with the unsteady Kutta condition for a sharp edge.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.3
• /
• pp.87-99
• /
• 1994

Reynolds Averaged Navier-Stokes equations are solved numerically for the computation of turbulent flow around a Wigley double model. A second order finite difference method is applied for the spatial discretization on the nonstaggered grid system and 4-stage Runge-Kutta scheme for the numerical integration in time. In order to increase the time step, residual averaging scheme of Jameson is adopted. Pressure field is obtained by solving the pressure-Poisson equation with the appropriate Neumann boundary condition. For the turbulence closure, 0-equation turbulence model of Baldwin-Lomax is used. Numerical computation is carried out for the Reynolds number of 4.5 million. Comparisons of the computed results with the available experimental data show good agreements for the velocity and pressure distributions.

• Journal of the Society of Naval Architects of Korea
• /
• v.33 no.1
• /
• pp.1-8
• /
• 1996

A potential-based panel method has been developed for numerical computation of wave resistance on a hybrid hydrofoil. Hybrid hydrofoil is composed of a main body, two struts and two hydrofoils. The main body, which is assumed to be an axisymmetric body for the present analysis, is normally used to support displacement of a body with its buoyancy. Normal dipoles and the sources are distributed on the body(main body, struts, hydrofoils) and the sources are distributed on the free surface. Linearized free surface and the radiation conditions are satisfied using the fourth order finite difference operator and the semi-linear pressure Kutta condition is used for the numerical computation of the hydrofoils. Poisson type free surface condition has been used for the numerical computation and hyperboloidal panel method has been used for better numerical accuracy. To verify this numeric method, model tests are performed in circulation water channel. From the comparison of experimental results with numeric ones, the present method can be used as a useful tool for the design of high speed vessels.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.50 no.6
• /
• pp.367-376
• /
• 2022

In this study, a steady and unsteady aerodynamic analysis program using a 3-dimensional subsonic unstructured panel method is developed and verified. Surfaces of bodies are modeled with the source and doublet distributions on triangular or quadrilateral panels. Geometry modeling of complex geometries and multi-body, therefore, can be easily accomplished. The Kelvin theory and the unsteady Kutta condition allow the doublet strength of the wake panels determined for unsteady flows. Various steady and unsteady flows in two and three dimensions are computed and compared with the analytical and the published computational results.

• Applied Chemistry for Engineering
• /
• v.9 no.7
• /
• pp.1079-1084
• /
• 1998

Reaction kinetics between 4,4'-dihexyl methane diisocyanate($H_{12}MDI$) and n-hexanol in toluene with dibutyltin dilaurate(DBTDL) as catalyst was studied by experimental measurements and mathematical modeling. Experiments were carried out at various temperatures, catalyst concentrations and [NCO]/[OH] ratios, and the reaction kinetics were described by two second-order reactions, the one between NCO and OH leading to urethane and the other between urethane and NCO leading to allophanate. The rate constants were estimated by the Runge-Kutta 4th-order method. Experiments and mathematical simulations showed a good agreement for various experimental conditions. The [allophanate]/[urethane] ratios at 90% conversion of initial NCO were estimated to be over 20% for most conditions employed in the present study, indicating that allophanate formation might significantly affect the properties of urethane polymers.