• 한국전산유체공학회:학술대회논문집
• /
• 2001.10a
• /
• pp.84-90
• /
• 2001

Secondary instability of flow past a circular cylinder is examined using direct numerical simulation at Reynolds number 220 and 250. The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. In x-y plane, the convection term is applied by the 5th order upwind scheme, and the pressure and viscosity terms are applied by the 4th order central difference. In spanwise, Navier-Stokes equation is distributed using Spectral Method. The critical Reynolds number for this instability is found to be about Re=190. The secondary instability leads re three-dimensionality with a spanwise wavelength about 4 cylinder diameters at onset (A-mode). Results of three-dimensional effect in wake of a circular cylinder are represented with spanwise and streamwise vorticity contours as Reynolds numbers.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1754-1759
• /
• 2003

Two-dimensional flow over a pair of circular cylinders in side-by-side arrangements at low Reynolds numbers has been numerically investigated in this study. Numerical simulations are performed, using the immersed boundary method, in the ranges of $40{\leq}Re{\leq}160$ and $g^{\ast}<5$, where Re and $g^{\ast}$ are, respectively, the Reynolds number and the spacing between the two cylinder surfaces divided by the cylinder diameter. Results show that total six kinds of wake patterns are observed over the ranges: antiphase-synchronized, inphase-synchronized, flip-flopping, single bluff-body, deflected, and steady wake patterns. It is found that the characteristics of the flow significantly depends both on the Reynolds number and gap spacing, with the latter much stronger than the former. Instantaneous flow fields, time traces, flow statistics and so on are presented to identify the wake patterns and then to understand the underlying mechanism. It is remarkable that, for the deflected wake pattern, the gap flow is deflected invariably to the cylinder of higher drag coefficient and the deflection way does not change at all. Moreover, the bifurcation phenomena where either of two wake patterns can occur are found at certain flow conditions.

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

Two-dimensional flow over a pair of circular cylinders in a side-by-side arrangement at low Reynolds numbers has been numerically investigated in this study Numerical simulations are performed, using the immersed boundary method, for the ranges of 40$\leq$Re$\leq$160 and $g^{*}$＜5, where Re and $g^{*}$ are, respectively, the Reynolds number and the spacing between the two cylinder surfaces divided by the cylinder diameter. Results show that a total of six kinds of wake patterns are observed over the ranges: antiphase-synchronized, inphase-synchronized, flip-flopping, single bluff-body, deflected, and steady wake patterns. It is found that the characteristics of the flow significantly depends both on the Reynolds number and gap spacing, with the latter much stronger than the former. Instantaneous flow fields, time traces, flow statistics and so on are presented to identify the wake patterns and then to understand the underlying mechanism. Moreover, the bifurcation phenomenon where either of two wake patterns can occur is found at certain flow conditions.ons.

• International Journal of Air-Conditioning and Refrigeration
• /
• v.10 no.1
• /
• pp.50-63
• /
• 2002

The Karman boundary-layer has been numerically investigated for the disturbance wave number, wave velocity, azimuth angle and radius (Reynolds number, Re). The disturbed flow over rotating disk can lead to transition at a much lower Re than that of the well-known Type I instability. This early transition is due to the excitation of the Type II. Presented are the neutral stability results concerning these instabilities by solving newly formulated stability equations with consideration of whole convective terms. When the present numerical results are compared with the previously known results, the value of critical Re corresponding to Type I is moved from ${Re}_{c.1}$=285.3 to 270.2 and the value corresponding to Type II from ${Re}_{c.2}$=69.4 to 36.9, respectively. Also, the corresponding wave number is moved fro)m $k_1$=0.378 to 0.386 for Type I; from $k_2$=0.279 to 0.385 for Type II. For Type II, the upped limit of wave number and azimuth angle is $k_u$=0.5872, $\varepsilon_u$=$-17.5^{\circ}$, while its lower limit is near $k_u$=0, $\varepsilon_u$=$-28.4^{\circ}$. This implies that the disturbances will be relatively fast amplified at small Re and within narrow bands of wave number compared with the previous results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.1
• /
• pp.118-125
• /
• 2004

The experimental study on flow and heat transfer characteristics was conducted to investigate and to compare the performance of swirl jet by a twisted tape as a swirl generator with the performance of impinging single circular jet in fully developed flow tube. The effects of jet Reynolds number(Re=8700, 13800, 20000, 26500), dimensionless distance of nozzle-to-plate(H/d=2, 4, 6, 8) and swirl ratio(S=0.11, 0.23, 0.30) of the jet on the local and average Nusselt number have been examined. Measurements of local heat transfer rate and flow patterns on the jet impinging plate were used naphthalene sublimation technique and flow visualization technique respectively. Mean velocity and turbulence intensity of the jet along the centerline were measured. With a twisted tape in the nozzle exit, average Nusselt number at the around area of stagnation point were higher than those without the twisted tape at H/d=2, 4 and with the increase of Reynolds number. With a twisted tape in the nozzle, in the case of H/d=2, Re=26500 and S=0.11, maximum local Nusselt number at the region of y/d=0 and x/d=0.44 was obtained.

• Journal of computational fluids engineering
• /
• v.17 no.3
• /
• pp.39-44
• /
• 2012

