• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.8
• /
• pp.1128-1138
• /
• 2000

Pumping action in ejector systems is generally achieved through the mixing of a high-velocity and high-energy stream with a lower-velocity and lower-energy stream within a duct. The design and performance evaluation of the ejector systems has developed as a combination of scale-model experiments, empiricism and theoretical analyses applicable only to very simplified configurations, because of the generic complexity of the flow phenomena. In order to predict the detailed performance characteristics of such systems, the flow phenomena throughout the operating regimes of the ejector system should be fully understood. This paper presents the computational results for the two-dimensional supersonic ejector system with a second throat. The numerical simulations are based on a fully implicit finite volume scheme of the compressible Reynolds-averaged Navier-Stokes equation in a domain that extends from the stagnation chamber to the diffuser exit. For a wide range of the operating pressure ratio the flow field inside the ejector system is investigated in detail. The results show that the supersonic ejector systems have an optimal throat area for the operating pressure ratio to be minimized.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.145-155
• /
• 1998

A computer code has been developed for the computation of the viscous flow around a ship model with the free surface. In this code, the incompressible Reynolds-averaged Navier-Stokes equations are solved numerically by a finite difference method which employes second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration of the governing equations. For the turbulence closure, a modified version of the Baldwin-Lomax model is exploited. The location of the free surface is determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and the boundary-fitted grid is generated at each time step so that one of the grid surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition is applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method and the computer code developed in the present study, the numerical computations are carried out for both Wigley parabolic hull and Series 60 $C_B=0.6$ ship model and the computational results are compared with the experimental data.

• Journal of computational fluids engineering
• /
• v.18 no.1
• /
• pp.1-6
• /
• 2013

Missile and fighter aircraft have been challenged by low restoring nose-down pitching moment at high angle of attach. The consequence of weak nose-down pitching moment can be resulting in a deep stall condition. Especially, the pressure oscillation has a huge effect on noise generation, structure damage, aerodynamic performance and safety, because the flow has strong unsteadiness at high angle of attack. In this paper, the unsteady aerodynamics coefficients were analyzed at high angle of attack up to 50 degrees around two dimensional NACA0012 airfoil. The two dimensional unsteady compressible Navier-Stokes equation with a LES turbulent model was calculated by OHOC (Optimized High-Order Compact) scheme. The flow conditions are Mach number of 0.3 and Reynolds number of $10^5$. The lift, drag, pressure, entropy distribution, etc. are analyzed according to the angle of attack. The results of average lift coefficients are compared with other results according to the angle of attack. From a certain high angle of attack, the strong vortex formed by the leading edge are flowing downstream as like Karman vortex around a circular cylinder. The primary and secondary oscillating frequencies are analyzed by the effects of these unsteady aerodynamic characteristics.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.484-489
• /
• 2003

The mass flow rate of gas flow through critical nozzle depends on the nozzle supply conditions and the cross-sectional area at the nozzle throat. In order that the critical nozzle can be operated at a wide range of supply conditions, the nozzle throat diameter should be controlled to change the flow passage area. This can be achieved by means of a variable critical nozzle. In the present study, both experimental and computational works are performed to develop variable critical nozzle. A cone-cylinder with a diameter of d is inserted into conventional critical nozzle. It can move both upstream and downstream, thereby changing the cross-sectional area of the nozzle throat. Computational work using the axisymmetric, compressible Navier-Stokes equations is carried out to simulate the variable critical nozzle flow. An experiment is performed to measure the mass flow rate through variable critical nozzle. The present computational results are in close agreement with measured ones. The boundary layer displacement and momentum thickness are given as a function of Reynolds number. An empirical equation is obtained to predict the discharge coefficient of variable critical nozzle.

• Tribology and Lubricants
• /
• v.30 no.2
• /
• pp.77-85
• /
• 2014

This study set out to predict the load capacity and rotordynamic coefficients of tilting-pad journal bearings, taking the pivot stiffness into account. The analysis uses rocker-back (cylindrical) and ball in socket (spherical) pivot models, both of which are based on Hertzian contact stress theory. The models ascertain the non-linear elastic deformation of the pivots according to the applied load, pivot geometry, and material properties. At present, the Reynolds equation for an isothermal, isoviscous, and incompressible fluid is used to calculate the film pressure by using the finite-element method, after which the Newton-Raphson method is used to simultaneously find the journal center location, pad angles, and pivot deflections. The bearing analysis, excluding the pivot models, is validated using predictions those are readily available in the literature. As the rotor speed increases, the predicted journal eccentricity and damping coefficients decrease, but the stiffness coefficients increase, as expected. Most importantly, the implementation of the pivot models increases the journal eccentricity but significantly decreases the stiffness and damping coefficients of the tilting-pad journal bearings.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2003.11a
• /
• pp.556-563
• /
• 2003

