• Journal of the Korean Society of Visualization
• /
• v.17 no.2
• /
• pp.10-16
• /
• 2019

Considering the role of viscosity in the hemorheology, the characteristics of non-Newtonian fluid are important in the pulsatile blood flows. Stenosis, with an abnormal narrowing of the vessel, contributes to block blood flows to downstream tissue and lead to plaque rupture. Therefore, systematic analysis of blood flow around stenosed vessels is crucial. In this study, non-Newtonian behaviors of blood analog fluids around the micro-stenosis with 60 % severity in diameter of $500{\mu}m$ was examined by using CFX under the pulsatile flow conditions with the period of 10 s. Viscosity information of two non-Newtonian fluids were obtained by fitting the value of normal blood and highly viscous blood. As the Newtonian fluid, the water at room temperature was used. During the pulsatile phase, wall shear stress (WSS) is highly oscillated. In addition, high viscous solution gives rise to increases the variation in the WSS around the micro-stenosis. Highly oscillating WSS enhance increasing tendency of plaque instability or rupture and damage of the tissue layer. These results, related to the influence on the damage to the endothelium or stenotic lesion, may help clinicians understand relevant mechanisms.

• Korea-Australia Rheology Journal
• /
• v.12 no.3_4
• /
• pp.165-173
• /
• 2000

Three-dimensional numerical simulation of isothermal/nonisothermal coextrusion process of two immiscible polymers through a rectangular channel has been done using the finite element method. The encapsulation phenomenon with the less viscous layer encapsulating the more viscous layer was investigated with the generalized Newtonian fluids. The interface position around the symmetric plane obtained by numerical simulation nearly coincided with the one observed in experiments, but the degree of encapsulation was less than the one observed experimentally. Open boundary condition method was found to be applied to the simulation of nonisothermal coextrusion process, however, the results are not far from those using the fully developed boundary condition, because the temperature development along the downstream direction is very slow in the case of convection dominated flow. When the inlet velocity is increased, the interface profile does not change in isothermal flow, while it moves upward in nonisothermal situation. The degree of encapsulation decreases along the downstream direction in nonisothermal flow. When the inlet temperature increases compared to the wall temperature, the outlet interface moves downward and the degree of encapsulation increases. The difference of degree of encapsulation between the simulation and the experiments seems to arise from the viscoelastic effect of the materials. It was concluded that the nonisothermal effect alone does not explain the complex coextrusion process and the viscoelastic effect needs to be considered.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.152-157
• /
• 1999

For the simulation of the blood flow in coronary artery, the system modeling of coronary hemodynamics is combined with CFD technique. The blood flow in coronary artery interacts with the global coronary circulation. Especially in case of the coronary artery with stenosis, the interaction plays an important role in the hemodynamics of the circulation. In this study we present a combined numerical approach using both the CFD technique for flow simulation and the global system model of coronary circulation. We use a lumped parameter model for the global simulation of coronary circulation whereas the finite element method is employed to compute the viscous flow field in stenosed coronary artery, The time variation of the pressure drop due to stenosis is obtained from the proposed numerical method. Numerical results shows that the flow resistance and pressure drop due to stenosis has a relatively large value in systole.

• Journal of the Semiconductor & Display Technology
• /
• v.17 no.4
• /
• pp.80-85
• /
• 2018

The present investigation on optical fiber mass manufacturing features the computational modeling and simulation on a double layer liquid coating process on glass fiber surface. The computational model employs a simplified geometry of typical fiber coating system which consists of primary and secondary coating dies along with secondary coating cup. The viscous dissipation in coating flow is incorporated into the double layer coating process simulations. Heavy temperature dependence of coating liquid viscosity is also considered in the model. The computational results found that the effects of viscous dissipation on both primary and secondary coating layer thicknesses are highly significant at higher drawing speed. Several important coating process parameters such as supply temperature and pressure of primary and secondary coating liquids are investigated and discussed in order to appreciate how those parameters affect the double layer coating layer thickness on fast moving glass fiber.

• 한국연소학회:학술대회논문집
• /
• 2002.11a
• /
• pp.49-56
• /
• 2002

A two-dimensional direct numerical simulation is performed to investigate the formation characteristics of a single vortex interacting with $CH_4/N_2$-Air counterflow nonpremixed flame. The numerical method was based on a predictor-corrector scheme for a low Mach number flow. The detailed transport properties and a 16-step augmented reduced mechanism are adopted in this calculation. The budgets of the vorticity transport equation arc examined to reveal the mechanisms leading to the formation, evolution and dissipation of a single vortex interacting with counterflow nonpremixed flame. It is found that the stretching term, which depends on the azimuthal component of vorticity, and radial velocity, mainly generates vortieitv in non-reacting and reacting flows. The viscous and baroclinic torque term destroy the vorticity in non-reacting flow. In addition, the baroclinic torque term due to density and pressure gradient generates vorticity, while viscous and the volumetric expansion terms due to density gradient destroy vorticity in reacting flow.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.4
• /
• pp.395-405
• /
• 2004

