• Journal of Mechanical Science and Technology
• /
• v.17 no.7
• /
• pp.1054-1062
• /
• 2003

The numerical simulations of flow and pollutant particle dispersion are described for two-dimensional bell shaped hills with various aspect ratios. The Reynolds-averaged incompressible Navier-Stokes equations with low Reynolds number $\kappa$-$\varepsilon$ turbulent model are used to simulate the flowfield. The gradient diffusion equation is used to solve the pollutant dispersion field. The code was validated by comparison of velocity, turbulent kinetic energy, Reynolds shear stress, speed-up ratio, and ground level concentration with experimental and numerical data. Good agreement has been achieved and it has been found that the pollutant dispersion pattern and ground level concentration have been strongly influenced by the hill shape and aspect ratio, as well as the location and height of the source.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.10
• /
• pp.1940-1954
• /
• 1992

Fine grid computations were attempted to analyze the turbulent flows in the near wall low Reynolds number region and the numerical analyses were incorporated by a finite-volume discretization with full find grid system and low Reynolds number k-.epsilon. model was employed in this region. For the improvement of low Reynolds number k-.epsilon. model, modification coefficient of eddy viscosity $f_{\mu}$ was derived as a function of turbulent Reynolds number $R_{+}$ and nondimensional length $y^{+}$ from the concept of two length scales of dissipation rate of turbulent kinetic energy. The modification coefficient $f_{\epsilon}$ in .epsilon. transport equation was also derived theoretically. In the turbulent kinetic energy equation, pressure diffusion term was added in order to consider low Reynolds number region effect. The main characteristics of this low Reynolds number k-.epsilon. model were founded as : (1) In high Reynolds number region, the present model has limiting behavior which approaches to the high Reynolds number model. (2) Present low Reynolds number k-.epsilon. model dose not need additional empirical constants for the transport equations of turbulent kinetic energy and dissipation of turbulent kinetic energy in order to consider wall effect. Present low Reynolds number turbulence model was tested in the pipe flow and obtained improved results in velocity profiles and Reynolds stress distributions compared with those from other k-.epsilon. models.s.s.

• Tribology and Lubricants
• /
• v.16 no.4
• /
• pp.295-301
• /
• 2000

The dynamic simulation of slider in hard disk drive is performed using Factored Implicit Finite Difference method. The modified Reynolds equation with Fukui and Kaneko model is employed as a governing equation. Equations of motion for the slider of three degrees of freedom are solved simultaneously with the modified Reynolds equation. The transient responses of the slider for disk step bumps and slider impulse forces are shown for various cases and are compared for the iteration algorithm and new algorithm.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.6 no.3
• /
• pp.211-221
• /
• 1982

Numerical analysis has been made on the heat transfer of two dimensional turbulent channel with a slat type blockage. Especially the effects of the height of slat and Reynolds number on the heat transfer characteristics of channel wall have been investigated. The methods of accelerating the convergence of the numerical solution of governing differential equation have been also examined. Line-by-line iterative method shows higher convergence rate than point-by-point iterative method for solution of both momentum equation and energy equation. The results show that the ratio of heat transfer coefficient of the wall near the blockage to that of the fully developed flow increase with increasing the ratio of blockage to channel height and decreasing the Reynolds number. These trends of variation of heat transfer coefficient with respect to the height of slat and Reynolds number agree with those of Sparrow's experiment on the pipe flow with slat type blockage.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2008.11a
• /
• pp.401-406
• /
• 2008

This paper proposes a method to calculate the static characteristics of the FDBs with the curved surface. The general Reynolds equations are derived for the curved surfaces in the ${\theta}s$ plane. And the Reynolds equation is transformed to the finite element equations by considering the continuity of pressure and flow at the interface between the curved, journal and the thrust bearings. It also includes the Reynolds boundary condition in the numerical analysis to simulate the cavitation phenomenon. The static characteristics of the coupled journal and conical bearings were investigated due to the variation of conical angle. It shows that the conical angle is one of the important design parameters affecting the static and dynamic characteristics of FBBs.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.5
• /
• pp.685-697
• /
• 1998

