首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 109 毫秒
1.
本文考虑了血液的菲牛顿流动特性,对血液在动脉管系中的脉动流建立了准二维流动模型。利用有限差分方法得到了动脉管系内非牛顿流体的准二维不定常流动的数值解。并以人体五根主要动脉所组成的动脉管系为例,进行了详细的数值计算。计算结果表明,在动脉的某些典型位置上。计算所得的理论波形与实测波形是相似的。  相似文献   

2.
曹骞  康灿  滕爽  焦侬  丁可金 《摩擦学学报》2022,42(6):1094-1104
为研究弯管内固体颗粒在液相夹带条件下的运动特性及颗粒对弯管内壁的磨损,采用计算流体动力学与离散元耦合的方法,建立数值模型,考虑固液两相之间的作用,对弯管内的固液两相流动进行数值模拟;通过软件的应用程序编程接口嵌入自编译磨损模型;借助试验结果,验证数值模型的有效性. 结果表明,所建立的数值模拟方案可以准确地模拟颗粒在管内的运动特征并能够预测弯管内壁的磨损位置以及磨损程度. 弯管内的二次流对颗粒运动有重要影响,弯管外侧壁面中心线附近的磨损较严重,磨损的形式以小角度划擦切削为主. 弯管磨损主要与颗粒对壁面的碰撞速度、碰撞角度及碰撞频率有关. 运动中的颗粒与壁面发生多次碰撞,碰撞角度逐渐减小. 随着颗粒球形度的增大,在相同碰撞条件下引起的磨损量变小,但是会降低颗粒的随流性. 颗粒形状影响颗粒在流场中的运动速度以及颗粒与壁面的碰撞. 随着颗粒球形度增大,严重磨损区域向弯管进口方向移动,壁面平均磨损量先减小后增大;当输送颗粒的球形度为0.91时,壁面磨损量最小.   相似文献   

3.
颈动脉分支的血流动力学数值模拟   总被引:7,自引:0,他引:7  
采用有限元法数值模拟颈动脉分支的血流动力学。根据在体测量的实际尺寸来构造颈动脉分支的几何模型,以保持模型的解剖精确度;利用在体测量的颈内动脉和颈外动脉流量波形以及主颈动脉的压力波形来确定数值计算的边界条件,以保持数值计算的生理真实性。关注的重点是颈动脉窦内的局部血流形态、二次流和壁面剪应力。在心脏收缩的减速期和舒张期的某些时刻,颈动脉窦中部外侧壁面附近产生了流动分离,形成了一个低速回流区。该流动分离是瞬态的,导致了壁面剪应力的振荡,其振荡范围在-2~6dyn/cm^2之间。同时,颈动脉窦中部横截面内的二次流存在于整个心动周期,最大的二次流速度为同时刻轴向速度平均值的1/3左右。  相似文献   

4.
局部狭窄血管中血液振荡流的速度和切应力   总被引:1,自引:1,他引:0  
孙辉  柳兆荣 《力学季刊》2002,23(2):148-156
本文建立一种分析局部缓慢狭窄血管中血液振荡流的数学模型,给出了血液的轴向流速,径向流速和切应力的包含压力梯度项的解析表达式,并讨论了血管内由局部狭窄引起的压力梯度沿轴向变化的规律。文章以局部余弦狭窄为例进行数值计算,详细讨论上游均匀管段压力梯度的定常部分和不同次谐波对狭窄管段内流速和切应力的影响。数值结果表明,与均匀管情况相比,在狭窄段内,血液振荡流轴向流速无论平均值还是脉动幅值均明显增大,且径向流速不再为零。但径向流速仍远小于轴向流速。同时,切应力也不再仅由轴向流速梯度提供,径向流速梯度也将产生切应力,但是在计算管壁切向上的切应力时,径向流速梯度的贡献仍相当大。与均匀管管壁切应力沿流运方向保持恒定不同。狭窄管管壁切应力(平均值和脉动值)将随着狭窄高度的增大而增大,在狭窄最大高度处达到最大,因而沿流动方向产生了较大的切应力梯度。  相似文献   

