• Journal of the Korean Society for Marine Environment & Energy
• /
• v.10 no.2
• /
• pp.107-111
• /
• 2007

A two-dimensional numerical method for inviscid two-fluid flows with evolution of density interface is developed, and an initially stationary two-dimensional fluid sheet surrounded by another fluid is studied. The interface between two fluids is modeled as a vortex sheet, and the flow field with the evolution of interface is solved by using vortex-in-cell/front-tracking method. The edge of the sheet is pulled back into the sheet due to surface tension and a blob is formed at the edge. This blob and fluid sheet are connected by a thin neck. In the inviscid limit, such process of the blob and neck formation is examined in detail and their kinematic characteristics are summarized with dimensionless parameters. The edge recedes at and the capillary wave propagating into the fluid sheet must be considered for better understanding of the edge receding.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.10
• /
• pp.757-766
• /
• 2009

In the present study, the effect of delta-wing vortex generators(DWVG) on the local heat transfer of the plate fin-oval tube was experimentally analyzed for Reynolds numbers for 2000, 2500 and 3200. The local heat transfer coefficient of the fin surface for four type DWVGs was measured by the naphthalene sublimation technique. As the results, the distribution of the heat transfer coefficient at rear of DWVGs showed longitudinal contours for common flow down DWVGs and wavy contours for common flow up DWVGs. The distribution showed many cell type contours at near wall and downstream for all DWVGs. Compared to case without DWVGs in present experimental tests, all DWVGs showed the best enhancement of heat transfer at Re=2000. Of 4 cases of DWVGs, D type showed the best enhancement of heat transfer.

• International Journal of Aeronautical and Space Sciences
• /
• v.10 no.2
• /
• pp.95-105
• /
• 2009

The flow fields around the HARTII rotor were numerically investigated using a viscous flow solver on adaptive unstructured meshes. An overset mesh and a deforming mesh technique were used to handle the blade motion including blade deflection, which was obtain from the HARTII experimental data. A solution-adaptive mesh refinement technique was also used to capture the rotor wake effectively. Comparison of the sectional normal force and pitching moment at 87% radial station between the two cases, with and without the blade deflection, showed that the blade loading is significantly affected by blade torsion. It was found that as the mesh was refined, the strength of tip vortex is better preserved, and the magnitude of high frequency blade loading, caused by blade-vortex interaction (BVI), is further magnified. It was also found that a proper time step size, which corresponds to the cell size, should be used to predict unsteady solutions accurately. In general, the numerical results in terms of the unsteady blade loading and the rotor wake show good agreement with the experimental data.

• Journal of computational fluids engineering
• /
• v.22 no.1
• /
• pp.51-58
• /
• 2017

A penalized VIC method offers an efficient hybrid particle-mesh algorithm to simulate an incompressible viscous flow passing a solid body in an infinite domain. In this manner, the computational domain needs to be restricted to a relatively small region to reduce computational cost which would be very high in case of using a large domain. In this paper, we present how to dispose of far-field particles to avoid an unnecessarily large computational domain. The present approach constraints expansion of the domain and thus prevents the incremental computational cost. To validate the numerical approach, a flow around an impulsively started sphere was simulated for Reynolds numbers of 100 and 1000.

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

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.407-413
• /
• 2011

VIC (Vortex-In-Cell) method for viscous incompressible flow is presented to simulate the wake behind a modified NACA16 foil. With uniform rectangular grid, the velocity in field is calculated using streamfunction from vorticity field by solving the Poisson equation in which FFT(Fast Fourier Transform) is combined with 2nd order finite difference scheme. Here, LES(Large Eddy Simulation) with Smagorinsky model is applied for turbulence calculation. Effective viscosity is formulated using magnitude of strain tensor(or vorticity). Then the turbulent diffusion as well as viscous diffusion becomes particle strength exchange(PSE) with averaged eddy viscosity. The well-established panel method is combined to obtain the irrotational velocity and to apply the no-penetration boundary condition on the body panel. And wall diffusion is used for no-slip condition numerical results of turbulent stresses are compared with experimental results (Bourgoyne, 2003). Before comparing process, LES(Large Eddy Simulation) SGS(Subgrid scale) stress is transformed Reynolds averaged stress (Winckelmans, 2001).

• International Journal of Fluid Machinery and Systems
• /
• v.3 no.3
• /
• pp.211-220
• /
• 2010

The objective of this study is to clarify the occurrence of the high-speed mode of unsteady swirling flows in straight tubes. The unsteady flows generated in the tube were measured by means of a semiconductor-type pressure transducer and an FFT analyzer. The high-speed mode measured has rotational speed which is approximately equal to or higher than the peripheral velocity of the swirling flow. The unsteady flow is due to cell rotation in the circumferential direction of the tube. The occurrence of the high-speed mode was confirmed, and the characteristics (rotational speed, pressure amplitude, and phase) of this mode were clarified. In order to understand the measured unsteady flows, the three dimensional vortex core profiles were discussed based on the distributions of the pressure amplitude and phase.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2011.11a
• /
• pp.649-652
• /
• 2011

A series of hybrid rocket combustion experiments were carried out with PMMA/GOx changing diameter and length of the disk installed at pre-chamber. The disk can generate vortex shedding flow and change flow conditions prior to entering the fuel grain which could also alter the combustion characteristics and pressure oscillations. Isolated dimple-like surface roughness patterns distributed all over the fuel surface, which can be thought of as a realization of the inherent flow instability. It is very likely that the formation of cell structures is originated from the modification of boundary layer characteristics of an entering oxidizer flow caused by a blowing effect mainly taking place near the wall. This coincided with our LES results. It would be a meaningful basis to understand combustion instability of hybrid rocket motor.

• Journal of Ship and Ocean Technology
• /
• v.10 no.3
• /
• pp.1-16
• /
• 2006

Despite the high cost of memory and CPU time required to resolve the boundary layer, a viscous unstructured grid solver has many advantages over a structured grid solver such as the convenience in automated grid generation and vortex capturing by solution adaption. In present study, an unstructured Cartesian grid solver is developed on the basis of the existing Euler solver, NASCART-GT. Instead of cut-cell approach, immersed boundary approach is applied with ghost cell boundary condition, which can be easily applied to a moving grid solver. The standard $k-{\varepsilon}$ model by Launder and Spalding is employed for the turbulence modeling, and a new wall function approach is devised for the unstructured Cartesian grid solver. Developed approach is validated and the efficiency of the developed boundary condition is tested in 2-D flow field around a flat plate, NACA0012 airfoil, and axisymmetric hemispheroid.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.7 s.94
• /
• pp.1805-1812
• /
• 1993

Two-dimensional spin-up in a rectangular domain is analysed by the numerical computation of the Navier-Stokes equations. The cells are in most cases generated by the vorticity developed near the uper and lower surfaces. Moreover, the movement and interaction of those vortices play a key role in establishing the quasi-steady state. The critical phenomena observed in the previous experiment turns out to be caused by the critical movement of the vortices.