• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.103-106
• /
• 2002

This paper presents an effective numerical method for analyzing three-dimensional unsteady conjugate heat transfer problems of a curved pipe subjected to infernally thermal stratification. In the present numerical analyses, the thermally stratified flows in the pipe are simulated using the standard $k-{\varepsilon}$turbulent model and the unsteady conjugate heat transfer is treated numerically with a simple and convenient numerical technique. The unsteady conjugate heat transfer analysis method is implemented in a finite volume thermal-hydraulic computer code based on a non-staggered grid arrangement, SIMPLEC algorithm and higher-order bounded convection scheme. Numerical calculations have been performed far the two cases of thermally stratified pipe flows where the surging directions are opposite each other i.e. In-surge and out-surge. The results show that the present numerical analysis method is effective to solve the unsteady flow and conjugate heat transfer in a curved pipe subjected to infernally thermal stratification.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.3 no.3
• /
• pp.186-196
• /
• 1991

Numerical analysis has been performed to investigate the effects of the turbulence intensity and Prandtl number on the local heat transfer around a circular cylinder in crossflow. The governing equations were reformulated in a non-orthogonal coordinate system with Cartesian velocity components and discretised by the finite volume method with a non-staggered variable arrangement. For laminar flow, the calculations were performed for the Reynolds numbers 26 and 200. The results showed good agreement with the experimental results. For turbulent flow of the Reynolds number $1{\times}10^5$ and $2{\times}10^6$, the results showed that with an increase in the turbulent intensity in the main stream, the local Nusselt number increases in the front region of the circular cylinder. But the effect of turbulent intensity on the local Nusselt number diminishes in the wake region. The influence of Prandtl numbers show similar trend to that of turbulent intensity.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.474-481
• /
• 2001

This paper addresses a numerical method for predicting transient temperature distributions in the wall of a curved pipe subjected to internally thermal stratification flow. A simple and convenient numerical method of treating the unsteady conjugate heat transfer in the non-orthogonal coordinate systems is presented. The proposed method is implemented in a finite volume thermal-hydraulic computer code based on a cell-centered, non-staggered grid arrangement, the SIMPLEC algorithm, a higher-order bounded convection scheme, and the modified version of momentum interpolation method. Calculations are performed for the transient evolution of thermal stratification in two curved pipes, where the one has thick wall and the other has so thin wall that its presence can be negligible in the heat transfer analysis. The predicted results show that the thermally stratified flow and transient conjugate heat transfer in a curved pipe with a finite wall thickness can be satisfactorily analyzed by the present numerical method, and that the neglect of wall thickness in the prediction of pipe wall temperature distributions can provide unacceptably distorted results.

• Proceedings of the KSME Conference
• /
• 2000.04a
• /
• pp.900-907
• /
• 2000

The three-dimensional turbulent flow in curved pipes susceptible to flow-accelerated corrosion has been analyzed numerically to predict the pressure and shear stress distributions on the inner surface of the pipes. The analysis employs the body-fitted non-orthogonal curvilinear coordinate system and a standard $ {\kappa}-{\varepsilon}$ turbulence model with wall function method. The finite volume method is used to discretize the governing equations. The convection term is approximated by a high-resolution and bounded discretization scheme. The cell-centered, non-staggered grid arrangement is adopted and the resulting checkerboard pressure oscillation is prevented by the application of a modified version of momentum interpolation scheme. The SIMPLE algorithm is employed for the pressure and velocity coupling. The numerical calculations have been performed for two curved pipes with different bend angles and curvature radii, and discussions have been made on the distributions of the primary and secondary flow velocities, pressure and shear stress on the inner surface of the pipe to examine applicability of the present analysis method. As the result it is seen that the method is effective to predict the susceptible systems or their local areas where the fluid velocity or local turbulence is so high that the structural integrity can be threatened by wall thinning degradation due to flow-accelerated corrosion.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.12
• /
• pp.1802-1812
• /
• 2001