5.
平行平板流动腔脉动流切应力的计算   总被引:5,自引:0,他引:5  
高度远小于横向和纵向几何尺寸的矩形平行平板流动腔是人们用以体外研究细胞在切应力作用下力学行为的主要工具之一。大多数研究者主要对定常层流情进行研究。本文通过对矩形平行平板流动腔内的层流脉动流进行详细分析,给出腔内速和腔室底部切应力的准确计算公式。当Womersley数较小时,准确公式简化为准定常公式。数值计算结果表明,在脉动流条件下,对于人们常用的流动腔几何尺寸,准定常公式具有相当高的精度。这为在脉  相似文献   

6.
90°弯管内流动的理论模型及流动特性的数值研究   总被引:27,自引:0,他引:27  
从三维不可压缩雷诺时均Navier-Stokes方程出发,对90°弯曲管道内湍流流动进行数值模拟。网格划分采用六面体网格,湍流模型为RNGk-ε模型,在近壁区采用两层壁面模型进行修正,流场的计算结果与实验数据吻合较好。在此基础上,本文数值研究了来流方向对流场结构和流动特性的影响。得出在弯管流场中发生了分离现象,且随着来流侧滑角的增大,分离区范围增大。此外,随着来流从同一侧滑角变换至同一攻角时,横截面的二次流图像中也从具有两个对称主涡变成只具有一个主涡的现象。  相似文献   

7.
在脉动流条件下,用计算机数值模拟的方法对低密度脂蛋白(LDL)在动脉狭窄血管段内的质量传输进行了研究。计算结果表明.无论是在定常流还是在脉动流条件下.LDL都将聚积于血管狭窄处峰口附近的流动分离点,LDL壁面浓度在此处最高。在脉动流条件下,LDL在血流受扰动区的聚积高于定常流的值;而且.流动分离点处LDL壁面浓度峰值覆盖的区域也宽于定常流。本文所揭示出的LDL在血管狭窄处的质量传输现象可能在动脉粥样硬化的局部性和动脉狭窄的形成中起着很重要的作用。  相似文献   

8.
Taylor补丁对新型动脉旁路移植流场影响的数值分析   总被引:1,自引:0,他引:1  
为了研究Taylor补丁对新型(S型)动脉旁路移植术中吻合口处流场的影响,使用数值方法研 究了采用Taylor补丁和未采用该补丁的两个S型旁路移植模型内流场的血流动力学差异. 对 流速、壁面切应力和切应力梯度等参数进行了比较分析. 结果表明,Taylor补丁对吻合口的 流场有显著影响. 采用Taylor补丁的模型其下游吻合口处的流场分布较未采用补丁的模型更 均匀,二次流平均流速减小约34.48%,壁面切应力梯度减小约52.22%,从壁面切应力 梯度方面分析,这将有助于改善血流动力学分布,抑制动脉粥样硬化. 但从壁面切应力值分 析,其动脉底部的壁面低切应力区明显增大,平均壁面切应力值减小30.33%,这又将促 使动脉粥样硬化. 因此,Taylor补丁是否对S型搭桥术具有治疗优越性,仅从血流动力学分 析尚不能定论,配合数值计算结果进行动物和临床实验研究是十分必要的.  相似文献   

9.
一种确定均匀动脉壁面切应力的非线性方法   总被引:4,自引:0,他引:4  
覃开蓉  姜宗来 《力学学报》2005,37(2):225-231
从Ling和Atabek提出的``局部流'理论出发,提出一种利用测量血液黏度、管轴上 的血流速度、压力和管径波形计算均匀动脉管壁切应力的非线性方法. 将这种方法与柳兆荣 等提出的利用测量血液黏度、管轴上的血流速度和平均管径计算切应力的线性方法比较,结 果表明,当管壁脉动幅度较小时,两种方法计算的压力梯度、流速剖面和管壁切应力差别较 小;而当管壁脉动幅度增大时,两种方法计算的压力梯度、流速剖面和管壁切应力差别增大. 对于小幅脉动均匀动脉,用线性方法计算管壁切应力有较高的精度;而对于大变形 均匀动脉,则需要考虑非线性因素对管壁切应力的影响. 由于作为输入量的血液黏度、轴心 血流速度、压力波形和管径波形可在活体上通过无损伤或微损伤的检测方法得到, 所提出的计算切应力的方法为在体或离体研究切应力与动脉重建的关系提供了方法学基础.  相似文献   

