• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.9
• /
• pp.1245-1252
• /
• 2001

The condensation pressure drop for pure refrigerants R-22, R-134a, and a binary refrigerant mixture R-410A flowing in a small diameter tube was investigated. The test section is a counterflow heat exchanger with refrigerant flowing in the inner tube and coolant flowing in the annulus. The test section consists of 1220 [mm] length with horizontal copper tube of 3.38 [mm] outer diameter and 1.77 [mm] inner diameter. The refrigerant mass fluxes ranged from 450 to 1050 [kg/(㎡$.$s)] and the average inlet and outlet qualities were 0.05 and 0.95, respectively. The main experimental results were summarized as follows : In the case of single-phase flow, the pressure drop of R-134a is much higher than that of R-22 and R-410A for the same Reynolds number. The friction factors for small diameter tubes are higher than those predicted by Blasius equation. In the case of two-phase flow, the pressure drop increases with increasing mass flux and decreasing quality. The pressure drop of R-134a is much higher than that of R-22 and R-410A for the same mass flux. Most of correlations proposed in the large diameter tube showed enormous deviations with experimental data. However, the correlation predicted by Honda et al showed relatively good agreement with experimental data for R-134.

• Tribology and Lubricants
• /
• v.23 no.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.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.7 no.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
• /
• v.11 no.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
• /
• v.11 no.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.

• Journal of the Society of Naval Architects of Korea
• /
• v.48 no.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.

• The KSFM Journal of Fluid Machinery
• /
• v.20 no.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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.34 no.10
• /
• pp.925-931
• /
• 2010

The paper describes a Large Eddy Simulation (LES) study of turbulent flow around a cavity. A series of three-dimensional cavities placed in a turbulent boundary layer are simulated at a Reynolds number of $1.0{\times}10^5$ by considering U and h, which represent the velocity at the top and the depth of the cavity, respectively. In order to obtain the appropriate solution for the filtered Navier-Stokes equation for incompressible flow, the computational mesh forms dense close to the wall of the cavity but relatively coarse away from the wall; this helps reduce computation cost and ensure rapid convergence. The Boussinesq hypothesis is employed in the subgrid-scale turbulence model. In order to determine the subgrid-scale turbulent viscosity, the Smagorinsky-Lilly SGS model is applied and the CFL number for time marching is set as 1.0. The results show the flow variations inside cavities of different sizes and shapes.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.2 no.3
• /
• pp.20-35
• /
• 1998

Characteristics of turbulent flow with wall transpiration is analyzed. The wall transpiration includes both of suction and injection and extends their range to 0~160 of absolute magnitude of Re$_{w}$ . Reynolds number based on inlet velocity also covers wide range of 3${\times}$$10^3$~8${\times}$$10^4$. The turbulent flow with wall transpiration induces change of wall boundary layer and rapid change of turbulent field. This, in turn, leads the change of whole flow field. For predicting this complicated flow field properly, newly modified $\kappa$-$\varepsilon$ model is utilized, which is formed by modifying dissipation rate equation. The modified $\kappa$-$\varepsilon$ model of Chien is also adopted for the comparison of model performance. Analysis shows the newly modified $\kappa$-$\varepsilon$ model is successfully able to reflect the characteristics of turbulent flow field with wall transpiration.ion.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.10 no.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.