• Journal of the Society of Naval Architects of Korea
• /
• v.54 no.5
• /
• pp.406-414
• /
• 2017

An experimental study in a recirculating water channel was carried out to investigate the effect of large coherent structures to the skin friction on a flat plate. Particle Image Velocimetry (PIV) technique was used to quantify characteristic features of coherent structures growing to the boundary layer. In the PIV measurement, it is difficult to calculate the friction velocity near the wall region due to laser deflection and uncertainty so that Clauser fitting method at the logarithmic region was adopted to compute the friction velocity and compared with the one directly measured by the dynamometer. With changing the free-stream velocity from 0.5 m/s to 1.0 m/s, the activity of coherent structures in the logarithmic region was increased over three times in terms of Reynolds stress. The flow field was separated by Variable Interval Time Averaging (VITA) technique into the weak and the strong structure case depending on the existence large coherent structures in order to validate its effectiveness. The stream-wise velocity fluctuation was scanned through at the boundary thickness whether it had a large deviation from background flow. With coherent structures connected from near-wall to the boundary layer, mean wall shear stress was higher than that of weak structure case. Proper Orthogonal Decomposition (POD) analysis was also applied to compare the energy budget between them at each free-stream velocity.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.9 no.2
• /
• pp.97-104
• /
• 2005

A simple but efficient grid generation technique by using the modified compressible form of stream function has been formulated. Transformation of a physical plane to a streamline plane, the Von Mises Transformation, has been widely used to solve the differential equations governing flow phenomena, however, limitation arises in low velocity region of boundary layer, mixing layer and wake region where the relatively large grid spacing is inevitable. Modified Von Mises Transformation with simple mathematical adjustment for the stream function is suggested and applied to solve the confined coaxial turbulent jet mixing with simple $\kappa-\epsilon$ turbulence model. Comparison with several experimental data of axial mean velocity, turbulent kinetic energy, and Reynolds shear stress distribution shows quite good agreement in the mixing layer except in the centerline where the turbulent kinetic energy distributions were somewhat under estimated. This formulation is strongly suggested to be utilized specially for free turbulent mixing layers in axisymmetric flow conditions such as the investigation of mixing behavior, jet noise production and reduction for Turbofan engines.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.5 s.236
• /
• pp.577-589
• /
• 2005

The present study has numerically investigated two-dimensional laminar flow over a steadily rotating circular cylinder with a uniform planar shear, where the free-stream velocity varies linearly across the cylinder. It aims to find the combined effect of rotation and shear on the flow. Numerical simulations using the immersed boundary method are performed for the ranges of $-2.5{\le}\alpha{\le}2.5$ and $0{\le}K{\le}0.2$ at a fixed Reynolds number of Re=100, where a and K are respectively the dimensionless rotational speed and velocity gradient. Results show that the positive shear, with the upper side having the higher free-stream velocity than the lower one, favors the effect of the counter-clockwise rotation $(\alpha<0)$ but countervails that of the clockwise rotation $(\alpha>0)$. Accordingly, the absolute critical rotational speed, below which vortex shedding occurs, decreases with increasing K for $(\alpha>0)$, but increases for $\alpha>0$. The vortex shedding frequency increases with increasing \alpha (including the negative) and the variation becomes steeper with increasing K. The mean lift slightly decreases with increasing K regardless of the rotational direction. However, the mean drag and the amplitudes of the lift- and drag-fluctuations strongly depend on the direction. They all decrease with increasing K for $\alpha>0$, but increase for $\alpha<0$. Flow statistics as well as instantaneous flow folds are presented to identify the characteristics of the flow and then to understand the underlying mechanism.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.3
• /
• pp.725-733
• /
• 1990

An analysis of laminar mixed convection flow adjacent to the inclined flat surface which is subjected to a uniform temperature in a uniform free stream is performed. Nonsimilar boundary layed equation are derived by using the mixed convection parameters such that smooth transition from the purely forced convection limit to the purely free convection limit is possible. The governing equations are solved by a finite difference method using the coupled box scheme of sixth order. Numerical results are presented for prandtl numbers of 0.7 and 7 with the angle of inclination ranging from 0 to 90 degree from the vertical. The velocity distributions for the buoyancy assisting flow exhibit a significant overshoot above the free stream value in the region of intense mixed convection and the velocity field is found to be more sensitive to the buoyancy effect than the temperature field. The separation point near the wall was obtained for the buoyancy opposing flow. The local Nusselt number increases for buoyancy assisting flow and decreases for opposing flow with increasing value of the local Grashoff number in the mixed convection parameter. For large Prandtl number, the Nusselt number and the friction factor decrease significantly near the separation point. Present numerical predictions are in good agreement with recent experimental results by Ramachandran.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.3
• /
• pp.305-314
• /
• 2001

