• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.7
• /
• pp.1675-1683
• /
• 1995

A numerical study of the fluid flow and heat transfer characteristics of the two-dimensional impingement jet with a confinement plate has been carried out. The fluid flow was calculated by solving the full Navier-Stokes equation. In doing that, the well known SIMPLER algorithm was used and the trouble making convection term was discretized according to QUICKER scheme. The energy equation was simply solved by using the SOR method. For the Reynolds number of 10000, two channel heights, say 1.5 and 3.0 times the jet exit width, and two thermal boundary conditions constant wall temperature and constant wall heat flux were considered. Discrete heat sources were flush mounted along the impingement plate at a distance of 0, 2, 3, 4, 5, 6, 10, 12, times the jet exit width from the stagnation point. The length of each heat source is 4 times the jet exit width long. The Nusselt number averaged over each heat source was compared with experiment. Comparison shows that both calculations and experiment have the secondary peak of Nusselt number at downstream of stagnation point, even though there is a little quantitative difference in between. The difference is believed due to abscure thermal boundary condition in experiment and also accuracy of turbulence model used. The secondary peak is shown to be caused by rigorous turbulent flow motion generated as the wall jet flow is retarded and developes into the channel flow without flow reversal.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.2
• /
• pp.261-268
• /
• 2002

Turbulent flow past a square cylinder confined in a channel is numerically investigated by Large Eddy Simulation(LES). The main objectives of this study are to verify the experimental results of Nakagawa et al.[Exp. in Fluids, Vol. 27, 3, pp. 284∼294, 1999] by LES and to obtain related flow information in detail. The LES results obtained are in excellent agreement with the experiment both qualitatively and quantitatively. The passive paticles numerically released into the flow field clearly show the barman vortex street. However, the vortices shed from the cylinder are significantly affected by the presence of the plates. Futhermore, periodic and alternating vortex-rollups are observed in the vicinity of the plates. The rolled-up vortex is convected downstream together with the corresponding Karman vortex forming a counter-rotating vortex pair. It is also revealed that the cylinder greatly enhances mixing process of the flow.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.6
• /
• pp.874-887
• /
• 1998

A 2-frame PTV (particle tracking velocimetry) system using the concept of match probability between two consequent image frames has been developed to obtain instantaneous velocity fields. The overall 2-frame PTV system including image pre-processing, tracking algorithm and post-processing routine was implemented to apply to real flows. The developed 2-frame PTV system has several advantages such as high recovery ratio of velocity vectors, low error ratio and small computational time compared with the conventional 4-frame PTV and the FFT-based cross-correlation PIV technique. The 2-frame PTV system was applied to a turbulent channel flow over a rectangular block to check its reliability and usefulness. Total 96 sequential image frames have been captured and processed to get both mean and fluctuating velocity vector fields over the recirculating region. The mean velocity and turbulent intensity profiles were well agreed with hte LDV measurements in the separated region behind the block. Time-averaged reattachment length is about 6.3 times of the block height.

• Journal of computational fluids engineering
• /
• v.13 no.3
• /
• pp.28-34
• /
• 2008

In engineering application of large eddy simulation, there are still questions as follows grid dependency on numerical results, the effect of upwind scheme against a calculation instability, appropriate boundary conditions dealing with turbulence fluctuation and the performance of SGS models. In this study, in order to develop the LES to the engineering application, large eddy simulation was carried out to investigate the effect of upwind scheme, turbulent subgrid model and the grid dependancy of the flow around a wall-mounted cube in a channel at Re=40,000 based on cubic height and bulk mean velocity. The computed velocities, turbulence quantities, separation and reattachment length were evaluated compared with the experimental results of R. Matinuzzi and C. Tropea.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.9
• /
• pp.1351-1359
• /
• 2003

This study presents a numerical procedure to optimize shape of streamwise periodic ribs mounted on both of the principal walls to enhance turbulent heat transfer in a rectangular channel flow. The response surface method is used as an optimization technique. The optimization is based on Navier-Stokes analysis of flow and heat transfer with $k-{\varepsilon}$ turbulence model. The width-to-height ratio of a rib, rib height-to-channel height ratio and rib pitch to rib height ratio are chosen as design variables. The object function is defined as a function of heat transfer coefficient and friction drag coefficient with weighting factor. Optimum shapes of the rib have been obtained for the range of 0.02 to 0.1 of weighting factor.

