• Journal of Mechanical Science and Technology
• /
• 제20권8호
• /
• pp.1256-1265
• /
• 2006

The supersonic flows around tandem cavities were investigated by two-dimensional and three-dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes (RANS) equation with the k- ω turbulence model. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split with van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge- Kutta method. The aspect ratios of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two- dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the first cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.

• Tribology and Lubricants
• /
• 제33권6호
• /
• pp.309-314
• /
• 2017

This paper presents a theoretical investigation of the dynamic characteristics of externally pressurized air pad bearings with closed loop grooves. These grooves are made on the surface of bearings to reduce the number of supply holes so that manufacturing costs can be reduced. The semi-implicit method is applied to calculate the time varying pressure profile on the air bearing surface owing to the advantages of numerical stability and fast time tracing characteristics. The static pressure of the groove bearings is much higher than that without grooves, so the groove bearings can provide high load carrying capacity. The equation of motion considering vertical motion and tilting motion are also solved using the Runge-Kutta 4th order method. By combining the semi-implicit method and the Runge-Kutta method, fast calculations of the dynamic behavior of the air bearing can be achieved. The variations of bearing reaction force, air film reaction moment, height, and tilting angle are investigated for the step force input, which is 20% higher than the bearing reaction, when the nominal clearance is 6 mm. The effect of the groove width and the groove depth are investigated by calculating the dynamic behavior. The possibility of the air hammering with the depth of the groove is found and discussed.

• 대한설비공학회:학술대회논문집
• /
• 대한설비공학회 2006년도 하계학술발표대회 논문집
• /
• pp.116-121
• /
• 2006

Numerical simulations for dendritic growth of crystals are conducted in this study by the level set method. The effect of order of difference is tested for reinitialization error in simple problems and authors founded in case of 1st order of difference that very fine grids have to be used to minimize the error and higher order of difference is desirable to minimize the reinitialization error The 2nd and 4th order Runge-Kutta scheme in time and 3rd and 5th order of WENO schemes with Godunov scheme are applied for space discretization. Numerical results are compared with the analytical theory, phase-field method and other researcher's level set method.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2002년도 학술대회지
• /
• pp.783-786
• /
• 2002

The unsteady turbulent flow characteristics of NACA4406 airfoil at close proximity to the pound are numerically investigated, especially focused on the noise generation mechanism near the blunt trailing edge. The unsteady two-dimensional compressible Wavier-Stokes equations with a Spalart-Allmaras turbulence closure model are solved by the 6th-order compact scheme and the 4th-order Runge-Kutta scheme. The computation shows a noise generation by a feedback mechanism at the blunt tailing edge, where the acoustic-fluidic coupling occurs between the wall-reflected sound waves and the periodically disturbed turbulent shear layer.

• 한국해양공학회지
• /
• 제20권4호
• /
• pp.31-36
• /
• 2006

Numerical schemes for spacing and time are tested to accurately simulate the wave propagation. The tested numerical schemesinclude 2nd-order central differencing, l-order upwind scheme, 2nd-order Leith scheme, 3rd-order MUSCLE, QUICK and QUICKEST schemes in spacing and the Euler and 4th-order Runge-Kutta(R-K) schemes in time. It is seen that more accurate results are expected when the higher-order schemes, especially the schemes combined with a TVD control limiter, are used for solving the wave equation. The 3rd-order upwind scheme with limiter and the 4th-order R-K scheme in tim￡ are finally applied to the wave-making simulation in a digital wave tank.

• Wind and Structures
• /
• 제16권4호
• /
• pp.323-340
• /
• 2013

An active mass damper system for flutter control of bridges is presented. Flutter stability of bridge structures is improved with the help of eccentric rotational actuators (ERA). By using a bridge girder model that moves in two degrees of freedom and is subjected to wind, the equations of motion of the controlled structure equipped with ERA are established. In order to take structural nonlinearities into consideration, flutter analysis is carried out by numerical simulation scheme based on a 4th-order Runge-Kutta algorithm. An example demonstrates the performance and efficiency of the proposed device. In comparison with known active mass dampers for flutter control, the movable eccentric mass damper and the rotational mass damper, the power demand is significantly reduced. This is of advantage for an implementation of the proposed device in real bridge girders. A preliminary design of a realization of ERA in a bridge girder is presented.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 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.

• Tribology and Lubricants
• /
• 제13권4호
• /
• pp.53-59
• /
• 1997

The camshaft locus at camshaft beating in a direct acting OHC valve train system has been investigated using the transient method. Forces applied to the camsfiaft are composed of two components, one is the transfer force between the cam and the tappet, the other is the frictional force. These forces have been calculated using the lumped mass model and the elastohydrodynamic lubrication theory. The alternating direction implicit method has been used for the numerical analysis of Reynolds equation, and 4th order Runge-Kutta method has been used for the transient journal locus analysis. The effects of various load conditions are presented in the form of journal locus. As a result of the analysis, it has been found that camshaft bearings were mainly in the hydrodynamic lubrication condition.

• 항공우주기술
• /
• 제6권1호
• /
• pp.19-28
• /
• 2007

본 연구에서는 올레오-뉴메틱(Oleo-pneumatic) 완충장치와 차등 제동 조향(Differential Braking Steering)을 위한 캐스터링 휠 포크(Castering Wheel Fork)가 적용된 소형항공기용 전방착륙장치의 충격하중 해석을 위한 비선형 2 자유도 수학적 모델을 생성하고, 수치해석 기법을 이용하여 시간에 따른 충격응답을 구하였다. 생성된 전방착륙장치의 수학적 모델은 타이어의 충격 시 동적 특성과 스트루트 기울기, 캐스터에 의한 모멘트, 타이어 접촉과 조향 등에 의해 발생하는 측면하중에 의한 올레오 스트루트 내부의 마찰 특성 등의 비선형 요소들을 포함하고 있다. 또 생성된 운동방정식의 수치해석은 4차 Runge-Kutta 방식을 이용하였으며, 운동방정식 및 수치해의 정확도를 평가하기 위한 타 항공기에 기 적용된 전방 착륙장치의 낙하시험 결과와의 비교 검증과정을 수행하였다.

• 한국추진공학회지
• /
• 제11권2호
• /
• pp.79-85
• /
• 2007

Polyethylene-GOX (PE-GOX) 조합의 추진제를 사용하는 하이브리드 모터의 성능 예측 모델을 제시하였다. 성능에 영향을 미치는 인자들에 대해 자세한 고찰이 이루어 졌으며, 연료의 후퇴율과 산화제 조절 밸브의 개폐 응답 특성을 고려한 비정상 지배방정식을 구성하고 연소 중 O/F(산화제/연료)비에 따라 변화하는 열역학 데이터를 고려한 4차 Runge-Kutta 수치기법으로 해를 구하여 실험결과와 비교하여 잘 일치함을 확인하였다. 또한 연소 중 그레인의 후퇴율과 O/F비등의 성능 변수들의 영향을 분석하였다.