• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.39 no.5
• /
• pp.415-423
• /
• 2015

In this study, we present a grid refinement model in the lattice Boltzmann method (LBM) for two-dimensional incompressible fluid flow. That is, the model combines the desirable features of the lattice Boltzmann method and stream function-vorticity formulations. In order to obtain an accurate result, very fine grid (or lattice) is required near the solid boundary. Therefore, the grid refinement model is used in the lattice Boltzmann method for stream function-vorticity formulation. This approach is more efficient in that it can obtain the same accurate solution as that in single-block approach even if few lattices are used for computation. In order to validate the grid refinement approach for the stream function-vorticity formulation, the numerical simulations of lid-driven cavity flows were performed and good results were obtained.

• Journal of the Korean Society of Visualization
• /
• v.20 no.2
• /
• pp.45-54
• /
• 2022

We analyzed the low wall-shear stress area in the intracranial aneurysm that occurred at an anterior communicating artery with a special emphasis on vortical structures close to the wall. We reconstructed the aneurysm model from patient CTA data. We assumed blood as an incompressible Newtonian fluid and treated the blood vessel as a solid wall. The pulsatile boundary condition was applied at the inlet of the anterior cerebral artery. From the instantaneous flow field, we computed the histogram of the wall-shear stress over the aneurysm wall and found the low wall-shear stress event (< 0.4 Pa). This extreme event was due to the low wall-shear stress area that occurred at the daughter sac. We found that the merging of two vortices induced the low wall-shear stress area; one arises from the morphological characteristics of the daughter sac, and the other is formed by a jet flow into the aneurysm sac. The latter approaches the daughter sac, which ultimately leads to the strong ejection event near the daughter sac.

• Journal of the Korean Society of Visualization
• /
• v.20 no.1
• /
• pp.52-63
• /
• 2022

Direct numerical simulation of blood flow in a stenosed, patient-specific carotid artery was conducted to explore the transient behavior of blood flow with special emphasis on the wall-shear stress distribution over the transition region. We assumed the blood as an incompressible Newtonian fluid, and the vessel was treated as a solid wall. The pulsatile boundary condition was applied at the inlet of the carotid. The Reynolds number is 884 based on the inlet diameter, and the maximum flow rate and the corresponding Womersley number is approximately 5.9. We found the transitional behavior during the acceleration and deceleration phases. In order to quantitatively examine the wall-shear stress distribution over the transition region, the probability density function of the wall-shear stress was computed. It showed that the negative wall-shear stress events frequently occur near peak systole. In addition, the oscillatory shear stress index was used to further analyze the relationship with the negative wall-shear stress appearing in the systolic phase.

• Journal of the Society of Naval Architects of Korea
• /
• v.32 no.1
• /
• pp.83-93
• /
• 1995

A numerical and experimental investigation on the flow characteristics in the rectangular duct of an MHD propulsion system has been carried out. In numerical analysis, three-dimensional, steady-state, viscous, incompressible electrically conducting fluid flow under the influence of uniformly applied magnetic and electric fields was treated using a finite-difference technique. It was found from the numerical study that when the Lorentz force is weak, the typical parabolic velocity profile under a laminar flow condition changes to an M shaped profile near the electrode region and that the pressure increases linearly from the inlet toward the outlet of the MHD duct under constant electro-magnetic field. In experiment, thrust of the MHD propulsion system can be controlled easily by varying electrode current. The measured pressure gradient along the MHD duct is proportional to the Lorentz force, which is in agreement with the numerical results.

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

A level-set method is used for analyzing the behaviors of gas bubbles in two fluids incompressible viscous flow domain. The governing equations are solved by using a finite volume method. The numerical results are verified by comparing with the experimental and other computational results. Computations for the deformations and motions of one or multi-bubbles in the flow domain with the initial undisturbed free interface are conducted. It can be seen that numerical results for different surface tension and density ratio arise very different behaviors of bubbles. When bubbles rise near the free interface, the free interface gives some great influence on the behaviors of bubbles. The present results computed by a level-set method give useful information about the properties of bubble motions and deformations.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.4
• /
• pp.265-273
• /
• 2014

Wave-body interaction is simulated using a developed code based on the flux-difference splitting scheme for immiscible and incompressible fluids and the hybrid Cartesian/immersed boundary method. A free surface is captured as a moving contact discontinuity within a fluid domain and an approximated Riemann solver is used to estimate the inviscid flux across the discontinuity. Immersed boundary nodes are identified inside an instantaneous fluid domain near a moving body, then dependent variables are reconstructed at those immersed boundary nodes based on interpolation along local normal lines to the boundary. Free surface flows around an oscillating cylinder are simulated and the computed wave elevations are compared with other reported results. The generation of a solitary wave by a moving wave-maker is simulated and the time histories of wave elevations at two different points are compared with other results. The developed code is applied to simulate body motion of an elastically mounted circular cylinder as a solitary wave passes the body. The force acting on an elastically mounted cylinder is compared with the force acting on a fixed cylinder. Grid independency of the computed body motion is established based on a comparison of results using three different-size grids.