• 한국전산유체공학회:학술대회논문집
• /
• 2002.05a
• /
• pp.128-133
• /
• 2002

The heat transfer characteristics of a turbulent flow in a ribbed two-dimensional channel have been investigated numerically. The fully elliptic governing equations, coupled with a four-equation turbulence model, $\kappa-\omega-\bar{t^2}-\epsilon_t$, are solved by a finite volume method of SIMPLE type. Calculations have been carried out for three rib cross-sections : square, triangular, and semicircular, with various rib pitches and Reynolds numbers. The procedure appears to be satisfactory as the results for the square rib compare favorably with available experimental data and earlier calculation. The optimal rib pitch that yields the maximum heat transfer has been identified. It is also found that the square rib is most effective in enhancing the heat transfer. The semicircular rib, on the other hand, incurs the least amount of pressure drop but the improvement in heat transfer is substantially lower.

• Journal of the Society of Naval Architects of Korea
• /
• v.40 no.3
• /
• pp.1-8
• /
• 2003

For the experiment of the friction drag reduction by microbubble injection, a drag reduction water tunnel was specifically designed and made. Experimental apparatus and procedures were devised and developed for measuring the change of wall friction drag with microbubble injection. For fully-developed channel flows. the change of friction drag with important parameters of microbubble injection is investigated and the experimental data and results obtained are presented. The amount of friction drag reduction up to 25% is observed in the present study.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.3 no.2
• /
• pp.1-9
• /
• 1999

Applicability of the pressure gradient method which is formulated based on pressure gradient is verified against turbulent flow analysis. In the pressure gradient method, pressure gradient instead of pressure itself is obtained using continuity constraint. Since correct pressure gradient is found only when mass conservation is satisfied, pressure gradient method can reflect physics of flow field properly The pressure gradient method is formulated with semi-staggered grid system which locates each primitive variables on the same grid point but evaluates pressure gradient in-between. This grid system ensures easy programming and reflection of correct physics in analysis. For verifying applicability of this method, the pressure gradient method is applied to turbulent flow analysis with low Reynolds number $\kappa$-$\varepsilon$ model. Turbulent flows include fully developed channel flow, backward-facing step flow, and conical diffuser flow. Prediction results show that the pressure gradient method can be applied to turbulent flow analysis. However, the pressure gradient method requires somewhat long computation time. Proper way to find optimum under-relaxation factor, $\gamma$, is also need to be developed.

• Journal of the Korean Society of Visualization
• /
• v.17 no.3
• /
• pp.12-18
• /
• 2019

Direct numerical simulation of a fully developed turbulent plane Couette-Poiseuille flow with a two-dimensional (2-D) rod-roughened wall is performed to investigate the impacts of the surface roughness. It is shown that the logarithmic region in the mean velocity profile over the rough wall Couette-Poiseuille flow is significantly shortened by the surface roughness compared to that over a turbulent Couette-Poiseuille flow with smooth wall. The Reynolds shear stress over the rough wall Couette-Poiseuille flow is decreased compared to that for a smooth case in the outer layer. These results are attributed to weakened turbulence activity or roll-cell mode over the rough wall Couette-Poiseuille flow near the channel centerline due to suppressed development of u'-structure on the top wall, as documented through spanwise energy spectra of the streamwise velocity fluctuations. Inspection of congregation motion near the bottom wall and time evolution of u'-structure reveal weakened co-supporting cycle for the rough wall case.

• Journal of Mechanical Science and Technology
• /
• v.20 no.1
• /
• pp.147-157
• /
• 2006

This paper represents a numerical study of the flow field due to the interactions between a pair of vortices produced by vortex generators in a rectangular channel flow. In order to analyze longitudinal vortices induced by the vortex generators, the pseudo-compressibility method is introduced into the Reynolds-averaged Navier-Strokes equations of a 3-dimensional unsteady, incompressible viscous flow. A two-layer $k-{\epsilon}$ turbulence model is applied to a flat plate 3-dimensional turbulence boundary to predict the flow structure and turbulence characteristics of the vortices. The computational results predict accurately the vortex characteristics related to the flow field, the Reynolds shear stresses and turbulent kinetic energy. Also, in the prediction of skin friction characteristics the computational results are reasonably close to those of the experiment obtained from other researchers.