• Journal of Korea Water Resources Association
• /
• v.50 no.2
• /
• pp.89-97
• /
• 2017

Turbulent flow structure in the high amplitude meandering channel is complex due to secondary recirculation with helicoidal motions and shear layers formed by flow separation from the curved sidewall. In this work, the secondary flow and the superelevation of the water surface produced in the high-amplitude Kinoshita channel are reproduced by the unsteady Reynolds-averaged Navier-Stokes (RANS) computations using the VOF technique for resolving the variation of water surface elevation and three statistical turbulence models ($k-{\varepsilon}$, RNG $k-{\varepsilon}$, $k-{\omega}$ SST). The numerical results computed by a second-order accurate finite volume method are compared with an existing experimental measurement. Among applied turbulence models, $k-{\omega}$ SST model relatively well predicts overall distribution of the secondary recirculation in the Kinoshita channel, while all three models yield similar prediction of water superelevation transverse slope. The secondary recirculation driven by the radial acceleration in the upstream bend affects the flow structure in the downstream bend, which yields a pair of counter-rotating vortices at the bend apex. This complex flow pattern is reasonably well reproduced by the $k-{\omega}$ SST model. Both $k-{\varepsilon}$ based models fail to predict the clockwise-rotating vortex between a pair of counter-rotating vortices which was observed in the experiment. Regardless of applied turbulence models, the present computations using the VOF method appear to well reproduce the superelevation of water surface through the meandering channel.

• Journal of the Korean Society of Visualization
• /
• v.9 no.4
• /
• pp.63-68
• /
• 2011

The flow-field of a liquid-metal system is very important for the safety analysis and the design of the steam generator of liquid-metal fast breeder reactor. Dynamic neutron radiography (DNR) is suitable for a visualization and measurement of a liquid metal flow and a two-phase flow in a metallic duct. However, the three dimensional DNR techniques is not enough to obtain the velocity information in the wide channel up to now. In this research, a high speed DNR technique was applied to visualize the heavy liquid-metal flow field in the narrow channel with the HANARO-beam facility. The images were taken with a high frame-rate neutron radiography at 250 fps and analyzed with a Particle Image Velocimetry(PIV) method. The images were compared with the results of the commercial CFX code to study the feasibility of DNR technique for the measuring the heavy liquid-metal flow field. The PIV images could discern the turbulent vortex flow in the two-dimensional narrow channel.

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

This paper presents the results of a detailed experimental examination of fully developed asymmetric flows between annular tubes with square-ribbed surface roughness. The main emphasis of the research has been on establishing the turbulence structure, particularly in the central region of the channel where the two dissimilar wall flows interact. Measurements have included profiles of time mean velocities, turbulence intensities, turbulent shear stresses, triple velocity correlations, skewness, and flatness. The region of greatest interaction is characterized by strong diffusional transport of turbulent shear stress and kinetic energy from rough toward the smooth wall region, giving rise to an appreciable separation between the planes of zero shear stresses depending on positions of roughness on the walls.

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

The shelter effect of a porous wind fence on a triangular prism was experimentally investigated in a circulating water channel. A porous fence of porosity .epsilon.=38.5% was installed in front of the prism model. The fence and prism model were embedded in a turbulent boundary layer. The instantaneous velocity fields around the fence and prism model were measured by using the instantaneous velocity fields around the fence and prism model were measured by using the 2-frame PTV(Particle Tracking Velocimetry) system. By installing the fence in front of the prism, the recirculation flow region decreases compared with that of no fence case. The porous fence also decreases the mean velocity, turbulent intensity and turbulent kinetic energy around the prism. Especially, at the top of the prism, the turbulent kinetic energy is about half of that without the fence.

• Ocean Science Journal
• /
• v.43 no.2
• /
• pp.81-90
• /
• 2008

A study of estuarine flows during a neap tide was performed using 13-hour roving acoustic Doppler current profiles (ADCP) and conductivity-temperature-depth (CTD) profiles in the Altamaha River estuary, Georgia, U.S.A. The least-squared harmonic analysis method was used to fit the tidal ($M_2$) component and separate the flow into two components: the tidal and residual ($M_2$-removed) flows. We applied this method to depth-averaged data. Results show that the $M_2$ component demonstrates over 95% of the variability of observation data. As the flow was dominated by the $M_2$ tidal component in a narrow channel, the tidal ellipse distribution was essentially a back-and-forth motion. The amplitude of $M_2$ velocity component increased slightly from the river mouth (0.45 m/sec) to land (0.6 m/sec) and the phase showed fairly constant values in the center of the channel and rapidly decreasing values near the northern and southern shoaling areas. The residual flow and transport calculated from depth-averaged flow shows temporal variability over the tidal time scale. Strong landward flows appeared during slack waters which may be attributed to increased baroclinic forcing when turbulent mixing decreases.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.7
• /
• pp.1866-1872
• /
• 1994

Turbulent flow in a channel with a square rib periodically mounted on one wall is studied by large-eddy simulation(LES). An efficient 3D Navier-Stokes solver has been written for this geometry using a fractional step method and a multi-grid technique. The Reynolds number considered is 82, 000 based on the mean velocity above the obstacle height. Near-wall turbulence is approximated by a wall-layer model based on the turbulence intensity at the grid point nearest a solid wall. The results show a good qualitative agreement with experiments currently available for a single rib, indicating that LES can be a useful tool in simulating complex turbulent flows.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.7
• /
• pp.973-983
• /
• 2003

The suitability of high-order accurate, centered and upwind-biased compact difference schemes is evaluated for large eddy simulation of turbulent flow. Two turbulent flows are considered: turbulent channel flow at Re = 23000 and flow over a circular cylinder at Re = 3900. The effects of numerical dissipation on the finite differencing and aliasing errors and the subgrid-scale stress are investigated. It is shown through the simulations that compact upwind schemes are not suitable for LES, whereas the fourth order-compact centered scheme is a good candidate for LES provided that proper dealiasing of nonlinear terms is performed. The classical issue on the aliasing error and the treatment of nonlinear terms is revisited with compact difference schemes.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.8
• /
• pp.1046-1055
• /
• 1997

An anisotropic k-.epsilon. turbulence model for predicting the rotating flows is proposed with the simple inclusion of a new parameter dealing with the extra straining effects in the .epsilon.-equation. This model is employed to compute the effects of Coriolis forces on fully-developed flow in a rotating channel. The predicted results indicate that the present model captures fairly well the striking rotational-induced effects on the Reynolds stresses and the mean flow distributions, including the argumentation of turbulent transport on the unstable side (pressure surface) of the channel and its damping on the stable side (suction surface).

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.17 no.9
• /
• pp.799-807
• /
• 2005

The flow characteristics and the heat transfer rate on a surface by the interaction of a pair of vortices are studied numerically. To analyze the common flow up produced by vortex generators in a rectangular channel flow, the pseudo-compressibility viscous method is introduced into the Reynolds-averaged Navier-Stokes equation for 3-dimensional unsteady, incompressible viscous flows. To predict turbulence characteristics, a two-layer $k-\varepsilon$ turbulence model is used on the flat plate 3-dimensional turbulence boundary The computational results predict accurately Reynolds stress, turbulent kinetic energy and flow field generated by the vortex generators. The numerical results, such as thermal boundary layers, skin friction characteristics and heat transfers, are also reasonably close to the experimental data.

• 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.