• 설비공학논문집
• /
• 제6권3호
• /
• pp.227-236
• /
• 1994

Steady flows of Newtonian and non-Newtonian fluids in the stenotic tubes with various stenotic shapes are numerically simulated. Validity of the modified power-law model as a constitutive equation for the purely viscous non-Newtonian fluid is discussed and the results of the power-law model are compared with those of the Carreau model, the Powell-Eyring model and experimental data for blood. Flow characteristics and reattachment lengths for non-Newtonian fluids in the stenotic tubes are presented extensively. Also, the analysis is extended to predict the influences of diameter ratio, stenosis spacing, number of stenosis and Reynolds number on the flow characteristics in the multiple stenotic tubes.

• 한국산학기술학회논문지
• /
• 제18권3호
• /
• pp.545-550
• /
• 2017

대동맥이나 협착된 경동맥에서는 심장수축기에 간헐적으로 난류현상이 발생하고 있으며, 혈액의 점성특성으로 인해 기존 난류모델로는 정확한 해석이 어려운 실정이다. 혈류는 점탄성 유체의 성질을 가지고 있어 유체의 전단 변형률 증가에 따라 점도가 감소하는 점탄성 유체이며, 이러한 점탄성 유체는 난류 유동시 저항 감소 현상이 발생한다. 기존의 난류해석 모델들은 점성변화가 없는 뉴턴 유체에 적합한 모델들이 대부분이기 때문에, 점탄성 유체의 저항 감소 현상을 고려한 비뉴턴 유체 해석에 적합한 난류 모델개발이 필요하다. 본 논문은 난류 모델 가운데 수렴성이 좋고 해석시간이 짧은 표준 $k-{\varepsilon}$ 모델을 기반으로 저항 완충 함수를 이용하여 비뉴턴 유체의 저항감소 현상을 해석할 수 있는 수정된 난류모델을 제시하였으며, 이를 기존 난류모델들과 비교하여 제시된 난류 모델을 검증하였다. 새로 제시된 수정된 난류모델은 벽함수 및 점성저층을 고려하지 않았기 때문에 해석시간이 대폭적으로 감소하였으며, 적은 격자수를 이용하여 효율적으로 비뉴턴 유체의 난류 현상을 해석할 수 있기 때문에 향후 혈류해석 및 점탄성유체 해석에 적용할 예정이다.

• 대한의용생체공학회:의공학회지
• /
• 제21권3호
• /
• pp.261-271
• /
• 2000

벽면운동(wall motion)과 임피던스 페이즈앵글(impedance phase angle; 압력파와 유랑파 기아의 위상차)을 고려하여 맥동유동하에 있는 복부대동맥 분기관모델에서 2차원 전산유체해석을 수행하였다. 해석결과 분기광 근처에서응 전단응력의 크기가 매우 급격한 변화를 보임을 관찰하였고 벽면운동은 전단응력의 진폭을 감소시키는 효과를 가져왔다. 임피던스 페이즈 앵글이 음의 값을 향해 갈수록 시간 평균된 벽면 전단응력(mean wall shear stress)의 값은 감소하였으나 진폭(amplitude of wall shear stress)은 오히려 증가하였다. 페이즈앵글의 영향은 평균 벽면전단응력이 영에 근접하는 외벽(outer wall or lateral wall)의 바같쪽으로 휘어지는 부분(curvature site)에서 상대적으로 크게 나타났는데 $-90^{\circ}$ 페이즈앵글(혈류파가 혈압파를 1/4주기 앞서는 경우)일 경우에 $0^{\circ}$의 경우에 비해 평균은 $50\%$정도 감소하였고 진폭은 $15\%$정도의 상승를 나타내었다. 그러므로 고혈압 환자와 같이 큰 음의 페이즈앵글을 갖는 경우, 벽면전단응력의 평균은 낮아지고 시간에 따라 변화량(진폭)은 증가하므로 low and oscillatory wall shear stress 이론에 의하면 동맥경화에 더 민감하게 된다. 비뉴턴유체로 모델링한 경우에는 뉴턴유체의 경우에 비해 벽면전단응력의 평균값이 증가하므로서 동맥경화에 덜 민감하게 된다.

