• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.11a
• /
• pp.760-763
• /
• 2010

The expansion of 1D numerical code to 2D or 3D is needed in order to improve the analysis accuracy of the interior ballistics. The cut cell method has been imposed for the code expansion to multi dimensions. The MUSCL-Hancock scheme as a high resolution method has been selected. A feasibility of the cut cell method has been verified by analyzing the free piston problem.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.8 s.227
• /
• pp.953-960
• /
• 2004

In this study, the effects of freestream turbulence intensity on laminar-turbulent transition of separated shear layers in the wake of a circular cylinder are investigated using an immersed boundary method and LES. It is shown that the present numerical results without freestream turbulence for Re=3,900 based on bulk mean velocity and the cylinder diameter are in good agreement with other authors' experimental observations and numerical results, verifying our numerical methodology. Then a 'prescribed power spectrum' method is imposed to generate isotropic turbulence at the inlet of the computational domain at each time step. The principal effects of freestream turbulence intensity on flow statistics are investigated for Re=3,900. Statistical study reveals that the Reynolds stresses in the near-wake region gradually increase, and transition occurs further upstream, as the turbulence intensity increases. On the other hand, the bubble size behind the cylinder decreases as the turbulence intensity increases, which indicates that the freestream turbulence helps mean velocity be quickly recovered.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.12 s.255
• /
• pp.1218-1227
• /
• 2006

In this paper a parametric study using an immersed boundary method has been carried out to investigate the effects of stable density stratification on the wakes past two identical three-dimensional hills aligned in tandem. The Reynolds number based on the uniform inlet velocity and twice the hill height was fixed at Re=300 while the Froude number based on the inlet velocity and the hill height was retained at Fr=0.2. Neutral flow without density stratification was also computed for comparison. Under a strong stratification, vertical motion of fluid particles over the three-dimensional hills is suppressed and the wake structures behind the hills become planar. Depending on the distance between the two hills, the flow pattern of each wake is significantly affected by the stratification. There is a critical hill distance at which flow characteristics drastically change. Qualitative and quantitative features of the wake interaction are reported.

• Journal of computational fluids engineering
• /
• v.22 no.1
• /
• pp.81-87
• /
• 2017

The characteristics of flow past a square cylinder submerged under the free surface have been numerically studied. An immersed boundary method was adopted for implementation of the cylinder cross-section in a Cartesian grid system. Also, a level-set method was used to capture the interface of the two fluids. The case for Reynolds number 150 was examined. At the specific Reynolds number, by varying the gap ratio(0.25, 0.40, 0.55, 0.70, 1.00, 1.50, 2.50, 5.00) the effects of the free surface on the force coefficients and Strouhal number of vortex shedding were identified. The presence of the free surface very close to the cylinder significantly affects the shedding pattern, resulting in considerable deviation of the force coefficients and Strouhal number from those of the single-phase flow. In addition, the influence of Froude number was considered in this study. By increasing Froude number(0.2-0.4), flow topology change was identified at the specific gap ratios(0.40, 0.70, 1.50, 5.00).

• Proceedings of the KSME Conference
• /
• 2005.11a
• /
• pp.21.1-21.1
• /
• 2005

• Journal of the Korean Mathematical Society
• /
• v.49 no.5
• /
• pp.927-945
• /
• 2012

An approximation of singular source terms in elliptic problems is developed and analyzed. Under certain assumptions on the curve where the singular source is defined, the second order convergence in the maximum norm can be proved. Numerical results present its better performance compared to previously developed regularization techniques.

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.3
• /
• pp.18-27
• /
• 2002

In order to improve the accuracy of the solution near the boundary in an analysis of viscous flow around an arbitrary boundary which move and be deformed using an Eulerian concept, a level-set based grid deformation method is introduced to concentrate grid points near the boundary. This paper presents a new monitor function which can easily control the level of the concentration of grid points along the boundary. Computations for steady flow around a semi-circular cylinder mounted on the bottom of the flow domain were carried out to check the improvement of the solution using the adaptive grid system with an immersed boundary method. The present numerical results show a good agreement with the solutions obtained by a body fitted grid system and more accurate solutions than those computed with non-adaptive grid system. For the validation of mechanical usefulness of the present method, an expanded analysis of flow around multi-body fixed in the flow domain was carried out. Finally, the present moving adaptive grid method was applied to a two-dimensional bubble rise problem. The computed results show well adapted grid points around the boundary of the bubble at every time and a good agreement with the result calculated by fixed grid system.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.39 no.4
• /
• pp.325-333
• /
• 2015

This study was devoted to elucidating the change in the flow characteristics of a laminar flow past a circular cylinder by modifying the cylinder shape with O-grooves. A numerical analysis was performed in a two-dimensional framework. The cylinder was represented using an immersed boundary method and marker particles on a Cartesian grid system. The number and locations of the O-grooves were the key parameters. An analysis of the flow pattern and flow-induced forces was performed at Re = 40 and 50. In addition, we calculated the critical Reynolds number depending on the number of O-grooves and their locations.

• Journal of the Korean Society of Visualization
• /
• v.15 no.3
• /
• pp.52-56
• /
• 2017

In this study, we investigate the flow characteristics of lift-type disk which behaves the up-down motion using the large eddy simulation (LES) and immersed boundary method (IBM). Also, we perform the noise analysis using pressure field at 1.35 m distance and reveal the cause of noise to observe the vortical structure analysis of flow result. It is observed that vortical structure and wind shear were generated at leading edge and tower with high velocity deficit and flow separation. High magnitude of flow noise was observed in low frequency range which is from 30 Hz to 60 Hz. It was observed that vortical structure at leading edge was generated in frequency range from 33.3 Hz to 41.6 Hz. Temporal characteristic in vortical structure at leading edge was similar to noise characteristics, having the similar frequency ranges.

• Journal of Ocean Engineering and Technology
• /
• v.24 no.5
• /
• pp.8-15
• /
• 2010

The flow around a two-dimensional, rectangular cylinder that is freely falling in a channel was simulated using the immersed boundary method with direct forcing to determine the interactions between the fluid and the structure. The results of the present study were in good agreement with previous experimental results. Regardless of the H/L ratio (where H and L are the height and width of the rectangular cylinder, respectively), the flow structures had essentially the same pattern as the two symmetrical circulations that form about the horizontal center of the cylinder, with those centers located at each lateral position near the wake. When the cylinder approaches very close to the bottom, a jet-like flow appeared between the bottom of the rectangular cylinder and the channel. When the jet-like flow goes through the channel, surrounding fluids are sucked into this jet, forming the secondary vortices.