In this study, a longitudinal pitch effect on in-line tube bank heat transfer has been analyzed numerically. To verify the accuracy of the solver model and boundary conditions, global Nusselt number(Nu) and pressure drop across the 2 row tube bank are compared with the existing experimental correlations under 500 ~ 2,000 Reynolds number(Re) range. By changing transverse pitch($S_T$) or longitudinal pitch($S_L$) separately in tube bank, we're trying to identify the each effect on heat transfer. We found that the effect of transverse pitch can be accounted for Reynolds number evaluated with maximum velocity($V_{max}$) at the smallest flow area similar to most existing correlations. Variation of the longitudinal pitch($S_L$) has a greater impact on the heat transfer compared to the transverse pitch($S_T$). Overall Nusselt number increases with larger longitudinal pitch($S_L$), however individual Nusselt number of the tube row has significant difference after the first row.

• Journal of computational fluids engineering
• /
• v.17 no.3
• /
• pp.59-67
• /
• 2012

Large Eddy Simulation(LES) of turbulent mass transfer in fully developed turbulent pipe flow has been performed to study the effect of Reynolds number on the concentration fields at $Re_{\tau}=180$, 395, 590 based on friction velocity and pipe radius. Dynamic subgrid-scale models for the turbulent subgrid-scale stresses and mass fluxes were employed to close the governing equations. Fully developed turbulent pipe flows with constant mass flux imposed at the wall are studied for Sc=0.71. The mean concentration profiles and turbulent intensities obtained from the present LES are in good agreement with the previous numerical and experimental results currently available. The effects of Reynolds number on the turbulent mass transfer are identified in the higher-order statistics(Skewness and Flatness factor) and instantaneous concentration fields. The budgets of turbulent mass fluxes and concentration variance were computed and analyzed to elucidate the effect of Reynolds number on turbulent mass transfer. Furthermore, to understand the correlation between near-wall turbulence structure and concentration fluctuation, we present an octant analysis in the vicinity of the pipe wall.

• Journal of Ocean Engineering and Technology
• /
• v.37 no.4
• /
• pp.137-144
• /
• 2023

This study investigates the experimental parameters of particle image velocimetry (PIV) to enhance the measurement technique for turbulent flow fields around a circular cylinder at a Reynolds number (Re) of 1.4×10⁵. At the Korea Research Institute of Ships & Ocean Engineering (KRISO), we utilized the cavitation tunnel and PIV system to capture the instantaneous flow fields and statistically obtained the mean flow fields. An aspect ratio and blockage ratio of 16.7% and 6.0%, respectively, were considered to minimize the tunnel wall effect on the cylinder wakes. The optimal values of the pulse time and the number of flow fields were determined by comparing the contours of mean streamlines, velocities, Reynolds shear stresses, and turbulent kinetic energy under their different values to ensure accurate and converged results. Based on the findings, we recommend a pulse time of 45 ㎲ corresponding to a particle moving time of 3-4 pixels, and at least 3,000 instantaneous flow fields to accurately obtain the mean flow fields. The results of the present study agree well with those of previous studies that examined the end of the subcritical flow regime.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.4 no.2
• /
• pp.151-161
• /
• 2012

The present study numerically investigates the aerodynamic characteristics of two-dimensional wings in the vicinity of the ground by solving two-dimensional steady incompressible Navier-Stokes equations with the turbulence closure model of the realizable k-${\varepsilon}$ model. Numerical simulations are performed at a wide range of the normalized ground clearance by the chord length ($0.1{\leq}h/C{\leq}1.25$) for the angles of attack ($0^{\circ}{\leq}{\alpha}{\leq}10^{\circ}$) in the prestall regime at a Reynolds number (Re) of $2{\times}10^6$ based on free stream velocity $U_{\infty}$ and the chord length. As the physical model of this study, a cambered airfoil of NACA 4406 has been selected by a performance test for various airfoils. The maximum lift-to-drag ratio is achieved at ${\alpha}=4^{\circ}$ and h / C = 0.1. Under the conditions of ${\alpha}=4^{\circ}$ and h / C = 0.1, the effect of the Reynolds number on the aerodynamic characteristics of NACA 4406 is investigated in the range of $2{\times}10^5{\leq}Re{\leq}2{\times}10^9$. As Re increases, $C_l$ and $C_d$ augments and decreases, respectively, and the lift-to-drag ratio increases linearly.

• Journal of Power System Engineering
• /
• v.11 no.4
• /
• pp.56-64
• /
• 2007

A numerical simulation is carried out mixed convection in horizontal channel with a heat source from below of rectangular cavity. Finite volume method was employed for the discretization and PISO algorithm was used for calculating pressure term. The parameters governing the problem are the Reynolds number ($10^{-2}{\leq}Re{\leq}50$), the Rayleigh number ($10^3{\leq}Ra{\leq}2.06{\times}10^5$), the Prandtl number ($0.72{\leq}Pr{\leq}909$), the aspect ratio ($0.5{\leq}AR=W/H{\leq}2$) and the angle of inclination ($0^{\circ}{\theta}60^{\circ}$). Mean Nusselt number distributions were obtained and effect of Reynolds number, Rayleigh number and Prandtl number on mixed convection in the horizontal channel with rectangular cavity were investigated.