• Journal of the Society of Naval Architects of Korea
• /
• v.32 no.3
• /
• pp.51-61
• /
• 1995

A numerical method was developed to solve the Navies-Stokes equations for unsteady laminar flow around a hydrofoil. The present method used the finite-difference scheme in the collocated grid system and the pressure-Poisson method was employed to obtain divergence-free velocity field each time step. The numerical method was applied at first to laminar flow around a circular cylinder to confirm capability of the code. In the next, calculations were carried out for a hydrofoil in an unbounded fluid at the Reynolds number of $10^4$ in order to investigate unsteady phenomena with vortex shedding. The calculate results showed reasonable features about laminar vortex shedding around a streamlined body.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.2
• /
• pp.139-152
• /
• 1998

In the present study, the flow characteristics of developing laminar oscillatory flows in a square -sectional 180 deg. curved duct are investigated experimentally. The experimental study using air in a square-sectional 180 deg. curved duct is carried out to measure velocity distributions with a data acquisition and LDV (Laser Doppler Velocimetry) processing system. In this system, Rotating Machinery Resolver (RMR) and PHASE program are used to obtain the results of unsteady flows. The major flow characteristics of developing oscillatory flows are found by analyzing velocity curves, mean velocity profiles, time-averaged velocity distribution of secondary flow, wall shear stress distributions, and entrance lengths. In a lower dimensionless angular frequency, the axial velocity distribution of laminar oscillatory flow in a curved duct shows a convex shape in a central part and axial symmetry. The maximum value of wall shear stress in a lower dimensionless angular frequency is located in an outside wall, but according to increasing the dimensionless angular frequency, the maximum of wall shear stress is moved to inner wall. The entrance lengths of laminar oscillatory flows in a square-sectional 180 deg. curved duct is obtained to 90 deg. of bended angle of duct in this experimental conditions.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2004.05a
• /
• pp.222-228
• /
• 2004

The effects of rotation on the unsteady laminar flow past a circular cylinder is numerically investigated in the present study. We obtained the numerical solutions for unsteady two-dimensional governing equation for the flow using two different numerical schemes. One is an accurate spectral method and another is finite volume method. Above all, the flow around a stationary circular cylinder is investigated to understand the basic phenomenon of flow separation, bluff body wake. Also, the validation of our own codes, expecially based on FVM, is carried out by the comparison of results obtained from our simulations using two different schemes and previous numerical and experimental studies. By the effect of rotation, the mean lift increases and drag deceases, which well represent the previous study.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.3
• /
• pp.350-357
• /
• 1997

The transient incompressible flow behind the axisymmetric bluff body is numerically simulated using the random vortex method(RVM). Based on the vorticity formulation of the unsteady Navier-Stokes equations, the Lagrangian approach with a stochastic simulation of diffusion using random walk technique is employed to account for the transport processes of the vortex elements. The numerical solutions for 2-dimensional recirculating flow behind a backward-facing step in the laminar range of Reynolds number are compared with experimental data. The present simulation focuses on the transitional flow regime where the recirculation zone behind the bluff body becomes highly unsteady and large-scale vortex eddies are shed from the bluff body wake due to intrinsic shear layer instabilities. The unsteady vertical flow structures and the mixing characteristics behind the bluff body are discussed in detail.

• Journal of computational fluids engineering
• /
• v.7 no.3
• /
• pp.35-43
• /
• 2002

This paper addresses a numerical method for predicting transient temperature distributions in the wall of a curved pipe subjected to internal laminar thermally-stratified flow. A simple and convenient numerical method of treating the unsteady conjugate heat transfer in non-orthogonal coordinate systems is presented. Numerical calculations are performed for the transient evolution of thermal stratification in two curved pipes, where one has thick wall and the other has so thin wall that its presence can be negligible in the heat transfer analysis. The predicted results show that the thermally stratified flow and transient conjugate heat transfer in a curved pipe with a finite wall thickness can be satisfactorily analyzed by the present numerical method, and that the neglect of wall thickness in the prediction of pipe wall temperature distributions can provide unacceptably distorted results for the cases of pipes with thick wall such as safety related-piping systems of nuclear power plant.

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

