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

• Geomechanics and Engineering
• /
• v.15 no.6
• /
• pp.1161-1171
• /
• 2018

The interaction between particles and fluid was investigated by IB-SEM numerical method which is a combination of combing the spectral/hp element method and the rigid immersed boundary method. The accuracy of this numerical method was verified based on the computed results with the traditional body-fitted mesh in numerical simulation of the flow through the cylinder. Then the governing equations of particles motion and contact in fluid are constructed. The movement of the particles and the interaction between the fluid and the particles are investigated. This method avoided the problem of low computational efficiency and error caused by the re-division of the grid when the solids moved. Finally, the movement simulation of multi particles in the fluid was carried out, which can provide a completely new numerical simulation method.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.131-136
• /
• 2008

Three-dimensional characteristics of fluid flow and heat transfer around a wavy circular cylinder having sinusoidal variation in cross sectional area along the spanwise direction are numerically investigated using the immersed boundary method. The three different wavelengths of ${\pi}/4$, ${\pi}/3$ and ${\pi}/2$ and at the fixed wavy amplitude of 0.1 have been considered to investigate the effects of waviness on especially the forced convection heat transfer around a wavy cylinder when the Reynolds and Prandtl numbers are 300 and 0.71, respectively. The present computational results for a wavy cylinder are compared with those for a smooth cylinder. The time- and total surface-averaged Nusselt number for a wavy cylinder with is larger than that for a smooth cylinder, whereas that with ${\lambda}={\pi}/4$ and ${\pi}/3$ is smaller than that for a smooth cylinder. However, because the surface area exposed to heat transfer for a wavy cylinder is larger than that for a smooth cylinder, the total heat transfer rate for a wavy cylinder with different wavelengths of ${\lambda}={\pi}/4$, ${\pi}/3$ and ${\pi}/2$ is larger than that for a smooth cylinder.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.37 no.2
• /
• pp.111-117
• /
• 2013

A valveless pump consisting of a pumping chamber with an elastic tube was simulated using an immersed boundary method. The interaction between the motion of the elastic tube and the pumping chamber generated a net flow toward the outlet through a full cycle of the pump. The net flow rate of the valveless pump was examined by varying the stretching coefficient, bending coefficient, and aspect ratio of the elastic tube. Photographs of the fluid velocity vectors and the wave motions of the elastic tube were examined over one cycle of the pump to gain a better understanding of the mechanism underlying the valveless pump. The relationship between the gap in the elastic tube and the average flow rate of the pump was analyzed.

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

• Bulletin of the Korean Mathematical Society
• /
• v.52 no.6
• /
• pp.1973-2000
• /
• 2015

In this study, we present a multidimensional open system for valveless pumping (VP). This system consists of an elastic tube connected to two open tanks filled with a fluid under gravity. The two-dimensional elastic tube model is constructed based on the immersed boundary method, and the tank model is governed by a system of ordinary differential equations based on the work-energy principle. The flows into and out of the elastic tube are modeled in terms of the source/sink patches inside the tube. The fluid dynamics of this system is generated by the periodic compress-and-release action applied to an asymmetric region of the elastic tube. We have developed an algorithm to couple these partial differential equations and ordinary differential equations using the pressure-flow relationship and the linearity of the discretized Navier-Stokes equations. We have observed the most important feature of VP, namely, the existence of a unidirectional net flow in the system. Our computations are focused on the factors that strongly influence the occurrence of unidirectional flows, for example, the frequency, compression duration, and location of pumping. Based on these investigations, some case studies are performed to observe the details of the ow features.

• Journal of computational fluids engineering
• /
• v.18 no.4
• /
• pp.17-24
• /
• 2013

In this study, we consider the characteristics of channel flow in the presence of an infinite streamwise array of equispaced identical rotating circular cylinders. This flow configuration can be regarded as a model representing a micro channel or an internal heat exchanger with cylindrical vortex generators. A numerical parametric study has been carried out by varying Reynolds number based on the bulk mean velocity and the cylinder diameter, and the gap between the cylinders and the channel wall for some selected angular speeds. An immersed boundary method was employed to facilitate implementing the cylinders on a Cartesian grid system. No-slip condition is employed at all solid boundaries including the cylinders, and the flow is assumed to be periodic in the streamwise direction. The presence of the rotating circular cylinders arranged periodically in the streamwise direction causes a significant topological change of the flow, leading to increase of mean friction on the channel walls. More quantitative results as well as qualitative physical explanations are presented to justify the effectiveness of rotating cylinders to modify flow topology, which might be used to enhance heat transfer on the channel walls.

• Journal of the Society of Naval Architects of Korea
• /
• v.48 no.3
• /
• pp.260-266
• /
• 2011

For long slender offshore structures, such as cables and pipe lines, their interaction with surrounding fluid flow becomes an important issue for global design of ocean systems. We employ a long circular cylinder as a representative case of slender offshore structure. A flexibly mounted cylinder in cross-flow generates complex vortex shedding and results in oscillation of the structure. In this paper, flow behind a circular cylinder at Re=100 is simulated. The vortex shedding pattern and flow induced motion are examined in the cross flow configuration as well as with various yaw-angled configurations. The "Lock-in" phenomenon is also observed when reduced velocity is approximately 4.0. The MAC Grid system, which is the typical grid system for Cartesian mesh and pressure correction methods, are used for solving the incompressible Navier-Stokes equations. Predictor/Corrector method is applied for obtaining a non-linear response of structure at the flexibly mounted. The existance and motion of the body is represented by the immersed boundary technique.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.22 no.11
• /
• pp.727-734
• /
• 2010

The present study numerically investigates the motion of a solid body suspended in the square enclosure with natural convection. A two-dimensional circular cylinder levitated thermally has been simulated by using thermal lattice Boltzmann method(TLBM) with the direct-forcing immersed boundary method. To deal with the ascending, falling or levitation of a circular cylinder in natural convection, the immersed boundary method is expanded and coupled with the TLBM. The circular cylinder is located at the bottom of a square enclosure with no restriction on the motion and freely migrates due to the Boussinesq approximation which is employed for the coupling between the flow and temperature fields. For different density ratio between the cylinder and the fluid, the motion characteristics of the circular cylinder for various Grashof numbers have been carried out. The Prandtl number is fixed as 0.7.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.2
• /
• pp.107-115
• /
• 2009

A Computational study was carried out in order to investigate the aerodynamic characteristics of circular cylinder moving near the wavy wall at a low Reynolds number of 50. Lattice Boltzmann method was used to simulate the flow field and immersed boundary method was combined to represent the moving cylinder and wavy wall regardless of the constructed grid in the domain. The aerodynamics characteristics of the cylinder moving near the wavy wall were represented by the comparing the lifting coefficients with various altitudes (H/D) and wave length and amplitudes of wavy wall. It indicated that the twice of increasing-decreasing variations of lifting coefficient are obtained while the cylinder moves near the wavy wall. The first variation is obtained where the cylinder locates near the peak of the wavy wall. Another variation occurs when the distance to the wavy wall becomes longer after passing the peak. It was also classified that three different patterns of relation between the lifting and drag coefficient of the cylinder. However, the classification is limited to the case of the same order of altitude, amplitude and wave length of the wavy wall.