• 대한기계학회논문집B
• /
• 제21권2호
• /
• pp.202-212
• /
• 1997

In this study, performance of the multi-stage explicit methods for numerical computation of the unsteady Navier-Stokes equations is investigated. Three methods under consideration are 1 st-, 2 nd-, and 4 th-order Runge-Kutta (R-K) methods. Compared in this estimation is stability, accuracy, and CPU time of each method. The computational codes developed are applied to the two-dimensional flow in a square cavity driven by an oscillating lid. It turned out that at Reynolds number 400, the 1 st-order R-K method is the best, while at 3200 the 2 nd-order R-K is recommended. At higher Reynolds numbers, it is conjectured that the 4 th-order R-K method will be the best algorithm among three due to its highest stability.

• 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.

• 한국윤활학회:학술대회논문집
• /
• 한국윤활학회 1997년도 제25회 춘계학술대회
• /
• pp.204-209
• /
• 1997

This paper reports on the theoretical study of the vane motions in a positive displacement vane pump. Vane detachment cause the pressure fluctuation, noise, wear in cam ring, and decrease the volumetric efficiency. Dynamic equation of vane motion and flow continuity equation have been modeled and solved simultaneously using 4th order Runge-Kutta method. As results of analysis, vane detachment occurs due to pressure overshoot by excess compression in the pumping chamber. Amount of vane detachment has been reduced by decreasing the pressure overshoot.

• 대한전기학회논문지
• /
• 제39권5호
• /
• pp.516-518
• /
• 1990

We have derived the equation which involves the variation of electron energy with time in a lowly ionized plasma when a low alternating electric field is applied. We consider only elastic collisions between electrons and neutral atoms. This equation is solved using the 4th-order Runge-Kutta method, and applied to argon gas discharge which is driven by source frequency of 100, 1K, 10K, 100K, and 1M (Hz). The results show that the variation of electron energy becomes flat with higher frequencies.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2006년도 PARALLEL CFD 2006
• /
• pp.311-314
• /
• 2006

The supersonic flow around tandem cavities was investigated by three- dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes(RANS) equation with the $\kappa-\omega$ thrbulence 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 using 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 ratio 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 fire cavity flow 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.

• 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.

• 대한토목학회논문집
• /
• 제31권2A호
• /
• pp.43-49
• /
• 2011

본 연구에서는 선교나 부선등에 실린 장비들과 같이 선박의 회전중심으로부터 이격된 선박구조와 교량의 상부구조가 충돌하는 경우의 충돌거동을 해석할 수 있도록 선박의 회전을 고려한 간이충돌모델을 개발하였다. 충돌모델은 선박과 교량의 질량과 강도 및 비선형성을 나타내도록 하였고, 모델의 비선형 운동방정식을 4차 Runge-Kutta방법을 이용하여 수치적으로 해석하였다. 현재 국내에서 운용되고 있는 부선의 제원을 선폭을 기준으로 표준화하고 작성된 충돌모델을 이용하여 각종 충돌조건이 교량의 충돌하중이 미치는 영향에 대한 매개변수연구를 수행하였다.

• 한국수자원학회논문집
• /
• 제38권6호
• /
• pp.429-438
• /
• 2005

본 연구에서는 Boussinesq 방정식을 이용하여, 불규칙 파랑의 직접적인 해석이 가능한 한 쌍의 상미분방정식을 유도하였다. 입사파랑은 TMA(TEXEL storm, MARSEN, ARSLOE) 천해 스펙트럼을 이용하여 재현하였으며, 지배방정식은 4차 Runge-Kutta 법을 이용하여 적분하였다. 새로 유도된 파랑 방정식을 이용하여, 일정 수심을 진행하는 파랑의 비선형 에너지 교환효과를 계산하였다. 또한, 일정 경사면의 정현파형 지형을 통과하는 불규칙파랑의 특성에 관해 수치적으로 검토하였다. 비선형성이 불규칙파랑의 통과와 반사에 큰 영향을 주었다.

• 한국해안·해양공학회논문집
• /
• 제36권2호
• /
• pp.70-79
• /
• 2024

본 연구에서는 3차원 공간에서 바닥의 움직임에 따른 선형파의 생성을 모의할 수 있는 스펙트럼 법을 소개한다. 지배방정식은 선형의 동역학적 및 운동학적 자유수면 경계조건이며, 두 식은 Fourier 공간에서 해석된다. 해석된 속도포텐셜 및 자유수면변위는 연속방정식과 운동학적 바닥경계조건을 항상 만족해야 한다. 수치해석에서 시간 적분은 4차 Runge-Kutta 법을 이용하여 해석하였다. Fourier 공간에서 해석한 결과는 Fourier 역변환을 통해 실제 공간에서의 속도포텐셜과 자유수면변위로 표현된다. 본 수치모델을 이용하여 다양한 형상의 바닥이 규칙적으로 움직이는 경우 생성되는 규칙파에 대해 모의하였다. 또한 바닥의 움직임을 이용하여 비스듬히 전파하는 규칙파의 생성도 모의하였다. 수치모델의 결과는 해석해와 비교하였으며, 거의 일치하는 결과를 보였다.