• Journal of the Society of Naval Architects of Korea
• /
• v.45 no.4
• /
• pp.379-388
• /
• 2008

A hybrid Cartesian/immersed boundary code is expanded to simulate flow field around a three-dimensional body which undergoes large dynamic deformation. Immersed boundary nodes are automatically distributed based on the edges crossing triangles on body boundary. Velocity vectors are reconstructed at those immersed boundary nodes along local normal lines to the boundary. The reconstruction of pressure is avoided using the hybrid staggered/non-staggered grid method. The developed code is validated through comparisons with other results on laminar flow over a sphere. The code is applied to simulate flow around a foil which is attached to a body of revolution and rotates under periodic deformation. The periodic variation of the tip vortex is observed and the effects of the deformation on hydrodynamic force acting on the body are investigated.

• Journal of computational fluids engineering
• /
• v.13 no.4
• /
• pp.39-44
• /
• 2008

Numerical simulations of laminar flow over a sphere are conducted to investigate the effect of the Reynolds number on the characteristics of vortex shedding. The Reynolds numbers considered are between 300 and 475, covering unsteady planar-symmetric and asymmetric flows. Results show that the unsteady planar-symmetric flow can be categorized into two different regimes: single-frequency regime and multiple-frequency regime. The former has a single frequency component due to periodic shedding of the vortices with the same strength in every shedding cycle, while the latter has multiple frequency components due to cycle-to-cycle variation in the strength of shed vortices with the shedding angle fixed. The multiple-frequency planar-symmetric flow, which is newly found in the present study, occurs at Re=330${\sim}$360 between the single-frequency planar-symmetric flow and the asymmetric flow. On the other hand, the asymmetric flow occurs at Re${\geq}$365, where the vortices shed from the sphere show variation both in strength and shedding angle unlike the planar-symmetric flow. Also, it is shown that the breaking of planar symmetry is closely related to the imbalance of vortical strength between a pair of streamwise vortices.

• Korean Chemical Engineering Research
• /
• v.53 no.5
• /
• pp.609-613
• /
• 2015

The onset of instability induced by impulsive spin-down of the rigid-body flow placed in the gap between two coaxial cylinders is analyzed by using the energy method. In the present stability analysis the growth rate of the kinetic energy of the base state and also that of disturbances are taken into consideration. In the present system the primary flow is a transient, laminar one. But for the Reynolds number equal or larger than a certain one, i.e. $Re{\geq}Re_G$ secondary motion sets in, starting at a certain time. For $Re{\geq}Re_G$ the dimensionless critical time to mark the onset of vortex instabilities, ${\tau}_c$, is here presented as a function of the Reynolds number Re and the radius ratio ${\eta}$. For the wide gap case of small ${\eta}$, the transient instability is possible in the range of $Re_G{\leq}Re{\leq}Re_S$. It is found that the predicted ${\tau}_c$-value is much smaller than experimental detection time of first observable secondary motion. It seems evident that small disturbances initiated at ${\tau}_c$ require some growth period until they are detected experimentally.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.560-565
• /
• 2001

The objective of the present study is to investigate the characteristics for flow and wall shear stress in the aneurysm which is a local dilatation of the blood vessel. The numerical simulation using the commercial software for the laminar and steady flow were carried out over the diameter ratios(ratio of maximum diameter of aneurysm to the diameter of blood vessel) ranging from 1.5 to 2.5 and Reynolds number ranging from 900 to 1800. It was shown that a recirculating vortex occupied the entire bulge with its core located closer to the distal end of the bulge and the strength of vortex increased with increase of the Reynolds number and diameter ratio. Especially, for the Reynolds number of 1800 and diameter ratio of 2.5, the very weak secondary recirculating flow was produced at the left upper of the aneurysm. The position of a maximum wall shear stress was the distal end of the aneurysm(z=18mm) regardless of the Reynolds number and diameter ratios. But the maximum values of the wall shear stress increased in proportion to the increase of Reynolds number and diameter ratio.

• Journal of computational fluids engineering
• /
• v.6 no.2
• /
• pp.22-31
• /
• 2001

