• Journal of computational fluids engineering
• /
• v.17 no.1
• /
• pp.44-53
• /
• 2012

Immersed boundary method(IBM) is a numerical scheme proposed to simulate flow field around complex objectives using simple Cartesian grid system. In the previous studies, the IBM has mostly been implemented to fractional step method based Navier-Stokes solvers. In this study, we implement the IBM to an incompressible Navier-Stokes solver which uses SIMPLE algorithm. The weight coefficients of the bi-linear and quadratic interpolation equations were formulated by using only geometric information of boundary to reconstruct velocities near IB. Flow around 2D circular cylinder at Re=40 and 100 was solved by using these formulations. It was found that the pressure buildup was not observed even when the bi-linear interpolation was adopted. The use of quadratic interpolation made the predicted aerodynamic forces in good agreement with those of previous studies. For an analysis of moving boundary, we smulated an oscillating circular cylinder with Re=100 and KC(Keulegan-Carpenter) number of 5. The predicted flow fields were compared with experimental data and they also showed good agreements.

• Journal of computational fluids engineering
• /
• v.17 no.1
• /
• pp.8-15
• /
• 2012

A code developed using the flux-difference splitting scheme on the hybrid Cartesian/immersed boundary method is applied to simulate three-dimensional internal waves. The material interface is regarded as a moving contact discontinuity and is captured on the basis of mass conservation without any additional treatment across the interface. Inviscid fluxes are estimated using the flux-difference splitting scheme for incompressible fluids of different density. The hybrid Cartesian/immersed boundary method is used to enforce the boundary condition for a moving three-dimensional body. Immersed boundary nodes are identified within an instantaneous fluid domain on the basis of edges crossing a boundary. The dependent variables are reconstructed at the immersed boundary nodes along local normal lines to provide the boundary condition for a discretized flow problem. The internal waves are simulated, which are generated by an pitching ellipsoid near an material interface. The effects of density ratio and location of the ellipsoid on internal waves are compared.

• Journal of computational fluids engineering
• /
• v.7 no.1
• /
• pp.10-19
• /
• 2002

The non-staggered(collocated) grid approach in which all the solution variables are located at the centers of control volumes is very popular for incompressible flow analyses because of its numerical efficiency on the curvilinear or unstructured grids. Rhie and Chow's paper is the first in using non-staggered grid method for SIMPLE algorithm, where pressure weighted interpolation was used to prevent decoupling of pressure and velocity. But it has been known that this non-staggered grid method has stability problems when pressure fields are nonlinear like in natural convection flows. Also Rhie-Chow scheme generates large numerical diffusion near curved walls. The cause of these unwanted problems is too large pressure damping term compared to the magnitude of face velocity. In this study the magnitude of pressure damping term of Rhie-Chow's method is limited to 1∼10% of face velocity to prevent physically unreasonable solutions. The wall pressure extrapolation which is necessary for cell-centered FVM is another source of numerical errors. Some methods are applied in a unstructured FV solver and analyzed in view of numerical accuracy. Here, two natural convection problems are solved to check the effect of the Rhie-Chow's method on numerical stability. And numerical diffusion from Rhie-Chow's method is studied by solving the inviscid flow around a circular cylinder.

• Journal of Mechanical Science and Technology
• /
• v.19 no.5
• /
• pp.1169-1181
• /
• 2005

Experimental and numerical studies on the unsteady wake field behind a square cylinder near a wall were conducted to find out how the vortex shedding mechanism is correlated with gap flow. The computations were performed by solving unsteady 2-D Incompressible Reynolds Averaged Navier-Stokes equations with a newly developed ${\epsilon}-SST$ turbulence model for more accurate prediction of large separated flows. Through spectral analysis and the smoke wire flow visualization, it was discovered that velocity profiles in a gap region have strong influences on the formation of vortex shedding behind a square cylinder near a wall. From these results, Strouhal number distributions could be found, where the transition region of the Strouhal number was at $G/D=0.5{\sim}0.7$ above the critical gap height. The primary and minor shedding frequencies measured in this region were affected by the interaction between the upper and the lower separated shear layer, and minor shedding frequency was due to the separation bubble on the wall. It was also observed that the position (y/G) and the magnitude of maximum average velocity $(u/u_{\infty})$ in the gap region affect the regular vortex shedding as the gap height increases.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.1
• /
• pp.262-273
• /
• 1991

