• Tribology and Lubricants
• /
• 제23권5호
• /
• pp.219-227
• /
• 2007

The dynamic performance of air foil bearings relies on a coupling between a thin air film and an elastic foil structure. A number of successful analytical techniques to predict dynamic performance have been developed. However, the evaluation of its dynamic characteristic is still not enough because of the mechanical complexity of the foil structure and strong nonlinear behavior of friction force. This work presents a nonlinear transient analysis method to predict dynamic performance of foil bearings. In this method, time dependent Reynolds equation is used to calculate pressure distribution and a finite element method is used to model the bump foil structure. The analysis is treated with a direct implicit integration technique that can handle nonlinear problems and the stick-slip algorithm is used to consider friction force. Using this method the response to the mass unbalance excitation is investigated for various design parameters and operating conditions. The results of analysis show that foil bearing is very effective on the restriction of vibration at the resonance frequency compared to the rigid surface bearings and the effectiveness depends on the operating conditions, static load and a amount of mass unbalance. In addition, there exist optimum values of friction coefficient, bump foil stiffness and number of circumferential slit with regards to minimizing dynamic response at the resonance frequency. These optimum values are system dependent.

• 한국자동차공학회논문집
• /
• 제7권8호
• /
• pp.91-98
• /
• 1999

This paper describes the air flow characteristics inside the throttle valve. Tow-dimensional steady incompressible Navier-Strokes equation are solved numerically with embedding the conceopt of the artificial compressibility and adopting the Baldwin-Lomax turbulence model. With varying the valve opening angles(the Reynolds number )such as 15$^{\circ}$(5000) , 45$^{\circ}$(3000) , 75$^{\circ}$(7000) and 90$^{\circ}$(10000), respectively. tow cases, with a valve shaft and without one, are analysed. The pressure loss between the entrance and exit is severe at 15$^{\circ}$, 100 times as larger as that of 90$^{\circ}$ case, which also depends much on the existece of the valve shaft. The counter rotating vortices are formed over the valve plate with the shaft at only 75$^{\circ}$. They are smally and very large scale in front and back of the valve shaft , respectively. The velocity profiles of 15$^{\circ}$ and 90$^{\circ}$ at the exit are almost symmetric to the horizontal center line, however, the symmetricity is no longer maintained at 45$^{\circ}$ and 75$^{\circ}$ , and in addition, the flow at 75$^{\circ}$ is enforced a lot below center line. The pressure distribution on the walls is largely changed near the valve shaft, and its magnitude becomes great as the valve angle decreases.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제11권2호
• /
• pp.883-898
• /
• 2019

This paper employs computational tools to predict power increase (or speed loss) and propulsion performances in waves of KVLCC2. Two-phase unsteady Reynolds averaged Navier-Stokes equations have been solved using finite volume method; and a realizable k-ε model has been applied for the turbulent closure. The free-surface is obtained by solving a VOF equation. Sliding mesh method is applied to simulate the flow around an operating propeller. Towing and self-propulsion computations in calm water are carried out to obtain the towing force, propeller rotating speed, thrust and torque at the self-propulsion point. Towing computations in waves are performed to obtain the added resistance. The regular short head waves of λ/LPP = 0.6 with 4 wave steepness of H/λ = 0.007, 0.017, 0.023 and 0.033 are taken into account. Four methods to predict speed-power relationship in waves are discussed; Taylor expansion, direct powering, load variation, resistance and thrust identity methods. In the load variation method, the revised ITTC-78 method based on the 'thrust identity' is utilized to predict propulsive performances in full scale. The propulsion performances in waves including propeller rotating speed, thrust, torque, thrust deduction and wake fraction, propeller advance coefficient, hull, propeller open water, relative rotative and propulsive efficiencies, and delivered power are investigated.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제11권2호
• /
• pp.736-749
• /
• 2019