A Navier-Stokes code with a modified low Reynolds number k-.epsilon. turbulence model was used to study the unsteady transitional boundary layer flow due to rotor-stator interaction. The modification, proposed by Launder, to improve prediction of stagnation flows was incorporated to the low Reynolds number k-.epsilon. turbulence model by Fan-Lakshminarayana-Barnett. Numerical solution is shown to capture well the calmed laminar flow as well as the wake induced transitional strip due to rotor-stator interaction and shows improvement, in terms of onset of transition and its length, over previous Euler/boundary layer solution. The turbulent kinetic energy shows local maximum along the upstream rotor wake in the wake induced transitional strip and this characteristics is observed untill the end of transition. The wake induced strip also shown apparent even in the laminar sublayer as the upstream rotor wake penetrates inside the boundary layer.

• Journal of Advanced Marine Engineering and Technology
• /
• v.15 no.2
• /
• pp.64-72
• /
• 1991

The unsteady, two-dimensional laminar flow and heat transfer of an incompressible, constant-property fluid flowing around a circular cylinder are numerically analyzed. The Navier-Strokes equation and the energy equation are solved by the finite difference method. The range of the Reynolds number is 10 to 100 and the Prandtl number considered is 0.7. The contours of the flow pattern, equi-vorticity line and isotherm pattern around a circular cylinder are shown. Also, numerical solutions of the surface vorticity, pressure coefficient, drag coefficient, local Nusselt number and mean Nusselt number are obtained. The numerical results for the final time fo calculation are compared with the other available experimental and theoretical results for the steady state and are found to be in good agreement with them.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.5
• /
• pp.652-659
• /
• 2000

The purpose of this study is to present the most effective numerical scheme to calculate the unsteady flows. In order to calculate the flow quantities of flow past a circular cylinder, Three-time level and five convective schemes are applied to unsteady and convective terms, respectively. The values obtained are compared with those from the existing experimental and numerical studies. At Reynolds numbers up to 160, time intervals can be expanded 10 times of Implicit Euler scheme using Three-time level method, and it is found that QUICK and CUI schemes work much stable than others even if less grid density conditions. The combination of Three-time level and QUICK scheme gives high resolutions for laminar unsteady problems with PC level.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2004.11a
• /
• pp.210-211
• /
• 2004

The effect of oscillating on the unsteady laminar flow past a circular cylinder is numerically investigated in the present study. Our study is to analyze the vortex formation behind a circular cylinder for different rotary oscillation conditions. And then we are study to portray the unsteady dynamics of wake flows. We decide lock-on region by observing the phase switching phenomena We classify the vortex formation patterns in the primary lock-on region The present study is to identify the quasi-periodic state around lock-on region. At the boundary between lock-on and non-lock-on the shedding frequency is bifurcated. After the bifurcation, one frequency follow the forcing frequency ($S_f$) and the other returns to the natural shedding frequency ($St_0$). In the quasi-periodic state, the variation of magnitudes and relevant phase changes of $C_L$ with forcing phase are examined.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.5 s.248
• /
• pp.480-488
• /
• 2006

Numerical simulations of flow over two spheres placed in a tandem arrangement are conducted to investigate the effect of the inter-sphere spacing on the flow characteristics. The Reynolds numbers considered are 100, 250, 300 and 425, corresponding to steady axisymmetric, steady planar-symmetric, unsteady planar-symmetric, and unsteady asymmetric flows, respectively, in the case of a single sphere. For small inter-sphere spacings, the flow past two spheres is more stable than that past a single sphere. For example, with the spacing of the sphere radius, the flow is steady axisymmetric up to Re=300. However, for relatively large spacings, the flow past two spheres becomes unstable and vortex shedding takes place even at Re=250. The drag coefficient of the rear sphere decreases significantly with decreasing inter-sphere spacing due to reduction of the stagnation pressure, thus being smaller than that of the front sphere. Also, the rear sphere shows large fluctuations of the lift force as compared to the front one in the case of unsteady flow.