This paper presents a method to investigate the characteristics of a hydrodynamic bearing of a HDD spindle motor due to elevated temperature considering the variation of the clearance as well as the lubricant viscosity. Iterative finite element analysis of the heat conduction and the thermal deformation is performed to determine the viscosity and clearance of a hydrodynamic bearing due to elevated temperature until the temperature of the bearing area converges. Proposed method is verified by comparing the calculated temperature with the measured one in elevated surrounding temperature as well as in room temperature. This research shows that elevated temperature changes the clearance as well as the lubricant viscosity of the hydrodynamic bearing of a HDD spindle motor. Once the viscosity and clearance of a hydrodynamic bearing of a HDD spindle motor are determined, finite element analysis of the Reynolds equation is performed to investigate the static and dynamic characteristics of a hydrodynamic bearing of a HDD spindle motor due to elevated temperature. It also shows that the variation of clearance due to elevated temperature is another important design consideration to affect the static and dynamic characteristics of a hydrodynamic bearing of a HDD spindle motor.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.11
• /
• pp.963-970
• /
• 2007

The objectives of this study are to examine a bubble trap mechanism of the turbidimeter for low turbidity and to acquire it's technology concerned. Reynolds-Averaged Wavier-Stokes equation and Laglangian discrete phase model were applied to analyze a flow field in the bubble trap. 3D hybrid grid system was used to simulate the flow field of bubble trap and numbers of it's node point are about 110,000. From the comparison between the standard $k-{\varepsilon}$ model and the laminar state, it was found that the former estimated less the velocity in the outlet of bubble trap than the latter did, and that the former estimated more the shear stress at the wall of bubble trap than the latter did. And, it was possible to visualize the path of bubbles in the bubble trap and to copy the removal process of bubbles out bubble trap. Also, it was found that nearly most of bubbles in the bubble trap disappeared.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.10 s.253
• /
• pp.966-972
• /
• 2006

A numerical simulation is performed fur developing turbulent flow in a strongly curved 180 deg pipe and its downstream tangent by a new solution code(PowerCFD) which adopts an unstructured cell-centered method. The governing equations are discretized as the full elliptic from of the equations of motion. Three typical two-equation turbulence models of low-Reynolds-number form are used to approximate the turbulent stress field. Solutions fur both streamwise and circumferential velocity components are compared with the experimental data by Azzola et at.(1986). The ${\kappa}-{\omega}$ model by Wilcox(1988) is found to give better prediction performance than the other two. Predicted secondary velocities and streamwise velocity component contours at sequential longitudinal stations are also presented in order to enable a detailed description of the complete flow. It is also found that, in the bend both mean streamwise and secondary velocities never achieve a fully-developed state and the code is capable of producing very well the complex nature of steady flow in a strongly curved pipe.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.39 no.3
• /
• pp.201-209
• /
• 2011

An unsteady RANS analysis study of the 3-D flow around the NREL Phase VI horizontal axis wind turbine(HAWT) was performed using sliding mesh approach. Two different analysis models such as rotor-only and rotor with tower/nacelle were constructed to investigate the blade/tower interaction. Analysis results for the rotor with tower/nacelle were compared with the corresponding NREL's experimental data which produced fairly good validation of the present CFD model. Comparison of flows around those two models also clearly showed the blade/tower interaction even it was small for upwind configuration. Other visualization results and integrated aerodynamic loads including torque of the blade demonstrated the effective unsteady flow simulation capability of the present CFD model.

• International Journal of Aeronautical and Space Sciences
• /
• v.14 no.3
• /
• pp.201-209
• /
• 2013

Much numerical and experimental research has been done for the flow around an oscillating airfoil. The main research topics are vortex shedding, dynamic stall phenomenon, MAV's lift and thrust generation. Until now, researches mainly have been concentrated on analyzing the wake flow for the variation of frequency and amplitude at a low angle of attack. In this study, wake structures and acoustic wave propagation characteristics were studied for a plunging airfoil at high angle of attack. The governing equations are the Navier-Stokes equation with LES turbulence model. OHOC (Optimized High-Order Compact) scheme and 4th order Runge-Kutta method were used. The Mach number is 0.3, the Reynolds number is, and the angle of attack is from $20^{\circ}$ to $50^{\circ}$. The plunging frequency and the amplitude are from 0.05 to 0.15, and from 0.1 to 0.2, respectively. Due to the high resolution numerical method, wake vortex shedding and pressure wave propagation process, as well as the propagation characteristics of acoustic waves can be simulated. The results of frequency analysis show that the flow has the mixed characteristics of the forced plunging frequency and the vortex shedding frequency at high angle of attack.