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

Three-Dimensional Euler equations are solved numerically for the analysis of contraction flows in wind or water tunnels. A second-order finite difference method is used for the spatial discretization on the nonstaggered grid system and the 4-stage Runge-Kutta scheme for the numerical integration in time. In order to speed up the convergence, the local time stepping and the implicit residual-averaging schemes are introduced. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. For the evaluation of the present Euler solver, numerical computations are carried out for the various contraction geometries, one of which was adopted in the Large Cavitation Channel for the U.S. Navy. The comparison of the computational results with the available experimental data shows good agreements.

• 한국전산유체공학회지
• /
• 제2권1호
• /
• pp.8-12
• /
• 1997

The three-dimensional Euler equations are solved numerically for the analysis of contraction flows in wind or water tunnels. A second-order finite difference method is used for the spatial discretization on the nonstaggered grid system and the 4-stage Runge-Kutta scheme for the numerical integration in time. In order to speed up the convergence, the local time stepping and the implicit residual-averaging schemes are introduced. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. For the evaluation of the present Euler solver, numerical computations are carried out for three contraction geometries, one of which was adopted in the Large Cavitation Channel for the U.S. Navy. The comparison of the computational results with the available experimental data shows good agreement.

• 항공우주시스템공학회지
• /
• 제10권1호
• /
• pp.1-7
• /
• 2016

During flight of the aircraft, the vortex merging phenomenon appears under the certain condition between co-rotating vortices which were generated at the wing tip and lifting-surface. And then these merged vortices at both sides show counter-rotating pattern to affect on the downstream of the aircraft. In this paper, the numerical simulations are conducted assuming this phenomenon in two-dimensional co-rotating or counter-rotating vortices pairs. Two-dimensional incompressible Navier-Stokes equations were converted into Vorticity-Streamfunction form and the Galerkin spectral method was adopted. The third order Runge-Kutta method was used for time integration. The effects on the vortex merger and degree of vortex merger were investigated according to time, Reynolds number, and changes in the distance between two vortices.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2004년도 춘계학술대회논문집
• /
• pp.894-899
• /
• 2004

In this paper, the track pull-in behavior analyses in an optical disk drive (ODD) using plane phase is treated and the parameters affecting it are discussed. Track pull-in is the track capture procedure to do track following control and it is key factor to increase data transfer rate. Simulation method, Runge-Kutta method to solve nonlinear equation, is used to evaluate the track pull-in conditions, and the real servo loop parameters are applied in this process to get the more real condition. Finally, the comments for the acquired results are discussed briefly.

• 대한기계학회논문집
• /
• 제10권5호
• /
• pp.751-757
• /
• 1986

본 연구에서는 준 3차원 유동에서 문제가 되는 B-B 유동만을 대상으로 하고 여기에 두개의 유동함수를 이용해서 수치계산한다. 이 방법의 정당성을 입증하기 위해 계산결과를 이론해와 비교하고 또 실제유동과의 관계를 보기 위해 실험결과와도 비교하여 또 실제유동과의 관계를 보기 위해 실험결과와도 비교하여 Wu에 의한 정상 비점성 유동의 준3차원적 수치계산방법의 유효성을 제시한다.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2006년도 춘계학술대회논문집
• /
• pp.396-401
• /
• 2006

Dynamic behaviors of the impact damper are studied experimentally and numerically. In order to investigate wide range of excitation frequencies and amplitudes, a simple but high amplifying and bias-free experimental setup is designed. Experiments focused on the harsh operation condition demonstrate Accelerated mass loading which not only deteriorates the performance of the impact damper but also involves the structural resonance which should be avoided for the stability of the system. In the previous studies, instability or deterioration of the performance was reported for the off resonance frequency region. But this paper shows that the performance deterioration and structural resonances can be predicted. Using finite element modeling and analysis, accurate system parameters were derived and used for the numerical modeling employing the conservation of the momentum. Numerical study of the transient responses using 4th-order Runge-Kutta method demonstrates general performance of the system, and shows that accelerated mass loading phenomenon is deeply related with the vibration amplitudes and the mass of the auxiliary system.

• 대한전기학회논문지:전기기기및에너지변환시스템부문B
• /
• 제51권10호
• /
• pp.560-567
• /
• 2002

This paper presents the calculation of opening characteristics for puffer GCB with the equations of the flow field and the motion of the driving mechanism. To obtain the stroke curve, the motion equation is solved simultaneously with the Euler equations. For a given Piston location, the flow field is solved. The pressure inside the Puffer chamber is then used to calculate the moving velocity and the new position of the piston. The FVFLIC method is employed to solve the axisymmetric Euler equations and the motion equation is solved by the Runge-Kutta method. The method is applied to the puffer GCB model and the stroke curve and the pressure rise in puffer chamber under no load condition are compared with the measured ones.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2002년도 봄 학술발표회 논문집
• /
• pp.115-122
• /
• 2002

The main purpose of this paper is to determine the static optimal shapes of tapered beams 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 Runge-Kutta method is used to integrate the differential equation and also Shooting method is used to calculate the unknown boundary condition. Then the static optimal shapes are determined by reading the minimum values of the deflection versus section ratio curves plotted by the deflection data. In numerical examples, the various tapered beams are analyzed and those numerical results of this study are shown in figures.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2008년 영문 학술대회
• /
• pp.597-601
• /
• 2008

Three-dimensional structures of detonation wave propagating in circular tube were investigated. Inviscid fluid dynamics equations coupled with a conservation equation of reaction progress variable were analyzed by a MUSCL-type TVD scheme and four stage Runge-Kutta time integration. Variable-$\gamma$ formulation was used to account for the variable properties between unburned and burned states and the chemical reaction was modeled by using a simplified one-step irreversible kinetics model. The computational code was parallelized based on domain decomposition technique using MPI-II message passing library. The computations were carried out using a home made Windows based PC cluster having 160 AMD AthloxXP and Athlon64 processor. The computational domain consisted of through a roundshaped tube with wall conditions. As an initial condition, analytical ZND solution was distributed over the computational domain with disturbances. The disturbances has circumferential large gradient. The unsteady computational results in three-dimension show the detailed mechanisms of multi-cell mode of detonation wave instabilities resulting diamond shape in smoked-foil record.

• Journal of Ship and Ocean Technology
• /
• 제5권1호
• /
• pp.38-54
• /
• 2001

A finite difference method for calculating turbulent flow and surface wave fields around a ship model is evaluated through the comparison with the experimental data of a Series 60 $C_B$=0.6 ship model. The method solves the Reynolds-averaged Navior-Stokes Equations using the non-staggered grid system, the four-stage Runge-Kutta scheme for the temporal integration of governing equations and the Bladwin-Lomax model for the turbulence closure. The free surface waves are captured by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and free-surface conforming grids are generated at each time step so that one of the grid surfaces coincides always with the free surface. The computational results show an overall close agreement with the experimental data and verify that the present method can simulate well the turbulent boundary layers and wakes as well as the free-surface waves.