In the present study, the steady, incompressible, isothermal, developing flow in a 90.deg. rectangular cross sectional strongly curved duct with aspect ratio 1:1.5 and Reynolds number of 9.4*10$^{4}$ has been investigated. Measurements of components of mean velocities, pressures, and corresponding components of the Reynolds stress tensor are obtained with a hot-wire anemometer and pitot tube. In general, flow in a curved duct is characterized by the secondary vortices which are driven mainly by centrifugal force-radial pressure gradient imbalance, and the stress field stabilizing effects near the convex wall and destablizing effects close to the concave wall. It was found that the secondary mean velocities attain values up to 39% of the bulk velocity and are largely responsible for the convections of Reynolds stress in the cross stream plane. Therefor upstream of the bend the Reynolds stress are low. Corresponding to the small boundary layer thickness. At successive planes, large values of Reynolds stress were observed near the concave surface and the side wall.

• Journal of computational fluids engineering
• /
• v.21 no.1
• /
• pp.110-116
• /
• 2016

The present study numerically investigated the deformation of the interface of two-phase fluids flow in a blast furnace. To simulate three-dimensional(3D) incompressible viscous two-phase flow in the furnace filled with the air and molten iron, the volume of fluid(VOF) method based on the finite volume method has been utilized. In addition, the porous medium with the porosity has been considered as the bed of the particles such as cokes and char etc. For the comparison, the single phase flow and the two-phase flow without the porosity have been simulated. The two-phase flow without porosity condition revealed the smooth parabolic profile of the free surface near the outlet. However, the free surface under the porosity condition formed the viscous finger when the free surface was close to the outlet. This viscous finger accelerated the velocity of the free surface falling and the outflow velocity of the fluids near the outlet.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.703-704
• /
• 2002

The two major problems related to the blood flow in a floating type polymer valve are thrombus formation and hemolysis. It is well known that the shear stress in the fluid and flow separation around the valve are blamed for such disastrous phenomena. In this viewpoint, through study of the flow field around the valve is imperative to improve design of the valve. The aim of this study is to investigate the fluid flow around a floating type polymer valve. The numerical method employed in this study is the finite element software called ADINA. Incompressible viscous flow is assumed for blood using the assumption of Newtonian fluid. In this study, two prominent features of the axisymmetric flow around the floating type polymer valve are observed: jet-like flows observed near the gap between the conduit and the valve, and recirculating flow downstream of the valve. We also provided a detailed description of shear stress field according to the variation of flow conditions. The shear stress in fluid has its maximum value near the gap between the valve and the conduit.

• Journal of computational fluids engineering
• /
• v.16 no.2
• /
• pp.49-55
• /
• 2011

Turbulent flows with wavy wall are simulated using Residual-based Variational Multiscale Method (RB-VMS) which is proposed by Bazilves et al(2007) as new Large Eddy Simulation methodology. Incompressible Navier-Stokes equations are integrated using Isogeometric analysis which adopt the basis function as NURBS. The Reynolds number is 6760 based on the bulk velocity and averaged channel height. And the amplitude (${\alpha}/{\lambda}$) of wavy wall is 0.05. The computational domain is $2{\lambda}{\times}1.05{\lambda}{\times}{\lambda}$ in the streamwise, wall normal and spanwise direction. Mean quantities and turbulent statistics near wavy wall are compared with DNS results of Cherukat et al.(1998). The predicted results show good agreement with reference data.

• Journal of Advanced Marine Engineering and Technology
• /
• v.24 no.1
• /
• pp.31-40
• /
• 2000

This paper describes the numerical computations of the interaction between the streamwise vortex and a flat plate 3-D turbulent boundary layer. In the present study, the main interest is in the behavior of the streamwise vortices introduced in turbulent boundary layers. The flow behind a vortex generator is modeled by the information that is avilable from studies on the dalta winglet. An algorithm of the solution of the incompressible Navier-Stokes equations for three-dimensional turbulent flows, together with a two layer turbulent model to resolve the near-wall flow, is based on the method of artificial compressibility. The present results show boundary layer distortion due to vortices, such as strong spanwise flow divergence and boundary thinning, and have a good agreement with the experimental data.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.195-199
• /
• 2005

Wave breaking phenomenon near the fore body of a ship is numerically simulated. The ship advance with uniform velocity in calm water. For the simulation, incompressible Navier-Stokes equations and continuity equation are adopted as governing equations. The simulation is carried out in staggered variable mesh system with finite difference method. Marker and Cell(MAC) method and Marker-Density method are employed to track the free surface. Body boundary conditions are satisfied with the adoption of porosity method and no-slip condition on the hull surface. The ship model has a wedge type fore-body, and the computational domain is an appropriate region around the fore-body. The computation results are compared with some experimental results. Also the difference of the free surface tracking methods are discussed.