• 한국전산유체공학회:학술대회논문집
• /
• 1996.05a
• /
• pp.68-73
• /
• 1996

The full incompressible Navier-Stokes equations are numerically integrated to solve the unsteady channel flow using a new numerical scheme of second-order accuracy developed by the authors. It is well known that in spite of the symmetry in the boundary condition and geometry, asymmetry can develop with time-dependency in a channel with sudden expansion. The instability of the shear flow and the cross-channel pressure contribute to such asymmetric flow. In this paper, we successfully generated a channel flow in which vortex waves were propagated downstream due to the harmonically oscillating inlet flow. The structure of the eddies and wall vorticity are parametrically investigated.

• Journal of Mechanical Science and Technology
• /
• v.20 no.12
• /
• pp.2250-2264
• /
• 2006

In micro- and nanoflows, the Boltzmann distribution is valid only when the electric double layers (EDL's) are not overlapped and the ionic distributions establish an equilibrium state. The present study has numerically investigated unsteady two-dimensional fully-developed electroosmotic flows between two parallel flat plates in the nonoverlapped and overlapped EDL cases, without any assumption of the Boltzmann distribution. For the study, two kinds of unsteady flows are considered: one is the impulsive application of a constant electric field and the other is the application of a sinusoidally oscillating electric field. For the numerical simulations, the ionic-species and electric-field equations as well as the continuity and momentum ones are solved. Numerical simulations are successful in accurately predicting unsteady electroosmotic flows and ionic distributions. Results show that the nonoverlapped and overlapped cases are totally different in their basic characteristics. This study would contribute to further understanding unsteady electroosmotic flows in micro- and nanofluidic devices.

• Journal of the Korean Society of Combustion
• /
• v.13 no.4
• /
• pp.8-15
• /
• 2008

When a circular cylinder is placed at the center of a slot burner nozzle, once stable Woflhard-Parker type laminar lean premixed flame is changed to an oscillating flame with self-induced noise. The wrinkled flame surface showed the same pattern and frequency of the Karman vortex street at the downstream of a circular cylinder. The interaction of flame with Karman vortex street is observed to be responsible for flame oscillation. The measured flame oscillation frequency is very similar to the estimated Karman vortex shedding frequency based on the St-Re relationship of the flow past circular cylinder, which could be considered as a strong evidence for the interaction between laminar pre-mixed flame and a Karman vortex street. As Reynolds number increases oscillation frequency decreases and the self-induced noise level increases as well as the flame front is more severly wrinkled. This result suggests that the flame/vortex interaction becomes more active at higher Re.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.358-361
• /
• 2010

A Reynolds-Averaged Navier-Stokes (RANS) code with transition prediction model is developed and the computational results on an oscillating airfoil are compared with the experimental data for OA209 airfoil. An approximated eN method that can predict transition onset points and the length of transition region is directly applied to the RANS code. The hysteresis loop in dynamic stall is compared for the computational results using transition prediction and fully turbulent models with the experimental data. Results with transition prediction show more correlation with the experimental data than the fully turbulent computation.

• Journal of computational fluids engineering
• /
• v.19 no.1
• /
• pp.80-86
• /
• 2014

Numerical investigation on the dynamic stall over an oscillating airfoil is presented. A Reynolds-Averaged Navier-Stokes (RANS) equations are coupled with transition transport equations for the natural transition. Computational results considering the turbulent transition are compared with the fully turbulent computations and the experimental data. Results with transition prediction show closer correlation with the experimental data than those with the fully turbulent assumption, especially in the reattachment region.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.950-955
• /
• 2003

The fluid propulsion mechanism of two oscillating flat plates is studied numerically using a discrete vortex method. Presently, the flat plates are assumed to be rigid. To analyze the closely coupled aerodynamic interference between the flat plates, a core addition scheme and a vortex core model are combined together. A calculated wake pattern for a flat plate in heaving oscillation motion is compared with the flow visualization. The effect of wake shapes on the aerodynamic characteristics of the flat plate in pitching oscillation is investigated. The velocity profiles behind the flat plates in pitching oscillations are plotted to investigate the possible thrust generation mechanism.

• Journal of the Korean Society of Visualization
• /
• v.4 no.2
• /
• pp.51-58
• /
• 2006

A boundary layer visualization was carried out in order to investigate the influence of Reynolds number on an oscillating airfoil. An NACA 0012 airfoil is sinusoidally pitched at the quarter chord point with oscillation amplitude of ${\pm}6^{\circ}$. A smoke-wire technique was employed to visualize the boundary layer and the near-wake. The freestream velocities are 1.98, 2.83 and 4.03m/s and corresponding chord Reynolds numbers are $2.3{\times}10^4,\;3.3{\times}10^4$, and $4.8{\times}10^4$, respectively. As the reduced frequency of K=0.1 is fixed, the corresponding frequency of an airfoil was adjusted in each case. The results reveal that the point at which the shear stress in an unsteady boundary layer separation disappears does not correspond with the position of the breakdown of the boundary layer, and that the breakdown of the boundary layer occurs further downstream.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.17 no.4
• /
• pp.13-20
• /
• 1980

The integral equation method to solve the boundary-value problem of a 2-dimensional body oscillating in the presence of a free surface generally breaks down at and near irregular frequencies due to the hypothetical flow inside the body. In this paper singularity distributions were extended to an inner free surface to remove the irregular frequency as Ohmatsu's work in 1978, and the solution for the above problem was found by using stream function. For various bodies including Lewis form cylinders, the hydrodynamic forces were calculated numerically at various wave numbers. From the results we concluded that the irregular frequencies can be removed even for the Lewis form cylinder as Ohmatsu done for circular cylinders, and calculated hydrodynamic forces by the present method are little higher than those of Ohmatsu's when the singularities are put on the inner free surface of the body. We specially point out that the solution for heaving motion converges in an oscillatory manner but not for swaying and rolling motions.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.8 s.239
• /
• pp.895-902
• /
• 2005

Lock-on characteristics of flow around a circular cylinder oscillating rotationally with a relatively high forcing frequency have been investigated experimentally. Dominant governing parameters are Reynolds number (Re), angular amplitude of oscillation (${\theta}_A$), and frequency ratio $F_R=f_f/f_n,\;where\;f_f$ is a forcing frequency and $f_n$ is a natural frequency of vortex shedding. Experiments were carried out under the conditions of $Re=4.14{\times}10^3,\;{\pi}/90{\leq}{\theta_A}{\leq}{\pi}/3,\;and\;F_R=1.0$. The effect of this active flow control technique on the lock-on flow characteristics of the cylinder wake was evaluated with wake velocity measurements and spectral analysis of hot-wire signals. The rotational oscillation modifies the flow structure of near wake significantly. The lock-on phenomenon always occurs at $F_R=1.0$, regardless of the angular amplitude ${\theta}_A$. In addition, when the angular amplitude is less than a certain value, the lock-on characteristics appear only at $F_R=1.0$,. The range of lock-on phenomena expands and vortex formation length is decreased, as the angular amplitude increases. The rotational oscillation create a small-scale vortex structure in the region just near the cylinder surface. At ${\theta}_A=60^{\circ}$, the drag coefficient was reduced about $43.7\%$ at maximum.

• Interaction and multiscale mechanics
• /
• v.5 no.1
• /
• pp.27-40
• /
• 2012

The two-dimensional incompressible flow of a linear viscoelastic fluid we considered in this research has rapidly oscillating initial conditions which contain both the large scale and small scale information. In order to grasp this double-scale phenomenon of the complex flow, a multiscale analysis method is developed based on the mathematical homogenization theory. For the incompressible flow of a linear viscoelastic Maxwell fluid, a well-posed multiscale system, including averaged equations and cell problems, is derived by employing the appropriate multiple scale asymptotic expansions to approximate the velocity, pressure and stress fields. And then, this multiscale system is solved numerically using the pseudospectral algorithm based on a time-splitting semi-implicit influence matrix method. The comparisons between the multiscale solutions and the direct numerical simulations demonstrate that the multiscale model not only captures large scale features accurately, but also reflects kinetic interactions between the large and small scale of the incompressible flow of a linear viscoelastic fluid.