10.
就我们所熟知,绝大部分正常动脉流,其血液的流动特性是属于层流范围,但随着弯曲和分 支部分会产生血液流之二次回流区,进而形成所谓近似非稳态流及紊流. 因此动脉流体的特 性会随动脉外形及条件的改变而改变. 在某些情形下,异常动脉的血液动力特性会造成动脉 的病变. 因此,近年来动脉血液流体的特性的研究,常着重于异常动脉的血液动力特性所形 成剪应力和病变部位动脉粥状硬化关系的探讨. 动脉血液流动经常包含分离流或二次回流运动,而这是流体力学的分析或数值模拟最困难的 部分. 有关分离流或二次回流的研究包括正常血管流和异常血管流,藉由二次回流的模拟与 测量可以观察血管病变的形成与演变,其中最受注目探讨题目是窄化血管如粥状斑块相关的 血液流动分析. 将回顾二维和三维、稳态、非稳态之动脉血流与窄化血管相关的几何外形作模拟研究和 实验. 并提供对血液动力学的研究方向,以作为未来医疗诊断与发展相关器材之参考.  相似文献   

11.
The incompressible laminar flow of air and heat transfer in a channel with a backward-facing step is studied for steady cases and for pulsatile inlet conditions. For steady flows the influence of the inlet velocity profile, the height of the step and the Reynolds number on the reattachment length is investigated. A parabolic entrance profile was used for pulsatile flow. It was found with amplitude of oscillation of one by Re=100 that the primary vortex breakdown through one pulsatile cycle. The wall shear rate in the separation zone varied markedly with pulsatile flows and the wall heat transfer remained relatively constant. The time-average pulsatile heat transfer at the walls was greater as with steady flow with the same mean Reynolds number.  相似文献   

12.
利用两相湍流KET模型对90°弯管内气固两相湍流流动进行了数值模拟,得到了弯管内两相流动的一些规律,并提出用颗粒动理学压力来定性表征弯管内磨损严重部位,为管道抗磨损设计提供了一定的理论依据。  相似文献   

13.
动脉瘤内流场以及瘤体尺寸的影响的数值研究   总被引:1,自引:0,他引:1  
采用计算流体动力学(CFD)数值模拟的方法,在周期性脉动速度入流条件下,建立刚性动脉瘤模型并研究了动脉瘤模型中流场的特征(速度、压力、壁面剪切应力)。得到了脉动入流一个周期内流场特征的变化规律,发现动脉瘤的后端有相当高的压力和壁面剪切应力,而且高压力和壁面剪切应力分布的位置几乎是固定的。探讨了不同动脉瘤尺寸对内部流场的影响,动脉瘤的直径与瘤体长度之比越大,瘤壁承受的剪切应力就越大,动脉瘤破裂的危险性就越高。  相似文献   

14.
A numerical analysis of the steady/pulsatile flow and macromolecular (such as LDL and Albumin) transport in curved blood vessels was carried out. The computational results predict that the vortex of the secondary flow is time-dependent in the aortic arch. The concentration of macromolecule concentrates at the region of sharp curve, and the wall concentration at the outer part is higher than that at the inner part. Atherosclerosis and thrombosis are prone to develop in such regions with sharp flow.  相似文献   

15.
The velocity field in a patient-specific abdominal aneurysm model including the aorto–iliac bifurcation was measured by 2D PIV. Phase-averaged velocities obtained in 14 planes reveal details of the flow evolution during a cycle. The aneurysm expanding asymmetrically toward the anterior side of the aorta causes the generation of a vortex at its entrance, covering the entire aneurysm bulge progressively before flow peak. The fluid entering the aneurysm impinges on the left side of its distal end, following the axis of the upstream aorta segment, causing an increased flow rate in the left (compared to the right) common iliac artery. High shear stresses appear at the aneurysm inlet and outlet as well as along the posterior wall, varying proportionally to the flow rate. At the same regions, elevated flow disturbances are observed, being intensified at flow peak and during the deceleration phase. Low shear stresses are present in the recirculation region, being two orders of magnitude smaller than the previous ones. At flow peak and during the deceleration phase, a clockwise swirling motion (viewed from the inlet) is present in the aneurysm due to the out of plane curvature of the aorta.  相似文献   

