• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.7
• /
• pp.920-928
• /
• 2003

The objective of the present study was to investigate the characteristics of flow and wall shear stress under steady and pulsatile flow in the aneurysm. The numerical simulation using the software were carried out for the diameter ratios ranging from 1.5 to 3.0, Reynolds number ranging from 900 to 1800 and Womersley number, 15.47. For steady flow, it was shown that a recirculating vortex occupied the entire bulge with its core located closer to the distal end of the bulge and the strength of vortex increased with increase of the Reynolds number and diameter ratio. The position of a maximum wall shear stress was the distal end of the aneurysm regardless of the Reynolds number and diameter ratios. For the pulsatile flow, a recirculating flow at the bulge was developed and disappeared for one period and the strength of vortex increased with the diameter ratio. The maximum values of the wall shear stress increased in proportion to the diameter ratio. However, the position of a maximum wall shear stress was the distal end of the aneurysm regardless of the diameter ratios.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.747-752
• /
• 2003

Many clinical studies have suggested that the blood flow in ideal geometry is involved in the development of atherosclerosis. This study simulated blood flow in the abdominal artery with real geometry to investigate MWSS(mean wall shear stress), AWSS(amplitude of wall shear stress) and OSI(oscillator shear index). The calculation grid for the real geometry was constructed by extracting the surface of arterial wall from CT(Computed Tomography) or MRI(Magnetic Resonance Imaging) sheets called as DICOM (Digital Imaging and Communications in Medicines). The calculated MWSS, AWSS and OSI are much different from those of ideal geometry calculation. The MWSS increased while the AWSS decreased. Many shear forces are related to shapes of gradient. This paper will give clinical datum where the MWSS, AWSS and OSI are strong or weak. The hemodynamic analysis based on real geometry can provide surgeons with more reliable information about the effect of blood flow.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.1
• /
• pp.39-44
• /
• 2002

The arterial stenosis related to the intimal thickening of the arterial wall is the main cause of many diseases in human cardiovascular system. Hemodynamic behavior of the blood flow is influenced by the presence of the arterial stenosis. In this study, effects of the pulsatile flow, caused by the periodic motion of the heart, on the blood flow and its interaction with the arterial stenosis are analyzed by the FEM-based computational fluid dynamics. As a result, it was found that the characteristics of the pulsatile flow in the artery with stenosis are quite different from those of the steady flow. And, the pulsatile flow condition affects the wall shear stress, which is one of the most important physiological parameters in the hemodynamics.

• 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 Korea Water Resources Association
• /
• v.49 no.12
• /
• pp.1035-1044
• /
• 2016

In joint portions of the levee and the barrier, complex 3-dimensional flow was generated and collapse of revetment occurred frequently. For these reasons, it is necessary to install the joint revetment with greater stability as compared with the general revetment at the joint portions. However, design criteria for joint revetment was not presented in River Design Criteria (KWRA, 2009). Therefore it is necessary to research for engineering design of the joint revetment. In this study, hydraulic experiments were performed under various flow conditions in order to realize the collapse conditions of riprap and carried out in 20.0 m straight open channel with one side levee and the width was 4.0 m. The diameter of riprap covered around joint revetment was 0.03 m and the inlet discharges were $0.5{\sim}0.8m^3/s$. The numerical simulations were performed under same conditions with experiment. as results of this numerical simulations, the influence range was confirmed from the distribution of flow characteristics and shear stress. As a result, the riprap diameter of the joint revetment was calculated from 4.1 to 6.9 times greater than that of general revetment. As the inlet discharge was large, the range of vulnerable area was developed long in the downstream direction despite of same withdrawal velocity of riprap. Through this study, the methods of calculating the riprap diameter and influence range were proposed according to hydraulic characteristics of flow around joint revetment. At a later study, if additional experiments about effect of flood plane and various types of barrier is applied, it is expected that rational design method with stability of joint revetment can be proposed.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.37 no.10
• /
• pp.895-900
• /
• 2013

A two-dimensional Stokes flow past a circular cylinder in a channel is analyzed. The circular cylinder is located at the center of the channel, and a plane Poiseuille flow exists upstream and downstream far from the circular cylinder. The Stokes approximation is used, and the flow is investigated analytically by using the eigenfunction expansion and the least square methods. From the analysis, the stream function and pressure distribution are obtained, and the pressure and shear stress distributions on the circular cylinder and channel wall are calculated. The additional pressure drop induced by the circular cylinder and the force exerted on it are calculated as functions of the length of the radius of the circular cylinder. For a typical length of the radius of the circular cylinder, the streamline pattern and pressure distribution are shown.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.6
• /
• pp.609-615
• /
• 2011