The filling pattern and an adaptive grid refinement based on the finite element method and Eulerian mesh advancement approach have been developed to analyze incompressible transient viscous flow with free surfaces. The governing equation for flow analysis is Navier-Stokes equation including inertia and gravity effects. The mixed FE formulation and predictor-corrector method are used effectively for unsteady numerical simulation. The flow front surface and the volume inflow rate are calculated using the filling pattern technique to select an adequate pattern among four filling patterns at each triangular control volume. By adaptive grid refinement, the new flow field that renders better prediction in flow surface shape is generated and the velocity field at the flow front part is calculated more exactly. In this domain the elements in the surface region are made finer than those in the remaining regions for more efficient computation. Using the proposed numerical technique, the collapse of a water dam has been analyzed to predict flow phenomenon of fluid and the predicted front positions with respect to time have been compared with the reported experimental results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.11
• /
• pp.1348-1358
• /
• 2004

The filling pattern and an adaptive grid refinement based on the finite element method and Eulerian mesh advancement approach have been developed to analyze incompressible transient viscous flow with free surfaces. The governing equation fur flow analysis is Navier-Stokes equation including inertia and gravity effects. The mixed FE formulation and predictor-corrector method are used effectively for unsteady numerical simulation. The flow front surface and the volume inflow rate are calculated using the filling pattern technique to select an adequate pattern among seven filling patterns at each tetrahedral control volume. By adaptive grid refinement, the new flow field that renders better prediction in flow surface shape is generated and the velocity field at the flow front part is calculated more exactly. In this domain the elements in the surface region are made finer than those in the remaining regions for more efficient computation. The collapse of a water dam and the filling of a fluidity spiral have been analyzed. The numerical results have been in good agreement with the experimental results and the efficiency of the adaptive grid refinement and filling pattern techniques have been verified.

• Ocean Systems Engineering
• /
• v.8 no.4
• /
• pp.361-380
• /
• 2018

The dynamic performance of the wave energy converter (WEC) rotor with different geometric parameters such as depth of submergence and beak angle has been assessed by considering the linear potential flow theory using WAMIT solver and along with the computational fluid dynamics (CFD). The effect of viscous damping is incorporated by conducting numerical free decay test using CFD. The hydrodynamic coefficients obtained from the WAMIT, viscous damping from the CFD and estimated PTO damping are used to solve the equation of motion to obtain the final pitch response, mean optimal power and capture width. The viscous damping is almost 0.9 to 4.6 times when compared to the actual damping. It is observed that by neglecting the viscous damping the pitch response and power are overestimated when compared to the without viscous damping. The performance of the pitch WEC rotor in the Jeju western coast at the Chagwido is analyzed using Joint North Sea Wave Project (JONSWAP) spectrum and square-root of average extracted power is obtained. The performance of WEC rotor with depth of submergence 2.8 m and beak angle $60^{\circ}$ found to be good compared to the other rotors.

• Journal of Navigation and Port Research
• /
• v.27 no.5
• /
• pp.473-478
• /
• 2003

Numerical study was made on the flow characteristics around a circular pipeline between parallel walls. The incompressible Navier-Stokes equations were solved by using a third-order upwind differential scheme. When the distance near a wall is small enough, the vortex shedding is almost completely suppressed because of the interaction with the wall boundary layer separation. This study aims to clarify the characteristics of the vortex shedding regime as the body approaches a wall as Reynolds number varies. The feature of separated vorticity dynamics is analyzed at different conditions with particular attention to the interaction between the pipeline wake and the induced separation on the plane walls.

• Nuclear Engineering and Technology
• /
• v.51 no.2
• /
• pp.573-578
• /
• 2019

Anisotropic distribution of the turbulent kinetic energy and the near-field excitations are the main causes of the steady state Flow-Induced Vibration (FIV) which could lead to fretting wear damage in vertically arranged supported slender rods. In this article, a combined Computational Fluid Dynamics (CFD) and Computational Structural Mechanic (CSM) approach named two-way Fluid-Structure Interaction (FSI) is used to investigate the modal characteristics of a typical rod's vibration. Performance of an Unsteady Reynolds-Average Navier-Stokes (URANS) and Large Eddy Simulation (LES) turbulence models on asymmetric fluctuations of the flow field are investigated. Using the LES turbulence model, any large deformation damps into a weak oscillation which remains in the system. However, it is challenging to use LES in two-way FSI problems from fluid domain discretization point of view which is investigated in this article as the innovation. It is concluded that the near-wall meshes whiten the viscous sub-layer is of great importance to estimate the Root Mean Square (RMS) of FIV amplitude correctly as a significant fretting wear parameter otherwise it merely computes the frequency of FIV.