• Journal of the Society of Naval Architects of Korea
• /
• v.41 no.4
• /
• pp.30-37
• /
• 2004

Flow characteristics of the inflow ahead of a rotating propeller attached to a container ship model were investigated using a two-frame PIV (Particle Image Velocimetry) technique. Ensemble-averaged mean velocity fields were measured at four different blade phases. The mean velocity fields show the acceleration of inflow due to the rotating propeller and the velocity deficit in the near-wake region. The axial velocity distribution of inflow in the upper plane of propeller is quite different from that in the lower plane due to the thick hull boundary layer. The propeller inflow also shows asymmetric axial velocity distribution in the port and starboard side. As the inflow moves toward the propeller, the effect of phase angle variation of propeller blade on the inflow becomes dominant. In the upper plane above the propeller axis the inflow has very low axial velocity and large turbulent kinetic energy, compared with the lower plane. The boundary layer developed along the bottom surface of stern hull forms a strong shear layer affecting vortex structure of the propeller near-wake.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.1
• /
• pp.21-28
• /
• 2000

Although commercial PIV systems have been widely used for the non-intrusive velocity field measurement of fluid flows, they are still under development and have considerable room for improvement. In this study, a single-frame double-exposure PIV system using a high-resolution CCD camera was developed. A pulsed Nd:Yag laser and high-resolution CCD camera were synchronized by a home-made control circuit. In order to resolve the directional ambiguity problem encountered in the single-frame PIV technique, the second particle image was genuinely shifted in the CCD sensor array during the time interval dt. The velocity vector field was determined by calculating the displacement vector at each interrogation window using cross-correlation with 50% overlapping. In order to check the effect of spatial resolution of CCD camera on the accuracy of PIV velocity field measurement, the developed PIV system with three different resolution modes of the CCD camera (512 ${\times}$ 512, lK ${\times}$ IK, 2K ${\times}$ 2K) was applied to a turbulent flow which simulate the Zn plating process of a steel strip. The experimental model consists of a snout and a moving belt. Aluminum flakes about $1{\mu}m$ diameter were used as scattering particles for the liquid flow in the zinc pot and the gas flow above the zinc surface was seeded with atomized olive oil with an average diameter of 1-$3{\mu}m$. Velocity field measurements were carried out at the strip speed $V_s$=1.0 m/s. The 2K ${\times}$ 2K high-resolution PIV technique was significantly superior compared to the smaller pixel resolution PIV system. For the cases of 512 ${\times}$ 512 and 1K ${\times}$ 1K pixel resolution PIV system, it was difficult to get accurate flow structure of viscous flow near the wall and small vortex structure in the region of large velocity gradient.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.5
• /
• pp.1032-1043
• /
• 1989

The heat transfer characteristics of the drag reducing polymer solutions are investigated experimentally in the thermal entrance region of circular tube flows. Fluids used in experiments are the aqueous solutions of high molecular polymer, polyacrylamide Separan AP-273 and the range of polymer concentrations is from 20 to 1000 wppm. Two stainless steel tubes with inside diameter 8.5mm(L/D=712) and 10.3mm(L/D=1160) are used for the heat transfer flow loops. The flow loop is set up to measure friction factors and heat transfer coefficients of test sections in two different modes; the recirculating flow system and once-through flow system. The test tubes are heated directly by electricity to apply the constant heat flux boundary conditions to the wall. Three different types of adaptors are used to observe the effects of the upstream flow conditions of the heat transfer test sections. The viscosity and characteristic relaxation time of the test fluids circulating in the flow system are measured by the capillary tube viscometer and falling ball viscometer at regular time intervals. The installed adaptors exhibit slight effect on the entrance heat transfer of Newtonian fluid. However, no noticeable effects are observed for the entrance heat transfer of the drag reducing fluids. The order of magnitude of the thermal entrance lengths of the drag reducing fluids which follow the minimum friction asymptote is much longer than that of Newtonian fluids in turbulent flows. A new dimensionless parameter, the viscoelastic Graetz number, is defined and all the experimental data are recasted in terms of the viscoelastic Graetz number. The local Nusselt number of the viscoelastic fluids is represented as a function of flow behavior index n and the viscoelastic Graetz number. As degradation continues the viscosity and the characteristic relaxation time of the testing fluids decrease. Weissenberg number defined by the relaxation time and D/V appears to be a proper dimensionless parameter in describing degradation effects on heat transfer of the viscoelastic fluids.