A modified model is proposed for calculation of transitional boundary layer flows. In order to develop the eddy viscosity model for the problem, the flow is divided into three regions; namely, pre-transition region, transition region and fully turbulent region. The pre-transition eddy-viscosity is formulated by extending the mixing length concept. In the transition region, the eddy-viscosity model employs two length scales, i.e., pre-transition length scale and turbulent length scale pertaining to the regions upstream and the downstream, respectively, and a universal model of stream-wise intermittency variation is used as a function bridging the pre-transition region and the fully turbulent region. The proposed model is applied to calculate three benchmark cases of the transitional boundary layer flows with different free-stream turbulent intensity (1%∼6%) under zero-pressure gradient. It was found that the profiles of mean velocity and turbulent intensity, local maximum of velocity fluctuations, their locations as well as the stream-wise variation of integral properties such as skin friction, shape factor and maximum velocity fluctuations are very satisfactorily predicted throughout the flow regions.

• 한국신재생에너지학회:학술대회논문집
• /
• 2009.11a
• /
• pp.478-481
• /
• 2009

As a first step for aerodynamic analysis of vertical axis wind turbine, dynamic stall characteristics of airfoil was investigated. Dynamic stall of wind turbine airfoil is caused by severe variation of angle of attack and relative velocity of flow around airfoil. Angle of attack and relative velocity can be expressed with tip speed ratio. Variation of angle of attack is strongly dependent on the tip speed ratio. For tip speed ratio, 1.4 and free stream velocity, 15m/s, dynamic stall characteristics of wind turbine airfoil is compared with those of oscillating airfoil having same angle of attack variation.

• Journal of the Society of Naval Architects of Korea
• /
• v.53 no.3
• /
• pp.210-216
• /
• 2016

Direct numerical simulations of a spatially developing turbulent boundary layer on a flat plate have been performed to verify the applicability of OpenFOAM and adapted mesh with prism layers to turbulent numerical simulation with high fidelity as well as provide a guideline on numerical schemes and parameters of OpenFOAM. Reynolds number based on a momentum thickness at inlet and a free-stream velocity was Re_θ=300. Time dependent inflow fields with near-wall turbulent structures were generated by a method of Lund et al. (1998), which was to extract instantaneous velocity fields from an auxiliary simulation with rescaled and recycled velocities at inlet. To ascertain the statistical characteristics of turbulent boundary layer, the mean profiles of streamwise velocity and turbulent intensities obtained from structured and adapted meshes were compared with the previous data.

• Journal of Korea Water Resources Association
• /
• v.56 no.spc1
• /
• pp.1037-1048
• /
• 2023

Fishways are constructed to ensure the fish migration because river-crossing structures such as dams and weirs cut off the stream longitudinal connectivity and influence on aquatic ecosystems. However, the passage efficiency of fishes varies depending on flow characteristics in the fishway and fish species. In this study, three-dimensional numerical simulations are carried out using a RANS model and the volume of fluid method for resolving free surface fluctuations to calculate the turbulent flow in the pool and weir type fishway. The Flow velocity and turbulent kinetic energy in the pool of fishway are analyzed according to variation of the upstream water level and the length of pool. The present numerical simulations reasonably well reproduce the stream flow and plunging flow characteristics in the pool. The simulation results show that the stream flow changes to the plunging flow as the length of the pool increases. When the upstream level increases, the stream flow becomes more evident. Key parameters related to the fish migration within the fishway such as the flow velocity and the turbulent kinetic energy are examined to assess the ascending possibility of fishes.

• Journal of Mechanical Science and Technology
• /
• v.16 no.5
• /
• pp.720-731
• /
• 2002

Flow with moving free surfaces is analyzed with an the Eulerian coordinate system. This study proposes a semi-implicit filling algorithm using VOF in which the PLIC (Piecewise Linear Interface Calculation) -type interface reconstruction method and the donor-acceptor-type front advancing scheme are adopted. Also, a new scheme using extrapolation of the stream function is proposed to find the velocity of the node that newly enters the computational domain. The effect of wall boundary conditions on the flow field and temperature field is examined by numerically solving a two-dimensional casting process.

• Journal of Ocean Engineering and Technology
• /
• v.23 no.5
• /
• pp.1-8
• /
• 2009

This paper is a continuation of the recent development on Hermite-based divergence free element method and deals with a non-isothermal fluid flow thru the buoyancy driven flow in a square enclosure with temperature difference across the two sides. The basis functions for the velocity field consist of the Hermite function and its curl while the basis functions for the temperature field consists of the Hermite function and its gradients. Hence, the number of degrees of freedom at a node becomes 6, which are the stream function, two velocities, the temperature and its x and y derivatives. This paper presents numerical results for Ra = 105, and compares with those from a stabilized finite element method developed by Illinca et al. (2000). The comparison has been done on 32 by 32 uniform elements and the degree of approximation of elements used for the stabilized finite element are linear (Deg. 1) and quadratic (Deg. 2). The numerical results from both methods show well agreements with those of De vahl Davi (1983).