Series of recent k-.epsilon. model modification have been carried out with the aid of DNS data to include the effect of near wall. Though these methods opened new way of turbulence modelings, newly developed turbulence models of its kind had yet shortcomings in prediction for the turbulent flows with various Reynolds numbers and various geometric conditions. As a remedy for these shortcomings, a new k-.epsilon. model proposed here by improving the dissipation rate equation and the damping function for eddy viscosity model. The new dissipation rate equation was modeled based on the energy spectrum and magnitude analysis. The damping function for eddy viscosity was also formulated on the ground of distribution of dissipation rate length scales near a wall and the DNS data. The new k-.epsilon. model was applied to the fully developed turbulent flows in a channel and a pipe with a wide range of Reynolds numbers. Prediction results showed that the present model represents properly the turbulence properties in all turbulent regions over a wide range of Reynolds numbers.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2007.11a
• /
• pp.840-846
• /
• 2007

This paper investigates the dynamic behavior of a HDD spindle system due to the imperfect roundness of a rotating shaft. The shaft of a spindle motor rotates with eccentricity by the unbalanced mass of the rotating part. The eccentricity generates the run-out of a spindle motor which results in the eccentric motion of a rotating part. Roundness of a shaft affects this motion which limits the memory capacity of a HDD. This research proposes a modified Reynolds equation for the coupled journal and thrust FDBs to include the variable film thickness due to the roundness. Finite element method is used to solve the Reynolds equation for the pressure distribution. Reaction forces and friction torque are obtained by integrating the pressure and shear stress, respectively. The dynamic behavior is determined by solving the equations of a motion of a HDD spindle system in six degrees of freedom with the Runge-Kutta method to characterize the motion of a rotating part. This research shows that the roundness of a rotating shaft causes the excitation frequency with integer multiple of a rotating frequency.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2007.05a
• /
• pp.852-858
• /
• 2007

This paper investigates the dynamic characteristics of a tilted HDD spindle system with fluid dynamic bearings (FDBs). Tilting motion of a HDD spindle system may be caused by improper manufacturing tolerance, such as imperfect cylindricity between shaft and sleeve of FDBs, imperfect perpendicularity between shaft and thrust as well as the gyroscopic moment of the unbalanced mass of the rotating part. Tilting motion may result in the instability of the HDD spindle system and it may increase the disk run-out to limit memory capacity. This research proposes a modified Reynolds equation for the coupled journal and thrust FDBs to include the variable film thickness due to the cylindricity and the perpendicularity. Finite element method is used to solve the Reynolds equation for the pressure distribution. Reaction forces and friction torque are obtained by integrating the pressure and shear stress, respectively. The dynamic behavior is determined by solving the equations of a motion of a HDD spindle system in six degrees of freedom with the Runge-Kutta method to study whirling and tilting motions. This research shows that the cylindricity and the perpendicularity increase the tilting angle and whirl radius of the rotor.

• Wind and Structures
• /
• v.19 no.6
• /
• pp.687-708
• /
• 2014

The performance of the $k-{\varepsilon}$ and $k-{\omega}$ two-equation turbulence models was investigated in computational simulations of the neutrally stratified atmospheric boundary layer developing above various terrain types. This was achieved by using a proposed methodology that mimics the experimental setup in the boundary layer wind tunnel and accounts for a decrease in turbulence parameters with height, as observed in the atmosphere. An important feature of this approach is pressure regulation along the computational domain that is additionally supported by the nearly constant turbulent kinetic energy to Reynolds shear stress ratio at all heights. In addition to the mean velocity and turbulent kinetic energy commonly simulated in previous relevant studies, this approach focuses on the appropriate prediction of Reynolds shear stress as well. The computational results agree very well with experimental results. In particular, the difference between the calculated and measured mean velocity, turbulent kinetic energy and Reynolds shear stress profiles is less than ${\pm}10%$ in most parts of the computational domain.

• Proceedings of the Korea Technical Association of the Pulp and Paper Industry Conference
• /
• 2002.11a
• /
• pp.81-90
• /
• 2002

Air entrainment of a winding roll with a nip roller was studied numerically. The amount of air entrainment between two rotating rollers was obtained by solving lubrication equation, Reynolds equation, which neglect the existence of a web. However, the numerical model of this study included the web existence, therefore it considered the two lubricating air films between a winding roll and a web and also between a nip roller and the web. The pressure profiles and gap profiles of the two films were obtained by solving lubrication equation for the two air films and force balance equation of the web. Ballooning phenomenon was examined in terms of nip force, wrap angle, web stiffness, web speed, and web tension. This ballooning phenomenon caused by the back flow of the air film blocked by the nip roller. Air entrainment of the two numerical models was compared.