• International Journal of Fluid Machinery and Systems
• /
• v.5 no.2
• /
• pp.60-64
• /
• 2012

The helical flow regime was investigated by using DPIV when the rotating Reynolds number is small. The wall slits were azimuthally located along the inner wall of outer cylinder and the slits number of each model was 9 and 18, another plain wall model was also studied for comparison purpose. The helical vortex flow regime can be observed in all the three models. The negative temperature gradients determine the direction of the rotation and movement of the helical vortex. But the helical wavy vortex flow can only be found in the plane and 9-slit models. And the result showed that the existence of slit wall accelerated the transition process.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.24 no.2
• /
• pp.242-254
• /
• 2016

This paper is the forth investigation on the evaluation methods of flow characteristics in a steady flow bench. In the previous works, it was concluded that the assumption of the solid rotation might cause serious problems and both of the eccentricity and the velocity profile distort the flow characteristics when using the ISM at 1.75B plane. Also particle image velocimetry (PIV) measurement at this position showed that the real velocity profile was far from the assumption of ISM evaluation. In this paper, the planar velocity profiles were measure from 1.75B to 6.00B position by PIV and the characteristics were examined according to the valve angles and lifts for further investigations about the effect of the position on the velocity profile. The results show that $26^{\circ}$ valve angle is always an unique exceptional case in all aspects. If the valve angle is $21^{\circ}$ and below, the planar velocity profiles according to the lift and the position are similar to each other, however, the tangential velocity curves along with the radial direction have common tendencies up to $16^{\circ}$ angle. Also the well arranged swirl behaviors are generally observed at the position above 3.00B and the velocity contour lines come closer to the concentric circle as the valve lift increases. In addition, the gradient of tangential velocity along with the radial direction from the swirl center becomes stable and constant as the position goes downstream. Concurrently the velocity gradient is larger to the eccentric direction of the center. In the meantime the tangential velocity curves along with the radial direction are irregular and various at 1.75B, however, they become regular and reach higher level as the evaluation position goes downstream. At this time the curves of 4.50B are the best fitted to the ideal one. On the other hand in an exceptional case, $26^{\circ}$, the velocity contours are very complicated over 6mm valve lift regardless the position and the gradient increases to the opposite direction of the eccentric center. Also, 6.00B is a best fitting position in the geometrical cylinder center base. With respect to the swirl center, the distribution range of centers for 1.75B is different to that for the other positions and the eccentricities of this plane are larger regardless the valve angle. After 1.75B, there is no certain tendency in the center position change according to the valve angle and lift. Additionally, the eccentricities are not sufficiently small to neglecting the effect on ISM measurement.

• Structural Engineering and Mechanics
• /
• v.20 no.4
• /
• pp.435-450
• /
• 2005

Here, the dynamic instability characteristics of aero-thermo-mechanically stressed functionally graded plates are investigated using finite element procedure. Temperature field is assumed to be a uniform distribution over the plate surface and varied in thickness direction only. Material properties are assumed to be temperature dependent and graded in the thickness direction according to simple power law distribution. For the numerical illustrations, silicon nitride/stainless steel is considered as functionally graded material. The aerodynamic pressure is evaluated based on first-order high Mach number approximation to the linear potential flow theory. The boundaries of the instability region are obtained using the principle of Bolotin's method and are conveniently represented in the non-dimensional excitation frequency-load amplitude plane. The variation dynamic instability width is highlighted considering various parameters such as gradient index, temperature, aerodynamic and mechanical loads, thickness and aspect ratios, and boundary condition.

• 한국해양학회지
• /
• v.30 no.4
• /
• pp.355-364
• /
• 1995