• Fire Science and Engineering
• /
• v.11 no.1
• /
• pp.21-35
• /
• 1997

The natural convection and combined heat transfer induced by fire in a rectangular enclosure is numerically studied. The model for this numerical analysis is partially opened right wall. The solution procedure includes the standard k-$\varepsilon$ model for turbulent flow and the discrete ordinates method (DOM) is used for the calculation of radiative heat transfer equation. In numerical study, SIMPLE algorithm is applied for fluid flow analysis, and the investigations of combustion gas induced by fire is performed by FAST model of HAZARD I program. In this study, numerical simulation on the combined naturnal convection and radiation is carried out in a partial enclosure filled with absorbed-emitted gray media, but is not considered scattering problem. The streamlines, isothermal lines, average radiation intensity and kinetic energy are compared the results of pure convection with those of the combined convection-radiation, the combined heat transfer. Comparing the results of pure convection with those of the combined convection-radiation, the combined heat transfer analysis shows the stronger circulation than those of the pure convection. Three different locations of heat source are considered to observe the effect of heat source location on the heat transfer phenomena. As the results, the circulation and the heat transfer in the left region from heating block are much more influenced than those in the right region. It is also founded that the radiation effect cannot be neglected in analyzing the building in fire. And as the results of combustion gas analysis from FAST model, it is found that O2 concentration is decreased according to time. While CO and CO2 concentration are rapidly increased in the beginning(about 100sec), but slowly decreased from that time on.

• Journal of the Korean Society of Industry Convergence
• /
• v.25 no.1
• /
• pp.61-70
• /
• 2022

Fire damage time in high-rise buildings and wildland fire increasing every year. The use of high-pressure fire pumps is required to effectively extinguish fires. Reflecting the curvature effect of the fire hose occurring at the actual fire fighting site, this study provides a database of pressure drop, discharge velocity and maximum discharge height through C FD numerical analysis and it can provide using standards for fire extinguishing. Two Reynolds numbers of 200000 and 400000 were numerically analyzed at 0° -180° bending with water of 25℃ as a working fluid in hoses with a diameter of 65mm, a length of 15m, and a radius of curvature of 130mm. Realizable k-ε turbulence model was used and standard wall function was used. The pressure drop increases as the bending angle increases, and the maximum value at 90° and then decreases. The increasing rate is greater than the decrease. The velocity of the secondary flow also decreases after having the maximum value at 90°. The decreasing rate is greater than the increase. The turbulent kinetic energy increases to 120° and decreases with the maximum value. Pressure drop, velocity of the secondary flow, and turbulence kinetic energy are measured larger in the second bending region than in the first bending region.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.32 no.3
• /
• pp.180-201
• /
• 2020

It has been widely known that the effect of diffracted waves at the tip of composite breakwater with finite length causes the change of standing wave height along the length of breakwater, the spatial change of wave pressure on caisson, and the occurrence of meandering damage on the different sliding distance in sequence. It is hard to deal with the spatial change of wave force on trunk of breakwater through the two-dimensional experiment and/or numerical analysis. In this study, two and three-dimensional numerical techniques with olaFlow model are used to approach the spatial change of wave force including the impulsive breaking wave pressure applied to trunk of breakwater, the effect of rear region, and the occurrence of diffracted waves at the tip of caisson located on the high crested rubble mound. In addition, it is thoroughly studied the mean wave height, mean horizontal velocity, and mean turbulent kinetic energy through the numerical analysis. In conclusion, it is confirmed that the larger wave pressure occurs at the front wall of caisson around the still water level than the original design conditions when it generates the shock-crushing wave pressure checked by not two-dimensional analysis, but three-dimensional analysis through the change of wave pressure applied to the caisson along the length of breakwater.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.32 no.6
• /
• pp.531-552
• /
• 2020