A two-dimensional Stokes flow past a vertical plate in a channel is analyzed. The vertical plate is located at the center of the channel, and plane Poiseuille flow exists far upstream and downstream of the vertical plate. The Stokes approximation is used, and the flow is investigated analytically using the method of eigenfunction expansion and the point collocation method. From the analysis, the stream function and pressure distribution are obtained, and the pressure and shear stress distributions on the plate and channel wall are calculated. The additional pressure drop induced by the vertical plate and the force exerted on it are calculated as functions of the length of the vertical plate. For a typical length of the vertical plate, the streamline pattern and pressure distribution are shown. In addition, numerical analysis of laminar flow with a small Reynolds number is carried out to analyze the effect of a small Reynolds number on the flow pattern.

• Journal of Biomedical Engineering Research
• /
• v.18 no.2
• /
• pp.187-196
• /
• 1997

The three-dimensional flow analysis using the finite volume method is presented to compare the steady flow characteristics of blood with those of blood substitutes such as water and aqueous polymer solution in an idealized double branching model. The model is used to simlllate the region of the abdominal aorta near the celiac and superior mesenteric branches. Apparent viscosities of blood and the aqueous Separan solution are represented as a function of shear rate by the Carreau model, Water and aqueoiu Separan AP-273 500wppm solution are frequently used as blood substitutes in vitro experiments. Water is a typical Newtonian fluid and blood and Separan solution are non-Newtonian fluids. Flow phenomena such as velocity distribution, pressure variation and wall shear stress distribution of water, blood and polymer solution are quite different due to differences of the rheological characteristics of fluids. Flow phenomena of polymer solution are qualitatively similar to those of blood but the phenomena of water are quite different from those of blood and polymer solution. It is recommended that a lion-Newtonian fluid which exhibits very similar rheological behavior to blood be used in vitro experiments. A non-Newtonian fluid whose rheological characteristics are very similar to those of blood should be used to obtain the meaninylll hemodynamic data for blood flow in vitro experiment and by numerical analysis

• The Korean Journal of Rheology
• /
• v.6 no.1
• /
• pp.30-40
• /
• 1994

고분자 물질의 신장점도를 측정하기 위하여 설계된 수렴관을 지나는 시험 유체 M1 에대하여 유한요소방법으로 수치모사를 수행하였다. 구성방정식은 세 개의 이완시간을 가진 적분형 K-BKZ모형을 사용하였다. 신장변형이 지배적이고 변형속도가 매우 큰 흐름에 대하 여 실험적 방법으로 측정이 가능한 범위까지 수치모사를 수행하였다. 두 개의 압력 측정꼬 지 사이의 벽면 압력차에 대하여 압력 신호로 측정한 실험값을 수치모사결과와 비교하였다. 걷보기 전단속도가 매우 큰 1300s-1에 이르는 높은 유속의 전 실험범위에 대하여 안정된 수 치해를 얻을수 있었다. 3$0^{\circ}C$에서는 모든 실험범위의 유속에서 압력차에 대한 수치모사 결과 가 실험값과 잘일치했다. 21$^{\circ}C$에서는 0.1$\times$10-3m3/s보다 낮은 유속범위에서 실험값과 일치하 는 결과를 얻었으나 그보다 높은 유속에서 실험값과 일치하는 결과를 얻었으나 그보다 높은 유속에서 실험값과 다른 경향의 결과를 얻었다. 이것은 낮은 온도 높은 유속 조건에서 M1 유체의 성질이 불안정하고 또한 그러한 조건의 실험에서 발생한 압력 측정꼭지 부근의 기포 들이 정확한 압력측정에 영향을 끼쳤기 때문이다. 수치모사 결과로부터 얻은 압력과 응력분 포로부터 수렴관 유변측정기의 유동특성을 밝힐수 있었다. 이는 실험적 방법을 통해서는 얻 기 어려운 결과들로서 중요한의미를 가진다. 특별한 모양을 갖도록 설계된 수렴관을 통과하 는 M1 유체가 중심부근에서 일정한 신장변형속도로 변형됨을 확인할 수 있었으며 수직응력 은 지수적으로 증가하다가 축소부분을 지난 후 매우 장점도를 얻기 위하여 신장변형속도가 일정한 구역이 두 배로 확장된 수렴관이 수치적으로 다루어졌고 이를 통하여 기존의 수렴관 에서 구한 값보다 큰 신장점도를 얻을 수 있었다.