A heavy fluid is injected to a rotating cylindrical container of flat or inclined bottom filled with homogeneous lighter fluid. Continuous flow-in and spreading patterns over the bottom of the container are observed and at the same time upper-layer motions induced by the movement of the heavy fluid are traced by thymol blue solution. Regardless of bottom geometry, the injected denser fluid is deflected toward "western wall" and continuous its path along the boundary with radial spreading which occurs in the bottom boundary layer to make a quite asymmetric flow. When the bottom contains a slope(${\beta}$-plane), increased pressure gradient causes the fluid move faster to produce a stronger Coriolis force. This makes the width of the flow narrower than that of f-plane. But, when the denser flow reaches the southern part of the container, a local-depth of denser fluid increases (much greater than the Ekman-layer depth) such that the spreading velocity along the wall is reduced and the interfacial slope increases to make the upper-layer adjust geographically to have oppositely directed upper-layer motion along the interfacial boundary. The role of the denser fluid in terms of vorticity generation in the upper-layer is such that it produces local topographic effect over the western half of the container and also induces vortex-tube stretching which is especially dominant in the f-plane.

• Journal of the korean Society of Automotive Engineers
• /
• v.14 no.4
• /
• pp.39-49
• /
• 1992

Three dimensional numerical calculation carried out to investigate the eccentricity effect of intake valve on the in-cylinder flow fields for the intake stroke and the compression stroke. During the intake stroke, a corner vortex in the vicinity of the valve exit interacted strongly with a toroidal vortex in the case of axisymmetric valve. But a weak interaction between the corner vortex and the toroidal vortex occurred due to the eccentricity of the valve in the narrow region between valve and cylinder wall in the case of offset valve. During the compression stroke, it was found that a solid body rotation was maintained in the radial-circumferential plane in the case of axisymmetric valve. But a weak secondary vortex was formed in the radial-circumferntial plane in the case of offset valve, because of the interaction between swirl flows and inward flows towards cylinder axis. The calculated turbulence intensity presented a similar trend with the experiental results but, in spite of using the modified k-.epsilon. model, it was found that the qualitative difference between the numerical results and experimental results was large in the region where the velocity gradient is substantial.

• The Korean Journal of Rheology
• /
• v.10 no.4
• /
• pp.202-216
• /
• 1998

Improved version of previous 'Orthotropic' closure approximation, termed 'ORW' has been numerically developed using new homogeneous flow data. Previous 'Orthotropic' closure approximation, i.e., ORF or ORL showed non-physical oscillation for interaction coefficient $C_1$<0.001 at simple shear flow. It also shows non-physcial oscillation and under-prediction compared with 'Distribution Function Calculation' at non-homogeneous flow of center-gated disk. These phenomena are mainly due to the flow data of 'Distribution Function Calculation' which were used for least-square optimization. ORW obtained by fitting flow data of low interaction coefficient does not show non-physical oscillation and results in reasonably good behaviors at non-homogeneous flows as well as homogeneous flows. Fitting function forms have not been found to improve overall behaviors. It has been found that considering all the eigenvalues of orientation tensor (including the third eigenvalues) might end up with a better closure approximation than just considering the first and second eigenvalues. It is, however, very important and yet difficult to select appropriate function forms of eigenvalues. Numerical simulation including coupling and in-plane velocity gradient effects were performed for injection mold filing process with a film-gated strip and a center-gated disk using ORW and various other closure approximations for comparisons. Although ORW is in excellent agreement with 'Distribution Function Calculation', the predicted results seem to have consistent error in comparison with experimental data. The diffusivity term with constant interaction coefficient might have to be further investigated in order to accurately describe orientation states.

• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.4 no.1
• /
• pp.46-56
• /
• 1994

