• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.5
• /
• pp.349-357
• /
• 2009

Effect of computational domain size in simulating of periodic obstacle flow has been investigated for the flow past tube banks. Reynolds number, defined by freestream velocity ($U_{\infty}$) and cylinder diameter (d), was fixed as 200, and center-to-center distance (P) as 1.5d. In-line square array and staggered square array were considered. Drag coefficient, lift coefficient and Strouhal number were calculated depending on domain size. Circular cylinders were implemented on a Cartesian grid system by using an immersed boundary method. Boundary condition is periodic in both streamwise and lateral directions. Previous studies in literature often use a square domain with a side length of P, which contains only one cylinder. However, this study reveals that the domain size is improper. Especially, RMS values of flow-induced forces are most sensitive to the domain size.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.810-815
• /
• 2003

Flow-induced noise propagated from flow over a sphere is numerically investigated for laminar flow at Re = 300 and 425, and for turbulent flow at Re = 3700 and $10^4$, where the Reynolds number is based on the freestream velocity and the sphere diameter. The numerical method used for obtaining the flow over a sphere is based on an immersed boundary method in a cylindrical coordinate system. The Curle’s solutions of the Lighthill’s acoustic analogy with and without the far-field and compact-source approximation are used in order to investigate the noise field from flow over a sphere. Since the drag and lift forces change irregularly in time at Re = 425, 3700 and $10^{4}$, the noise propagates in a complicated manner. At Re = 300, 425 and $10^{4}$, the noise from dipole sources is much larger than that from quadrupole sources. On the other hand, at Re = 3700, the quadrupole source becomes dominant. The temporal variation of the flow-induced noise around a sphere is obtained at some observation points, which shows that the peak frequency corresponds to the Strouhal number associated with the wake instability.

• International Journal of Fluid Machinery and Systems
• /
• v.10 no.1
• /
• pp.76-85
• /
• 2017

Wall pressure measurements and flow visualization were conducted for a 90 degree elbow with an axis curvature radius the same as its inner diameter (125 mm). Reynolds numbers 320,000 and 500,000, based on the inner diameter and bulk velocity, were examined. A deflected inflow, having an almost constant velocity slope and a faster velocity at the inside, was introduced. Ensemble averaged pressure distributions showed that no difference of normalized pressure could be found in both the Reynolds number cases. Power spectral density functions of pressures exhibited that the fluctuation having the Strouhal number (based on the inner diameter and bulk velocity) of 0.6 existed in the downstream region of the elbow, which was 0.1 larger than that of the uniform inflow case [1]. Results of numerical calculations qualitatively coincided with the experimental ones.

• Structural Engineering and Mechanics
• /
• v.75 no.5
• /
• pp.621-631
• /
• 2020

A laboratory experiment is conducted is to investigate the behaviour of a low-mass-ratio and high aspect ratio flexible cylinder under vortex-induced vibration (VIV). A flexible cylinder with aspect ratio of 100 and mass ratio of 1.17 is towed horizontally to generate uniform flow profile. The range of Reynolds number is from 1380 to 13800. Vibration amplitude, in-line and cross-flow frequency response, amplitude trajectory, mean tension variation and hydrodynamic force coefficients are analyzed based on the measurement from strain gauges, load cell and CCD camera. Experimental results indicate that broad-banded lock-in region is found for the cylinder with a small Strouhal number. The frequency switches in the present study indicates the change of the VIV phenomenon. The hydrodynamic force responses provide more understanding on the VIV of a low mass ratio cylinder.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.5
• /
• pp.1023-1031
• /
• 1989

An experimental study is conducted to investigate the turbulent wake flow around tubes in vertical single row tube banks. All measurements are performed at Reynolds number(Re$_{max}$) 4.2*10$_{3}$ - 2.5*10$_{4}$ with varying tube spacings from the wide pitch ratio(H/D=3.07) to the very narrow one(H/D=1.23). Flow patterns are visualized using the smoke-wire method. Mean static pressures, velocity components, and various statistical quantities of turbulence are obtained by the computer on-line technique. In the case of wide tube spacings, the near wakes of tube show similar trends to those of a single tube, and their flow indicats an anisotropic turbulence. However, as the pitch ratio decreases, wide and narrow wakes appear alternately behind adjacent tubes due to the deflected flow. Also, in the case of H/D .leq. 1.54, Karman vortex is not formed at the side of relatively wide wake.e.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.13 no.2
• /
• pp.54-63
• /
• 2009

Large eddy simulation was conducted for flow development in a chamber with wall injection which simulates the cold flow in an idealized hybrid rocket motor. It was found that a peculiar timescale, roughly corresponding to St~0.5, resides in the flowfield resulting from the interaction between the main oxidizer and wall injected flows. However, the fact that this time characteristics is absent in the temperature field in the vicinity of the wall indicates that even a small regression rate renders the passive scalar, such as temperature, dissimilar to the velocity field. This implies that a classical approach, which assumes that constant turbulent Prandtl number, should be replaced by a more sophisticated turbulence models to accurately predict the temperature field in the hybrid motor.

• Wind and Structures
• /
• v.17 no.4
• /
• pp.379-397
• /
• 2013