• 대한기계학회논문집B
• /
• 제20권11호
• /
• pp.3686-3694
• /
• 1996

A numerical investigation of natural convection flow along irregular vertical surfaces is reported. A transformation method is applied to the problem of natural convection under the assumption of a large Grashof number. A vertical wavy surface is used as an example to demonstrate the advantages of the transformation method, and to show the heat transfer mechanism near such surfaces. Surface non-uniformities on the boundary layer flow induced by a constant was temperature, semi-infinite surface are investigated. Also the effects of Prandtl number, flow index, and surface amplitude in Non-Newtonian fluids are discussed. When possible, the comparison of the numerical results shows a good agreement. The amplitude is proportional to the amplitude of a wavy surface. The results demonstrate that the local heat flux along a wavy surface is smaller than that of a flat surface. The frequency of the wavy surface is half that of the local heat transfer rate. The amplitude of the local Nusselt number gradually decreases downstream where the natural convection boundary layer grows thick.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2010년 춘계학술대회논문집
• /
• pp.214-220
• /
• 2010

This article described that a high Reynolds number version of a turbulence model was modified by using drag reduction to analyze the turbulent flows of non-Newtonian fluid with visco-elastic viscosity and it was applied hemodynamics which was representative of visco-elastic fluid. The turbulence characteristics of visco-elastic fluid was expanded viscous sublayer region and buffer layer region by drag reduction phenomenon and also Newtonian turbulence models does not predict because viscosity was related with shear rate of fluid flow. Hence numerical simulation using a modified turbulence model was conducted under the same conditions that were applied to obtain the experiment results and previous turbulence models and then the numerical investigation of turbulent blood flow in the stenosed artery bifurcation under periodic acceleration of the human body.

• 유변학
• /
• 제6권2호
• /
• pp.157-164
• /
• 1994

인체 동맥혈관내 혈액의 유동현상을 수치적으로 해석하기 위해서는 혈액의 유변학 적 성질을 구성방정식으로 나타내어야한다. 본 연구에서는 혈액의 점성계수를 표현하기 위 하여 비뉴턴유체의 점성을 나타내는 식으로서 Carreau 모델, 수정 Cross 모델, 수정 Powell-Eyring 모델과 수정멱법칙모델을 사용하였고 원형관내 혈액의 정상유동을 수치모사 하기 위하여 겉보기점성계수를 이용하는 구성방정식을 운동량방정식에 적용하였다. Carreau 모델과 수정멱법칙모델을 적용할 때 레이놀즈수의 변화가 중심선상의 속도와 길이방향의 압 력변화에 미치는 영향을 고찰하였다. 전단율이 높은영역에서 혈액의 겉보기점성계수를 효과 적으로 나타낼수 있는 수정멱법칙모델을 제시하였다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년도 춘계학술대회논문집
• /
• pp.385-388
• /
• 2008

The objective of this paper is simulating blood flow through the branched and stenotic tube numerically. SC-Tetra, which is one of the commercial code using FVM method, was utilized for this analysis. The flow is assumed as an incompressible laminar flow with the additional condition of non-Newtonian fluid. As the constitutive equation for the fluid viscosity, the following models were solved with governing equations ; Cross Model, Modified Cross Model, Carreau Model and Carreau-Yasuda Model. Final goal was achieved to get analytic data about shear stress, at specific points, changing the geometry with various factors like the bifurcation angle, diameter of the branches, the ratio of stenosis, and etc. The material property of blood was referred from the related papers. Furthermore, to verify results they were compared with those of the published papers. There were some discrepancies based on the different solver and the different data post-processing method. However, many parameters like the location of low shear stress, which arised from bifurcation or stenosis, and the tendency of various factors were found to be very similar.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년도 학술대회
• /
• pp.177-180
• /
• 2008