The purpose of this study is to analyze numerically the movement of particles included in turbulent fluid flow characteristics of metallic surfaces. To describe fluid flew, the incompressible Navier-Stokes equation discretized by the finite volume method were solved on the non-orthogonal coordinates with non-staggered variable arrangement, and the k-$\xi$ turbulence model was adapted. After fluid flow was calculated, particle movement was predicted from the Lagrangian approaches. Non-essential complexities were avoided by assuming that the particles had spherical shapes and the Stoke's drag formula only consisted of external farces acting upon them. In order to validate the numerical calculations, the results were compared with the experimental data reported in literature and agreed well with them. The drag force coefficient equation showed better agreement with the experimental data in the prediction of particle movement than the correction factor equation. Impact velocity and impact angle increased as inlet turbulence intensity decreased, relative jet height was lower. or the Reynolds number was larger.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.7
• /
• pp.849-861
• /
• 1997

In this study, the effect of hole geometry of the cooling system on the flow and temperature field was numerically calculated. The finite volume method was employed to discretize the governing equation based on the non-orthogonal coordinate with non-staggered variable arrangement. The standard k-.epsilon. turbulence model was used and also the predicted results were compared with the experimental data to validate numerical modeling. The predicted results showed good agreement in all cases. To analyze the effect of the discharge coefficient for slots of different length to width, the inlet chamfering and radiusing holes were considered. The discharge coefficient was increased with increment of the chamfering ratio, radiusing ratio and slot length to width and also the effect of radiusing showed better result than chamfering in all cases. In order to analyze the difference between the predicted results with plenum region and without plenum region, the velocity profiles of jet exit region for a various flow conditions were calculated. The normal velocity components of jet exit showed big difference for the low slot length to width and high blowing rate cases. To analyze the flow phenomena injected from a row of inclined holes in a real turbine blade, three dimensional flow and temperature distribution of the region including plenum, hole and cross stream with flow conditions were numerically calculated. The results have shown three-dimensional flow characteristics, such as the development of counter rotating vortices, jetting effect and low momentum region within the hole in addition to counter rotating vortex structure in the cross stream.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.12
• /
• pp.4013-4026
• /
• 1996

A numerical simulation has been carried out for the jet impinging on a flat plate and a semi-circular concave surface. In this computation finite volume method was employed to solve the full Navier-Stokes equation based on a non-orthogonal coordinate with non staggered variable arrangement. The standard k-.epsilon. turbulent model and low Reynolds number k-.epsilon. model(Launder-Sharmar model) with Yap's correction were adapted. The accuracy of the numerical calculations were compared with various experimental data reported in the literature and showed good predictions of centerline velocity decay, wall pressure distribution and skin friction. For the jet impingement on a semi-circular concave surface, potential core length was calculated for two different nozzle(round edged nozzle and rectangular edged nozzle) to consider effects of the nozzle shape. The result showed that round edged nozzle had longer potential core length than rectangular edged nozzle for the same condition. Heat transfer rate along the concave surface with constant heat flux was calculated for various nozzle exit to surface distance(H/B) in the condition of same jet velocity. The maximum local Nusselt number at the stagnation point occurred at H/B = 8 where the centerline turbulent intensity had maximum value. The predicted Nusselt number showed good agreement with the experimental data at the stagnation point. However heat transfer predictions along the downstream were underestimated. This results suggest that the improved turbulence modeling is required.

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

The numerical simulation has been conducted for the investigation of flow and heat transfer characteristics of an oblique impingement jet injected to a flat plate. The finite volume method was used to discretize the governing equations based on the non-orthogonal coordinate with non-staggered variable arrangement. The $textsc{k}$-$\varepsilon$-ν(sup)'2 turbulence model was employed to consider the consider the anisotropic flow characteristics generated by the impingement jet flow. The predicted results were compared with the experimental data and those of the standard $textsc{k}$-$\varepsilon$ turbulence model. The results of the $textsc{k}$-$\varepsilon$-ν(sup)'2 model showed better agreement with the experimental data than those of the standard $textsc{k}$-$\varepsilon$ model. In order to get the optimum condition, the flow and temperature fields were calculated with a variation of inclined angle($\alpha$=30$^{\circ}$~90$^{\circ}$) and the distance between the jet exit and impingement plate-to-diameter (L/D=4~10) at a fixed Reynolds number(Re=20,000). For a small L/D, the near-peak Nusselt numbers were not significantly effected by the inclined angle. The near-peak Nusselt numbers were not significantly affected by the L/D in the case of a large $\alpha$. The overall shape of the local Nusselt numbers was influenced by both the jet orifice-to-plate spacing and the jet angle.