16.
 A laminar wall jet undergoing transition is investigated using the particle image velocimetry (PIV) technique. The plane wall jet is issued from a rectangular channel, with the jet-exit velocity profile being parabolic. The Reynolds number, based on the exit mean velocity and the channel width, is 1450. To aid the understanding of the global flow features, laser-sheet/smoke flow visualizations are performed along streamwise, spanwise, and cross-stream directions. Surface pressure measurements are made to correlate the instantaneous vorticity distribution with the surface pressure fluctuations. The instantaneous velocity and vorticity field measurements provide the basis for understanding the formation of the inner-region vortex and the subsequent interactions between the outer-region (free-shear-layer region) and inner-region (boundary-layer region) vortical structures. Results show that under the influence of the free-shear-layer vortex, the local boundary layer becomes detached from the surface and inviscidly unstable, and a vortex is formed in the inner region. Once this vortex has formed, the free-shear-layer vortex and the inner-region vortex form a vortex couple and convect downstream. The mutual interactions between these inner- and outer-region vortical structures dominate the transition process. Farther downstream, the emergence of the three-dimensional structure in the free shear layer initiates complete breakdown of the flow. Received: 8 November 1995/Accepted: 6 November 1996  相似文献   

17.
Numerical Simulation of the Flow in the Carotid Bifurcation   总被引:3,自引:0,他引:3  
Pulsatile flow through the three-dimensional carotid artery bifurcation has been studied using the artificial-compressibility method. The part of the flow with large inertia bifurcates and creates a very steep velocity gradient on the divider walls. The flow near the nondivider walls slows down because of dilation of the cross section and strong adverse pressure gradient. The secondary flow in the bifurcation region, which is similar to the Dean vortex in a curved pipe, is strong and very complex. The region of separation is not closed for the cases of steady and pulsatile flow. The extent of this region is small and the streamlines are smooth except in the decelerating phase of systole. The change of common-internal bifurcation angle (25°± 15°) for fixed internal–external bifurcation angle of 50° has more effect on the shear on the bifurcation-internal carotid wall and less effect on the shear on the common-internal carotid wall. The mean wall shears are not sensitive to the input flow-rate waveform for constant mean flow, but the maximum wall shears are. Received 3 January 1997 and accepted 11 April 1997  相似文献   

18.
If the torque exerted on a fluid element and the source of streamwise vorticity generation are analyzed, a turbulence-driven secondary flow is found to be possible in a curved pipe. Based on this analysis, it is found that the secondary flow is primarily induced by high anisotropy of the cross-stream turbulent normal stresses near the outer bend (furthest from the center of curvature of the bend). This secondary flow appears as a counterrotating vortex pair embedded in a Dean-type secondary motion. Recent hot-wire measurements provide some evidence for the existence of this vortex pair. To verify the formation and extent of this turbulence-driven vortex pair further, a near-wall Reynolds-stress model is used to carry out a detailed numerical investigation of a curved-pipe flow. The computation is performed specifically for a U-bend with a full developed turbulent flow at the bend entrance and a long straight pipe attached to the exit. Numerical results reveal that there are three vortex pairs in a curved pipe. The primary one is the Dean-type vortex pair. Another pair exists near the pipe core and is a consequence of local pressure imbalance. A third pair is found near the outer bend and is the turbulence-driven secondary flow. It starts to appear around 60° from the bend entrance, grows to a maximum strength at the bend exit, and disappears altogether at about seven pipe diameters downstream of the bend. On the other hand, calculations of developing laminar curved-pipe flows covering a range of pipe-to-bend curvature ratios, Reynolds number, and different inlet conditions fail to give rise to a third cell near the outer bend. Therefore, experimental and numerical evidence together lend support to the formation of a pair of turbulence-driven secondary cells in curved-pipe flows.Research supported by the Office of Naval Research under Grant No. N0014-81-K-0428 and by the David Taylor Research Center, Annapolis, Maryland, under Contract No. N00167-86-K0075.  相似文献   

19.
In this research, the developing turbulent swirling flow in the entrance region of a pipe is investigated analytically by using the boundary layer integral method. The governing equations are integrated through the boundary layer and obtained differential equations are solved with forth-order Adams predictor-corrector method. The general tangential velocity is applied at the inlet region to consider both free and forced vortex velocity profiles. The comparison between present model and available experimental data demonstrates the capability of the model in predicting boundary layer parameters (e.g. boundary layer growth, shear rate and swirl intensity decay rate). Analytical results showed that the free vortex velocity profile can better predict the boundary layer parameters in the entrance region than in the forced one. Also, effects of pressure gradient inside the boundary layer is investigated and showed that if pressure gradient is ignored inside the boundary layer, results deviate greatly from the experimental data.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号