Large Eddy Simulations (LES) were carried out to investigate the aerodynamic characteristics of a rectangular cylinder with side ratio B/D=5 at Reynolds number Re=22,000 (based on cylinder thickness). Particular attention was devoted to the effects of velocity shear in the oncoming flow. Time-averaged and unsteady flow patterns around the cylinder were studied to enhance understanding of the effects of velocity shear. The simulation results showed that the Strouhal number has no significant variation with oncoming velocity shear, while the peak fluctuation frequency of the drag coefficient becomes identical to that of the lift coefficient with increase in velocity shear. The intermittently-reattached flow that features the aerodynamics of the 5:1 rectangular cylinder in non-shear flow becomes more stably reattached on the high-velocity side, and more stably separated on the low-velocity side. Both the mean and fluctuating drag coefficients increase slightly with increase in velocity shear. The mean and fluctuating lift and moment coefficients increase almost linearly with velocity shear. Lift force acts from the high-velocity side to the low-velocity side, which is similar to that of a circular cylinder but opposite to that of a square cylinder under the same oncoming shear flow.

• Ocean Systems Engineering
• /
• v.4 no.3
• /
• pp.185-200
• /
• 2014

The objective of the present simulations is to evaluate the applicability of the standard $k-{\varepsilon}$ turbulence model in engineering practice in the subcritical to supercritical flow regimes. Two-dimensional numerical simulations of flow around a circular cylinder at $Re=1{\times}10^5$, $5{\times}10^5$ and $1{\times}10^6$, had been performed using Unsteady Reynolds-Averaged Navier Stokes (URANS) equations with the standard $k-{\varepsilon}$ turbulence model. Solution verification had been studied by evaluating grid and time step size convergence. For each Reynolds number, several meshes with different grid and time step size resolutions were chosen to calculate the hydrodynamic quantities such as the time-averaged drag coefficient, root-mean square value of lift coefficient, Strouhal number, the coefficient of pressure on the downstream point of the cylinder, the separation angle. By comparing the values of these quantities of adjacent grid or time step size resolutions, convergence study has been performed. Solution validation is obtained by comparing the converged results with published numerical and experimental data. The deviations of the values of present simulated quantities from those corresponding experimental data become smaller as Reynolds numbers increases from $1{\times}10^5$ to $1{\times}10^6$. This may show that the standard $k-{\varepsilon}$ model with enhanced wall treatment appears to be applicable for higher Reynolds number turbulence flow.

• Journal of the Korean Society of Combustion
• /
• v.17 no.1
• /
• pp.48-57
• /
• 2012

Laminar lifted propane coflow-jet flames diluted with nitrogen were experimentally investigated to determine heat-loss-related self-excitation regimes in the flame stability map and elucidate the individual flame characteristics. There exists a critical lift-off height over which flame-stabilizing effect becomes minor, thereby causing a normal heat-loss-induced self-excitation with O(0.01 Hz). Air-coflowing can suppress the normal heat-loss-induced self-excitation through increase of a Peclet number; meanwhile it can enhance the normal heat-lossinduced self-excitation through reducing fuel concentration gradient and thereby decreasing the reaction rate of trailing diffusion flame. Below the critical lift-off height. the effect of flame stabilization is superior, leading to a coflow-modulated heat-loss-induced self-excitation with O(0.001 Hz). Over the critical lift-off height, the effect of reducing fuel concentration gradient is pronounced, so that the normal heat-loss-induced self-excitation is restored. A newly found prompt self-excitation, observed prior to a heat-loss-induced flame blowout, is discussed. Heat-loss-related self-excitations, obtained laminar lifted propane coflow-jet flames diluted with nitrogen, were characterized by the functional dependency of Strouhal number on related parameters. The critical lift-off height was also reasonably characterized by Peclet number and fuel mole fraction.

• 유체기계공업학회:학술대회논문집
• /
• 2006.08a
• /
• pp.349-352
• /
• 2006

The conventional segregated finite element formulation produces a small and simple matrix at each step than in an integrated formulation. And the memory and cost requirements of computations are significantly reduced because the pressure equation for the mass conservation of the Navier-Stokes equations is constructed only once if the mesh is fixed. However, segregated finite element formulation solves Poisson equation of elliptic type so that it always needs a pressure boundary condition along a boundary even when physical information on pressure is not provided. On the other hand, the conventional integrated finite element formulation in which the governing equations are simultaneously treated has an advantage over a segregated formulation in the sense that it can give a more robust convergence behavior because all variables are implicitly combined. Further it needs a very small number of iterations to achieve convergence. However, the saddle-paint-type matrix (SPTM) in the integrated formulation is assembled and preconditioned every time step, so that it needs a large memory and computing time. Therefore, we newly proposed the P2PI semi-segregation formulation. In order to utilize the fact that the pressure equation is assembled and preconditioned only once in the segregated finite element formulation, a fixed symmetric SPTM has been obtained for the continuity constraint of the present semi-segregation finite element formulation. The momentum equation in the semi-segregation finite element formulation will be separated from the continuity equation so that the saddle-point-type matrix is assembled and preconditioned only once during the whole computation as long as the mesh does not change. For a comparison of the CPU time, accuracy and condition number between the two methods, they have been applied to the well-known benchmark problem. It is shown that the newly proposed semi-segregation finite element formulation performs better than the conventional integrated finite element formulation in terms of the computation time.