In the previous study, both the wave characteristics at the tip of composite breakwater and on caisson were investigated by applying olaFlow numerical model of three-dimensional regular waves. In this paper, the same numerical model and layout/shape of composite breakwater as applied the previous study under the action of one directional irregular waves were used to analyze two and three-dimensional spatial change of wave force including the impulsive breaking wave pressure applied to trunk of breakwater, the effect of rear region, and the occurrence of diffracted waves at the tip of caisson located on the high crested rubble mound. In addition, the frequency spectrum, mean significant wave height, mean horizontal velocity, and mean turbulent kinetic energy through the numerical analysis were studied. In conclusion, the larger wave pressure occurs at the front wall of caisson around the still water level than the original design conditions when it generates the shock-crushing wave pressure in three-dimensional analysis condition. Which was not occurred by two-dimensional analysis. Furthermore, it was confirmed that the wave pressure distribution at the caisson changes along the length of breakwater when the same significant incident wave was applied to the caisson. Although there is difference in magnitude, but its variation shows the similar tendency with the case of previous study.

• Journal of Energy Engineering
• /
• v.2 no.3
• /
• pp.268-275
• /
• 1993

A series of parametric investigations are performed in order to resolve the flow characteristic of a Quebec city stoker incinerator. The parameters considered in this study are five internal configurations of the Quebec city stoker itself and its modified ones, primary air velocity, the injection velocity and angle of the secondary air, and the reduction of the stoker exit area. A control-volume based finite-difference method by Patankar together with the power-law scheme is employed for discretization. The resolution of the pressure-velocity coupling is made by the use of SIMPLEC algorithm. The standard, two equation, k-$\varepsilon$ model is incorporated for the closure of turbulence. The size of recirculation region, turbulent viscosity, the mass fraction of the secondary air and pressure drop are calculated in order to analyze the characteristics of flow field. The results are physically acceptable and discussed in detail. The flow field of the Quebec city stoker shows the strong recirculation zone together with the high turbulence intensity over the upper part of the incinerator.

• The KSFM Journal of Fluid Machinery
• /
• v.15 no.5
• /
• pp.60-66
• /
• 2012

One of commonly physical phenomena encountered in pump sump systems in which its significant influence to the hydraulic performance of pump system plays an important role in the field of fluid engineering, is the appearance of free surface and submerged vortices. In this paper, a study of the vortices behavior and their formative mechanism of asymmetry is considered in this paper by using numerical approach. The Reynolds-Averaged Navier-Stokes (RANS) equations and k-omega Shear Stress Transport turbulence model used to describe the properties of turbulent flows, in company with VOF multiphase model, are implemented by Fluent code with multi-block structured grid system. In the numerical simulation, the calculated elevation of air-water interface and vortex core contours are used to classify visually surface vortices as well as submerged vortices. It is shown that the free surface vortex is identified by the concavity of liquid region from the free surface and swirling flow at that own plane. To investigate the distinctive behavior of these vortices corresponding to each given flow rate at the same water level, some numerical testing of them are considered here in such a manner that the flow pattern of surface vortex are obtained similarly to the obtained results from experiment. Furthermore, the influence due to the change of grid refinement and the variation of depth of the concavity are also considered in this paper. From that, these influential factors will be implemented to design a good pump sump with higher performance in the future.

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.3
• /
• pp.48-56
• /
• 2002

The flow characteristics around KRISO 3600TEU container ship model have been experimentally investigated in a circulating water channel. The instantaneous velocity vectors were measured using 2-frame PIV measurement system. The mean velocity fields and turbulent statistics including turbulent kinetic energy and vorticity were obtained by ensemble-averaging 400 instantaneous velocity fields. The free stream velocity was fixed at 0.6m/s and the corresponding Reynolds number was $9{\times}10^5$. The test sections were divided into two regions, three transverse sections of the wake region(Station -0.5767, -1, -3) and five longitudinal sections of the wake((Z/(B/2)=0, 0.1, 0.2, 0.4, 0.6). In the wake region, large-scale longitudinal vortices of nearly same strength are symmetric with respect to the wake centerline and a relatively weak secondary vortex is formed near the waterline. With going downstream, the strength of longitudinal vortex is decreased and the wake region expands.