This paper employs computational tools to investigate the diffraction effects in regular head waves on the added resistance and wake on the propeller plane. The objective ships are a 66,000 DWT bulk carrier and a 3,600 TEU container ship. Fixed and free to heave and pitch conditions at design speed have been taken into account. Two-phase unsteady Reynolds averaged Navier-Stokes equations have been solved using the finite volume method; and a realizable k-ε model has been applied for the turbulent closure. The free surface is obtained by solving a VOF equation. The computations are carried out at the same scale of the model tests. Grid and numerical wave damping zones are applied to remove unwanted wave reflection at the boundaries. The computational results are analyzed using the Fourier series. The added resistances in waves at the free condition are higher than those at the fixed condition, which are nearly constant for all wavelengths. The wake velocity in waves is higher than that in calm water, and is accelerated where the wave crest locates on the propeller plane. When the vertical motion at the stern goes upward, the wake velocity also accelerated.

• 대한조선학회논문집
• /
• 제48권2호
• /
• pp.147-157
• /
• 2011

A RANS(Reynolds averaged Navier-Stokes) based numerical method is developed for the evaluation of ship resistance and self-propulsion performances. In the usability aspect of CFD for the hull form design, the field grid around practical hull forms is generated by solving a grid Poisson equation based on the hull surface grid generated from station offsets and centerline profile. A body force technique is introduced to model the effects of the propeller in which the propeller loads are obtained from potential flow analysis using an unsteady lifting surface method. The free surface is captured by using a two-phase level-set method and the realizable $k-{\varepsilon}$ model is used for turbulence closure. The hull attitude in vertical plane, i.e., trim and sinkage, is calculated by using a quasi-steady method and then considered in the computation by translating and rotating the grid system according to the values. For the validation of the proposed method, the numerical results of resistance tests for KCS, KLNG, and KVLCC1 and of self-propulsion test for KCS are compared with experimental data.

• 한국유체기계학회 논문집
• /
• 제20권1호
• /
• pp.15-20
• /
• 2017

The small turbocharger for the automotive application is designed to operate up to 200,000 rpm to increase system efficiency. Because of high rotation speed of turbocharger, floating ring bearing are widely adopted due to its low friction loss and high rotordynamic stability. This paper presents a linear and nonlinear analysis model for a turbocharger rotor supported by a semi-floating ring bearing. The rotordynamic model for the turbocharger rotor was constructed based on the finite element method and fluid film forces were calculated based on the infinitely short bearing assumption. In linear analysis, we considered fluid film force as stiffness and damping element and in nonlinear analysis, the fluid film force was calculated by solving the time dependent Reynolds equation. We verified the developed theoretical model by comparing to modal test results of test rotors. The analysis results show that there are two unstable modes, which are conical and cylindrical modes. These unstable modes appear as sub-synchronous vibrations in nonlinear analysis. In nonlinear analysis, frequency jump phenomenon demonstrated when vibration mode is changed from conical mode to cylindrical one. This jump phenomenon was also demonstrated in the test. However, the natural frequency measured in the test differs from those obtained using nonlinear analysis.

• 대한기계학회논문집B
• /
• 제34권10호
• /
• pp.925-931
• /
• 2010

본 논문은 공동 주위 난류유동특성을 LES 기법으로 수치해석을 수행하여 알아보았다. 본 연구에 적용된 레이놀즈수는 공동 깊이만큼의 높이 h 에서의 유속을 기준으로 $1.0{\times}10^5$ 이며 3 차원 공동에서의 유동특성을 알아보았다. 적절한 비압축성 Filtered Navier-Stokes 방정식을 적용하기 위해, 계산격자를 공동 표면 근처에는 조밀하게 멀어질수록 성기게 생성하였으며, 이는 계산시간을 단축시키며 빠른 수렴을 도와준다. 또한, Boussinesq 가설을 subgrid-scale 난류모델에 적용하였고, Subgrid-scale 난류점성을 얻기 위해 smagorinsky-Lilly SGS 모델을 적용하였으며, 그 때의 CFL 수는 1.0 이다. 또한, 본 논문은 서로 다른 4 가지 형상의 공동의 및 입구조건의 변화에 따른 유동 특성도 함께 연구되었다.