The characteristics of flows, temperatures, concentrations of the boron are numerically studied when uniform axial magnetic fields are applied in the Czechralski crucible. The to governing factors to the flow regimes are buoyancy, thermocapillarity, centrifugal forces, magnetic forces, diffusion coefficient and segregation coefficient of the boron. Since the concentration of the boron is so low that buoyancy effects are negligible, it cannot affect the flow and temperature fields. From the fact that the flow fields are rotationally symmetric, two velocity components in the meridional plane and the circumferential velocity are calculated together with the temperature in the steady state. Based on the known velocity and temperature distributions the unsteady concentration distributions of the boron are calculated. As the strength of the magnetic is increased, the flow velocities are decreased. Circumferential velocities are large near the crucible side-wall and in the region below the rotating crystal. Steep temperatures gradient near the edge of the rotating crystal causes the Marangoni convection. It has been found out that the convection characteristics affects the unsteady transport phenomena of the boron.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.11 s.242
• /
• pp.1265-1276
• /
• 2005

The elliptic conceptual second moment models for turbulent heat fluxes, which are proposed on the basis of elliptic-blending and elliptic-relaxation equations, are applied to calculate the combined forced and natural turbulent convection in a vertical plane channel. The models satisfy the near-wall balance between viscous diffusion, viscous dissipation and temperature-pressure gradient correlation, and also have the characteristics of approaching its respective conventional high Reynolds number model far away from the wall. Also the models are closely linked to the elliptic blending model which is used for the prediction of Reynolds stress. In order to calibrate the heat flux models, firstly, the distributions of mean temperature and scala flux in fully developed channel flow with constant wall difference temperature are solved by the present models. The buoyancy effect on the turbulent characteristics including the mean velocity and temperature, the Reynolds stress tensor, and the turbulent heat flux vector are examined. In the opposing flow, the turbulent transport is greatly enhanced with both the Reynolds stresses and the turbulent heat fluxes being remarkably increased; whereas, in the aiding flow, the opposite change is observed. The results of prediction are directly compared to the DNS to assess the performance of the model predictions and show that the behaviors of the turbulent heat transfer in the whole flow region are well captured by the present models.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.1
• /
• pp.262-273
• /
• 1991

In the present study, the steady, incompressible, isothermal, developing flow in a 90.deg. rectangular cross sectional strongly curved duct with aspect ratio 1:1.5 and Reynolds number of 9.4*10$^{4}$ has been investigated. Measurements of components of mean velocities, pressures, and corresponding components of the Reynolds stress tensor are obtained with a hot-wire anemometer and pitot tube. In general, flow in a curved duct is characterized by the secondary vortices which are driven mainly by centrifugal force-radial pressure gradient imbalance, and the stress field stabilizing effects near the convex wall and destablizing effects close to the concave wall. It was found that the secondary mean velocities attain values up to 39% of the bulk velocity and are largely responsible for the convections of Reynolds stress in the cross stream plane. Therefor upstream of the bend the Reynolds stress are low. Corresponding to the small boundary layer thickness. At successive planes, large values of Reynolds stress were observed near the concave surface and the side wall.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.9
• /
• pp.688-698
• /
• 2009

Using the multi-dimensional, multi-phase, nonisothermal Polymer Electrolyte Fuel Cell (PEFC) model presented in Part I, the effects of land/channel flow-field on temperature and liquid saturation distributions inside PEFCs are investigated in Part II. The focus is placed on exploring the coupled water transport and heat transfer phenomena within the nonisothermal and two-phase zone existing in the diffusion media (DM) of PEFCs. Numerical simulations are performed varying the land and channel widths and simulation results reveal that the water profile and temperature rise inside PEFCs are considerably altered by changing the land and channel widths, which indicates that oxygen supply and heat removal from the channel to the land regions and liquid water removal from the land toward the gas channels are key factors in determining the water and temperature distributions inside PEFCs. In addition, the adverse liquid saturation gradient along the thru-plane direction is predicted near the land regions by the numerical model, which is due to the vapor-phase diffusion driven by the temperature gradient in the nonisothermal two-phase DM where water evaporates at the hotter catalyst layer, diffuses as a vapor form and then condenses on the cooler land region. Therefore, the vapor phase diffusion exacerbates DM flooding near the land region, while it alleviates DM flooding near the gas channel.