• Journal of ILASS-Korea
• /
• v.15 no.2
• /
• pp.74-82
• /
• 2010

In order to understand the behavior of transient gaseous injection used in an LPG (Liquefied Petroleum Gas) engine, turbulent incompressible transient jets with butane and propane were measured and analyzed at pressures of 1.5 bar and 2.0 bar with injector diameters of 3 mm and 5 mm. Mie-scattering method with a tracer was used, and images were processed to investigate the behavior of butane and propane jets. Distances from the nozzle to transition region were measured as $L_e/d_{inj}$=4.35~19.4, where $L_e$ and $d_{inj}$ indicate respectively a distance from nozzle to transition point and nozzle diameter. Slits and tubes around jet at near-field were introduced to measure the effect of entrainment and the diameter of jet, which revealed that the entrainment of surrounding air is significant for developing jet diameter. When the entrainment is restricted, the behavior of jet became deviating from the baseline. It was found that the virtual origin located outside of a nozzle towards jet tip within the conditions of this work, and its location was estimated as $x_o/d_{inj}$=0.56~7.25, where $x_o$ is a distance from nozzle to virtual origin.

• Journal of computational fluids engineering
• /
• v.17 no.3
• /
• pp.68-74
• /
• 2012

Near-wall effect on wakes behind particles is one of the important factors in precise tracking of particles in turbulent flows. However, most aerodynamic force models for particles did not fully consider the wall effect. In the present study, we focused on changes of hydrodynamic forces acting on a particle depending on wall proximity. To this end, we developed an immersed boundary method with multi-direct forcing incorporated to a fully implicit decoupling procedure for incompressible flows. We validate the present immersed boundary method through two-dimensional simulations of flow over a circular cylinder. Comprehensive parametric studies on the effect of the wall proximity on the drag and lift forces acting on an immersed circular cylinder in a channel flow are performed in order to investigate general flow patterns behind the circular cylinder for a wide range of Reynolds number (0.01 ${\leq}$ Re ${\leq}$ 200). As the cylinder is closer to the wall, the drag coefficient decreases while the lift coefficient increases with a local maximum. Maximum drag and lift coefficients for different wall proximities decrease with increment of Reynolds number. Normalized drag and lift coefficients by their maximum values show universal correlations between the coefficients and wall proximity in a low Reynolds number regime (Re ${\leq}$ 1).

• Transactions of the Korean Society of Automotive Engineers
• /
• v.7 no.8
• /
• pp.91-98
• /
• 1999

This paper describes the air flow characteristics inside the throttle valve. Tow-dimensional steady incompressible Navier-Strokes equation are solved numerically with embedding the conceopt of the artificial compressibility and adopting the Baldwin-Lomax turbulence model. With varying the valve opening angles(the Reynolds number )such as 15$^{\circ}$(5000) , 45$^{\circ}$(3000) , 75$^{\circ}$(7000) and 90$^{\circ}$(10000), respectively. tow cases, with a valve shaft and without one, are analysed. The pressure loss between the entrance and exit is severe at 15$^{\circ}$, 100 times as larger as that of 90$^{\circ}$ case, which also depends much on the existece of the valve shaft. The counter rotating vortices are formed over the valve plate with the shaft at only 75$^{\circ}$. They are smally and very large scale in front and back of the valve shaft , respectively. The velocity profiles of 15$^{\circ}$ and 90$^{\circ}$ at the exit are almost symmetric to the horizontal center line, however, the symmetricity is no longer maintained at 45$^{\circ}$ and 75$^{\circ}$ , and in addition, the flow at 75$^{\circ}$ is enforced a lot below center line. The pressure distribution on the walls is largely changed near the valve shaft, and its magnitude becomes great as the valve angle decreases.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.5 no.3
• /
• pp.333-347
• /
• 2013

The present study numerically models the interaction between a regular wave and the roll motion of a rectangular floating structure. In order to simulate two-dimensional incompressible viscous two-phase flow in a numerical wave tank with the rectangular floating structure, the present study used the volume of fluid method based on the finite volume method. The sliding mesh technique is adopted to handle the motion of the rectangular floating structure induced by fluid-structure interaction. The effect of the wave period on the flow, roll motion and forces acting on the structure is examined by considering three different wave periods. The time variations of the wave height and the roll motion of the rectangular structure are in good agreement with experimental results for all wave periods. The present response amplitude operator is in good agreement with experimental results with the linear potential theory. The present numerical results effectively represent the entire process of vortex generation and evolution described by the experimental results. The longer wave period showed a different mechanism of the vortex evolution near each bottom corner of the structure compared to cases of shorter wave periods. In addition, the x-directional and z-directional forces acting on the structure are analyzed.