• 한국추진공학회지
• /
• 제2권3호
• /
• pp.20-35
• /
• 1998

벽면을 통한 유체의 유동이 있는 난류운동장의 특성을 분석하였다. 벽면을 통한 유체의 유동은 유출과 유입의 경우를 포함하여 Re$_{w}$ 의 절대크기가 0~160으로 다양하며 입구에서의 $Re_{in}$도 3${\times}$$10^3$에서 8${\times}$$10^4$까지 넓은 범위를 대상으로 하였다. 벽면을 통한 유체의 유동은 벽면 경계층의 변화를 초래하고 이에 따라 난류 유동장의 급격한 변화를 일으켜 벽면을 포함한 전 유동장의 특성을 변화시킨다. 이러한 복잡한 유동장을 올바로 예측하기 위하여 난류운동에너지 소산율 방정식의 자 항에 대한 모델을 개선한 수정 $\kappa$-$\varepsilon$ 모델을 도입하였으며 기존의 난류모델과의 비교를 통하여 성능 검증을 시도하였다. 해석을 통한 분석결과로부터 수정 $\kappa$-$\varepsilon$ 모델은 벽면을 통한 유체유동이 있는 복잡한 유동장을 올바로 예측할 수 있음을 알 수 있었다. 따라서 수정 $\kappa$-$\varepsilon$ 모델을 이용하여 다양한 경우의 벽면 및 입구 조건을 갖는 난류유동장을 해석할 수 있을 것으로 판단된다.단된다.

• 설비공학논문집
• /
• 제10권3호
• /
• pp.324-333
• /
• 1998

Measurements of the local heat transfer coefficients on a hemispherically concave surface with a round oblique impinging jet were made. The liquid crystal transient method was used for these measurements. This method, which is a variation of the transient method, suddenly exposes a preheated wall to an impinging jet while video recording the response of liquid crystal for the surface temperature measurements. The Reynolds number used was 23,000 and the nozzle -to -jet distance was L/d=2, 4, 6, 8 and 10 and the jet angle was $\alpha$=0$^{\circ}$, 15$^{\circ}$, 30$^{\circ}$and 40$^{\circ}$. In the experiment, the maximum Nusselt number at all region occurred at L/d(equation omitted)6 and Nusselt number decreases as the inclined jet angle increases. For the normal jet the contours of constant Nusselt number are circular and as the jet is inclined closer and closer to the surface the contours become elliptical shape. The decreasing rate of the Nusselt number at X/d＞ 0(upstream) on a surface curvature are higher than those on a flate plate and the decreasing rate of the Nusselt number at X/d ＜0(downstream) on a surface curvature are lower than those on a flate plate. And also, the decreasing rate of local Nusselt number distribution at X/d ＜0(upstream) exhibit lower than with X/d ＜0(downstream) as jet angle increases. The second maximum Nusselt number occurred at long distance from stagnation point as jet angle increases.

• 한국태양에너지학회 논문집
• /
• 제35권2호
• /
• pp.1-10
• /
• 2015

Promising onshore wind farm sites in Gangwon province of Korea were investigated in this study. Gangwon province was divided into twenty five simulation regions and a commercial program based on Reynolds averaged Navier-Stokes equation was used to find out wind resource maps of the regions. The national wind atlas with a period 2007-2009 developed by Korea institute of energy research was used as climatologies. The wind resource maps were combined to construct a wind resource map of Gangwon province with a horizontal spatial resolution of 100m. In addition to the wind resource, national environmental zoning map, distance from substation, residence and automobile road, Beakdudaegan mountain range, terrain slope, airport and military reservation district were considered to find out promising wind farm sites. A commercial wind farm design program was used to find out developable wind farm capacities in promising wind farm site with and without excluding environmental protection regions. The total wind farm capacities with and without excluding the protection regions were estimated to be 46MW and 598MW, respectively, when a 2MW commercial wind turbine was employed.