The two-dimensional unsteady incompressible viscous flow behind rectangular bluff bodies between two parallel plates was analyzed numerically. The steady state flow and the vortex flow behind rectangular bluff bodies in the channel were investigated for two regimes i.e., the laminar(Re = 100, 300, 500) and the turbulent flows(Re = 10⁴∼10/sup 6/). The vortex shedding was generated by a physical disturbance(6%) numerically imposed at the rear of the bluff bodies for a short time. It was observed that the perturbed flow became periodic after a transient period. And in the case of unsteady inflow, the sinusoidal pulsatile flow was applied as the inlet condition in the turbulent flow of Reynolds number of 1.0×10/sup 5/. FLUENT code was employed to solve the problems. The power-law scheme was used to get stable linearized equations and the PISO algorithm was applied to finding the solution of them.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2009.05a
• /
• pp.516-523
• /
• 2009

The artificial compressibility (AC) method for the incompressible Navier-Stokes equations in the generalized curvilinear coordinates using the primitive form is implemented. The main advantage of the AC approach is that the resulting system of equations resembles the system of compressible N-S equations and can thus be integrated in time using standard, well-established time-marching methods. The errors, which are the odd-even oscillation, for pressure field in using the artificial compressibility can be eliminated by using the $4^{th}$ order artificial dissipation term which is explicitly included. Even though this paper focuses exclusively on 2D laminar flows to validate and assess the performance of this solver, this numerical method is general enough so that it can be readily extended to carry out 3D URANS simulation of engineering flows. This algorithm yields practically identical velocity profiles and primary vortex and secondary vortices that are in excellent overall agreement with the results of the vorticity-stream function formulation (Ghia et al., 1982). However, the grid resolution have to be required to be large enough to express the various vortices.

• Wind and Structures
• /
• v.30 no.5
• /
• pp.475-483
• /
• 2020

This work investigates Reynolds number Re (= 50 - 200) effects on the flows around a single cylinder and the two tandem (center-to-center spacing L^⁎= L/D = 4) cylinders, each of a diameter D. Vorticity structures, Strouhal numbers, and time-mean and fluctuating forces are presented and discussed. For the single cylinder, with increasing Re in the range examined, the vorticity magnitude, Strouhal number and fluctuating lift all monotonically rise but time-mean drag, vortex formation length, and lateral distance between the two rows of vortices all shrink. For the two tandem cylinders, the increase in Re leads to the formation of three distinct flows, namely reattachment flow (50 ≤ Re ≤ 75), transition flow (75 < Re < 100), and coshedding flow (100 ≤ Re ≤ 200). The reattachment flow at Re = 50 is steady. When Re is increased from 75 to 200, the Strouhal number of the two cylinders, jumping from 0.113 to 0.15 in the transition flow regime, swells to 0.188. The two-cylinder flow is more sensitive to Re than the single cylinder flow. Fluctuating lift is greater for the downstream cylinder than the upstream cylinder while time-mean drag is higher for the upstream cylinder than for the other. The time-mean drags of the upstream cylinder and single cylinder behaves similar to each other, both declining with increasing Re.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.2
• /
• pp.127-134
• /
• 2004

Control of drag force on a circular cylinder using multiple detached splitter plates is numerically studied for laminar flow Two splitter plates with the same length as the cylinder diameter (d) are placed horizontally in the upstream of the cylinder and in the near-wake region, respectively. Their positions are described by the gap ratios (G$_1$/d, G$_2$/d), where G$_1$ represents the gap between the cylinder stagnation point and the rear edge of the upstream splitter plate, and G$_2$ represents the gap between the cylinder base point and the leading edge of the rear splitter plate. The drag varies with the two gap ratios; it has the minimum value at a certain set of gap ratios for each Reynolds number The upstream splitter plate decreases the stagnation pressure, while the rear splitter plate increases the base pressure by suppressing vortex shedding. This combined effect causes a significant drag reduction on the cylinder Particularly, the drag sharply increases past an optimum G$_2$/d; this seems to be related to a sudden change in bubble size in the wake region.

• 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.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.253-258
• /
• 2001

Control of drag force on a circular cylinder using a detached splitter plate is numerically studied for laminar flow. A splitter plate with the same length as the cylinder diameter(d) is placed horizontally in the wake region. Its position is described by the gap ratio(G/d), where G represents the gap between the cylinder base point and the leading edge of the plate. The drag varies with the gap ratio; it has the minimum value at a certain gap ratio for each Reynolds number. The drag sharply increases past the optimum gap ratio; this seems to be related to the sudden change in the bubble size in the wake region. This trend is consistent with the experimental observation currently available in case of turbulent flow. It is also found that the net drag coefficient significantly depends on the variation of base suction coefficient.