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

A two-dimensional Navier-Stokes solver using hybrid meshes is parallelized with a domain decompostion method. The focus of this paper is placed on minimizing the amount of effort in parallelizing the serial version of the solver, and this is achieved by adding an additional layer of cells to each decomposed domain. Most subroutines of the serial solver are used without modification, and the information exchange between neighboring domains is achieved using MPI(Message Passing Interface) library. Load balancing among the processors and scheduling of the message passing are implemented to reduce the overhead of parallelization, and the speed-up achieved by parallelization is measured on the transonic invisicd and turbulent flow problems. The parallelization efficiencies of the explicit Runge-Kutta scheme and the implicit point-SGS scheme are compared and the effects of various factors on the results are also studied.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.392-394
• /
• 2011

This paper is concerned about the hydrodynamic coefficients of hydrofoil. We discretized the incompressible Navier-Stokes equation with second order Runge-kutta for the time in the second order compact scheme for the spatial. The three-dimensional CFD code based on hybrid mesh on the finite volume method is used to simulate flow around NACA series foils. Lift and drag coefficient is calculated for several NACA series foils using different mesh types. Our aim is to obtain the lift and drag coefficient to evaluate the robustness of the solver and to shaw the advantage of using trim tab at the trailing edge. It concludes with a discussion of results and recommendations for future work.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2003.04a
• /
• pp.215-219
• /
• 2003

Buckling and post-buckling Analysis of Ludwick type and modified Ludwick type elastic materials was carried out. Because the constitutive equation, or stress-strain relationship is different from that of linear elastic one, a new governing equation was derived and solved by $4^{th}$ order Runge-Kutta method. Considered as a special case of combined loading, the buckling under both point and distributed load was selected and researched. The final solution takes distinguished behavior whether the constitutive relation is chosen to be modified or non-modified Ludwick type as well as linear or non-linear. We also derived strain energy function for non-linear elastic constitutive relationship. By doing so, we calculated the criterion function which estimates the stability of the equilibrium solutions and determines critical buckling load for non-linear cases. We applied this theory to the constitutive relationship of fabric, which also is the non-linear equation between the applied moment and curvature. This results has both technical and mathematical significance.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.2
• /
• pp.336-347
• /
• 1992

The three dimensional inviscid transonic cascade flow was investigated numerically, incorporation a four stage Runge-Kutta integration method proposed by Jameson. Time marching to the steady state was accelerated by using optimum time step and enthalpy damping. In describing the boundary conditions at inlet and outlet, Riemann invariants are considered. By adding a second and a fourth order artificial viscocities, the numerical instability due to the propagation of undamped disturbance or the rapid change of state near the shock has been prevented. The numerical results for are bump cascade, cambered two dimensional turbine cascade and three dimensional stator cascade agreed reasonably well with previous results. It has been known that the accuracy of the solution depended a lot on the modeling of the leading or trailing edge.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.813-818
• /
• 2000

The objective of this study is to elucidate the characteristic performance of Roots-type vacuum pump with Cassini oval lobe shape. The modelled lobe shape of Roots-type vacuum pump is two-lobe spur gear. The numerical analyses are performed fer leak flows, using 4th-order Runge-Kutta method and are compared with experimental results. Results show that for the case of involute lobe shape the total amount of the leak flow conductance is greater than that of cycloid and Cassini oval lobe shapes.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.1 no.3
• /
• pp.22-33
• /
• 1993

A theoretical analysis of predicting the detailed motion of a piston-crank mechanism within piston-guide clearance is presented, and the analysis is applied to the piston motion in a gasoline engine. A piston movement program is developed to calculate the piston attitude relative to the bore, the piston to bore impact velocity and kinetic energy loss and the net transverse force acting on the piston. This paper presents the formulation of a set of differential equations governing the transverse and rotational motion of a piston. These equations of motion were solved by well established Runge-Kutta method. As a result of this study, it is possible to predict the effects of piston geometry and piston pin offset on a piston motion and kinetic energy loss.

• Journal of Korean Society for Atmospheric Environment
• /
• v.23 no.6
• /
• pp.708-717
• /
• 2007

In the present study, one of the widely applied equations for gas-solid cyclones, Leith and Licht model, was evaluated based on the 3-D CFD technique. The initial and boundary values of radial position and tangential velocity obtain-ed from the CFD simulation enabled complete calculation of the nonlinear second differential equation. This approach showed about 30% errors between calculations with and without the second order differential term. The calculation by using the simple first order equation presented shorter times to migrate up to the inner wall of the cyclone than by the second order, which theoretically implies higher separation efficiency. Further comparison is now under evaluation in terms of the detailed grade efficiency.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2003.11a
• /
• pp.906-910
• /
• 2003

This paper numerically analyzes the dynamic characteristics of a spindle system supported by two identical journal bearings considering bearing span that has dynamic load due to its mass unbalance. The Reynolds equation is transformed to solve a herringbone grooved journal bearing. The Reynolds equations are solved using FEM in order to calculate the pressure distribution in a fluid film. Reaction forces and friction torque are obtained by integrating the pressure and shear stress along the fluid film, respectively. Dynamic behaviors, such as whirl radius or angular displacement of a rotor, are determined by solving its nonlinear equations of motion with the Runge-Kutta method. This research shows that the same bearing spans of upper and lower journal bearings produce the minimum runout and friction torque of a spindle system.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.21 no.2
• /
• pp.137-145
• /
• 2011

This study is proposed to develop a numerical interaction model of the magnetically levitated(maglev) train and guideway. For this purpose, equation of motion for 6-DOF vehicle model, EMS, guideway and guideway irregularity are derived as the state-space equation. In order to solve the state space equations, the present work was performed via matlab simulation using Runge-Kutta method. Through the simulation, the effect of dynamic response of maglev system to different vehicle speeds, guideway rigidity(EI) and masses is investigated.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2002.05a
• /
• pp.417-422
• /
• 2002

The main purpose of this paper is to determine the dynamic optimal shapes of tapered column with constant volume. The linear, parabolic and sinusoidal tapers with the regular polygon cross-section are considered, whose material volume and span length are always held constant. The ordinary differential equation including the effect of axial load is applied to calculate the natural frequencies. The Runge-Kutta method and Regula-Falsi methods are used to integrate the differential equation and compute the frequencies, respectively. Then the dynamic optimal shape whose lowest natural frequency is highest is determined by reading the critical value of the frequency versus section ratio curve plotted by the frequency data. In the numerical examples, the tapered columns are analysed and the numerical result of this study are shown in table and figures.