The objective of this study is to get simulation data about pulsatile flow of a non-Newtonian fluid through a bifurcated tube. All the process was based on CFD method, with a commercial FVM code, SC/Tetra ver. 6.0 for solving, and with CATIA R16 for generating geometries. To define a non-Newtonian fluid, the following viscous models are used; the Powell-Eyring model, the modified Powell-Eyring model, the Cross model, the modified Cross model, the Carreau model, the Carreau-Yasuda model and the modified Power Law model. The flow calculation data using each model were compared with the other data of a existing paper. Finally, the Carreau model was recognized to give the best result with the SC/Tetra code, and the succeeding simulations are made with the model. For the pulsating flow condition, the sine wave type velocity profile is given as the inlet boundary condition. To investigate the effect of geometries and mesh, the pre-test is carried out with various curvature conditions of the bifurcated corner, and then with various mesh conditions. The final process is to calculate flow variables such as the wall shear stress (WSS) and the wall shear stress gradient (WSSG). To validate all the result, the simulation is compared with the existing data of the other papers. Generally speaking, there is a noticeable difference in the maximum and minimum value of WSS. It is not sure that the values in each data are on the exactly same location. However, the overall trend is similar. The next study needs to investigate the same situation by experimental method. Furthermore, if the flow is simulated with more pulsatile conditions, more data of flow field through a bifurcated tube could be achieved.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년 추계학술대회논문집
• /
• pp.177-180
• /
• 2008

The objective of this study is to get simulation data about pulsatile flow of a non-Newtonian fluid through a bifurcated tube. All the process was based on CFD method, with a commercial FVM code, SC/Tetra ver. 6.0 for solving, and with CATIA R16 for generating geometries. To define a non-Newtonian fluid, the following viscous models are used; the Powell-Eyring model, the modified Powell-Eyring model, the Cross model, the modified Cross model, the Carreau model, the Carreau-Yasuda model and the modified Power Law model. The flow calculation data using each model were compared with the other data of a existing paper. Finally, the Carreau model was recognized to give the best result with the SC/Tetra code, and the succeeding simulations are made with the model. For the pulsating flow condition, the sine wave type velocity profile is given as the inlet boundary condition. To investigate the effect of geometries and mesh, the pre-test is carried out with various curvature conditions of the bifurcated corner, and then with various mesh conditions. The final process is to calculate flow variables such as the wall shear stress (WSS) and the wall shear stress gradient (WSSG). To validate all the result, the simulation is compared with the existing data of the other papers. Generally speaking, there is a noticeable difference in the maximum and minimum value of WSS. It is not sure that the values in each data are on the exactly same location. However, the overall trend is similar. The next study needs to investigate the same situation by experimental method. Furthermore, if the flow is simulated with more pulsatile conditions, more data of flow field through a bifurcated tube could be achieved.

• 대한기계학회논문집B
• /
• 제25권1호
• /
• pp.18-28
• /
• 2001

The present study investigates flow dynamics between two dimensional compliant plates under sinusoidal flow conditions in order to understand influence of wall motion, impedance phase angle (time delay between pressure and flow waveforms), and non-Newtonian fluid on wall shear stress using computational fluid dynamics. The results showed that wall motion induced additional terms in the streamwise velocity profile and the pressure gradient. These additional terms due to wall motion reduced the amplitude of wall shear stress and also changed the mean wall shear stress. The trend of the changes was very different depending on the impedance phase angle. As the impedance phase angle was changed to more negative values, the mean wall shear stress decreased while the amplitude of wall shear stress increased. As the phase angle was reduced from 0°to -90°under $\pm$4% wall motion, the mean wall shear stress decreased by 12% and the amplitude of wall shear stress increased by 9%. Therefore, for hypertensive patients who have large negative phase angles, the ratio of amplitude and mean of the wall shear stress is raised resulting in a more vulnerable state to atherosclerosis according to the low and oscillatory shear stress theory. We also found that non-Newtonian characteristics of the blood protect atherosclerosis by decreasing the oscillatory shear index.