动脉狭窄内低密度脂蛋白传输的数值研究：LDL的浓度极化现象

用计算机数值模拟的方法 ,对低密度脂蛋白 ( LDL)在动脉狭窄血管段内的质量传输进行了研究。计算结果表明 ,由于血管壁渗流的存在 ,LDL这样的脂质大分子会聚积在血管的内壁表面 ,发生一种工程上称为浓度极化的现象。LDL浓度在动脉狭窄口后的流动分离点出现峰值。该浓度峰值随雷诺数和动脉狭窄度的增加而呈逐渐下降的趋势。作者认为 ,该区域 LDL浓度的局部升高是引发动脉粥样硬化局部性和动脉狭窄产生的一个非常重要的原因。     

动脉分岔血管内膜增生过程的数值模拟

内膜增生从发生到阻塞血管是一个复杂的变化过程,在这个过程中,内膜的增生、血管腔体形状的改变和血流动力学之间是相互影响的。为了研究这些变化,本文提出一种单元填充方法数值模拟了三维颈动脉分岔血管在低切应力作用下血管内膜增生的过程。该方法既可以克服节点移动方法所不可避免的内膜增生的不连续性,也可以避免网格重划分的困难。结果发现,如果单纯以切应力阈值作为内膜增生的判据,低切应力的作用将无法导致血管完全阻塞,但内膜增生和血流动力学之间的相互影响是可以通过数值方法进行模拟的。在本数值模拟中,内膜增生的过程分为"增厚"(先)和"扩展"(后)两个阶段,最大狭窄率为34.4%,发生在距血管分岔5mm处动脉窦的外侧壁面。其发生位置和形状与临床观察吻合。     

飞秒激光作用下硅载流子浓度的超快变化过程数值模拟

根据硅材料的各参数和载流子浓度的速率方程,对超短脉冲激光辐照硅后栽流子浓度随时间的变化过程进行了数值计算。从硅载流子浓度驰豫速率方程出发,在单一载流子温度假设的基础上,推导出简化后的硅载流子浓度驰豫速率方程,并在Matlab平台上用4阶龙格-库塔法对微分方程进行了数值计算。数值计算和实验结果相似性较好。     

内爆炸下钢管混凝土结构变形与破坏的试验分析及数值模拟

钢管混凝土结构由于具有良好的力学性能和低经济成本而受到广泛应用,因此对它的研究很多,但对它在爆炸作用下尤其是内爆炸作用下的研究很少,故分析钢管混凝土结构在内爆炸下的变形与破坏具有重要意义.本文通过试验与数值模拟方法,分析了内爆炸下钢管混凝土结构在具有二条预制狭缝时的变形与破坏情形.结果表明,钢管变形主要发生在靠近狭缝处,变形大小随药量的增加而增大,最终导致钢管在狭缝两端出现约45°方向的裂纹,引起钢管混凝土剪切破坏.     

明渠条件下单丁坝绕流特征的数值模拟Numerical simulations of flows around a single spur dike in an open channel

丁坝结构广泛应用于水利工程中,用以调整水流和护滩固堤,维护优良的水道通航条件。针对单丁坝局部水流,基于非结构网格,采用有限体积法（FVM ）,建立了三维自由表面水流模型,湍流模型采用S‐A一方程模型。针对非淹没、正挑单丁坝,开展了系列 Fr数条件下的水流模拟。重点分析丁坝局部流动结构,探讨丁坝坝根处局部涡系演化等特征,研究了丁坝下游回流区长度和宽度的变化,总结了丁坝引起的剪切流的沿程变化特征。     

血流动力学与动脉粥样硬化

动脉粥样硬化引发的心血管疾病已成为疾病死亡的主要原因，动脉粥样硬化虽是一多因素疾病，但该疾病的局部性现象却与心血管系统的血流动力学特性有关。目前，人们在这一领域的研究重点主要集中于弄明局部应力，特别是流动剪切应力变化对动脉血管生理和病理的影响，本在综述血流动力学因素对动脉粥样硬化的影响的同时，根据实验研究提出了一些新的见解，认为血管内壁脂质浓度极化也是造成动脉粥样硬化的重要因素。     

颈动脉分支的血流动力学数值模拟

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

Computational technique for flow in blood vessels with porous effects

IntroductionThehemodynamicsofflowsthroughbloodvesselsisofgreatinterest,becausethesevesselspresentasubstantialhealthriskandareamajorcauseofmortalityandmorbidityintheindustrializedworld .Researchpapersonthebloodflowhaveappearedbutmostofthemhaveneglectedtheporosityeffectsduetovesselwalls.Inthisstudyweareinterestedintheflowthroughabloodvesseltakingintoaccounttheporosityeffectsofthevessels.Fluidflowthroughaporousmediumisoffundamentalimportancetowiderangeofdisciplinesinthevariousbranchesofnaturalsci…     

Computer Simulation of Non-Newtonian Flow and Mass Transport Through Coronary Arterial Stenosis

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｍｏｓｔｏｂｖｉｏｕｓｃｈａｒａｃｔｅｒｏｆａｔｈｅｒｏｓｃｌｅｒｏｓｉｓｉｓｔｈｅｌｏｃａｌｉｎｔｉｍａｌａｃｃｕｍｕｌａｔｉｏｎｏｆｌｏｗ＿ｄｅｎｓｉｔｙｌｉｐｏｐｒｏｔｅｉｎ(ＬＤＬ) ,ｔｈｅｉｎｃｒｅａｓｅｏｆｔｈｅｆｉｂｅｒａｌｔｉｓｓｕｅａｎｄｔｈｅｓｔｅｎｏｓｉｓ.Ｏｎｏｎｅｈａｎｄ ,ｉｔｉｓｆｏｕｎｄｔｈａｔｔｈｅｄｉｓｅａｓｅｏｆｔｅｎｏｃｃｕｒｓａｔｔｈｅｃｏｍｐｌｅｘｇｅｏｍｅｔｒｙｒｅｇｉｏｎ ,ｓｕｃｈａｓｂｉｆｕｒｃａｔｉｏｎｚｏｎｅ,ｃｕｒｅ…     

Some applications of perturbation method in thin plate bending problems

In this paper, problems of bending of thin plates under the combined action of lateral loading and in-plane forces are studied by means of perturbation method.     

Numerical simulation of turbulent flow through series stenoses

The flow fields in the neighbourhoods of series vascular stenoses are studied numerically for the Reynolds numbers from 100 to 4000, diameter constriction ratios of 0.2–0.6 and spacing ratios of 1, 2, 3, 4 and ∞. In this study, it has been further verified that in the laminar flow region, the numerical predictions by k–ω turbulence model matched those by the laminar‐flow modelling very well. This suggests that the k–ω turbulence model is capable of the prediction of the laminar flow as well as the prediction of the turbulent stenotic flow with good accuracy. The extent of the spreading of the recirculation region from the first stenosis and its effects on the flow field downstream of the second stenosis depend on the stenosis spacing ratio, constriction ratio and the Reynolds number. For c₁ = 0.5 with c₂ ≤ c₁, the peak value of wall vorticity generated by the second stenosis is always less than that generated by the first stenosis. However, the maximum centreline velocity and turbulence intensity at the second stenosis are higher than those at the first stenosis. In contrast, for c₁ = 0.5 with c₂ = 0.6, the maximum values at the second stenosis are much higher than those at the first stenosis whether for centreline velocity and turbulence intensity or for wall vorticity. The peak values of the wall vorticity and the centreline disturbance intensity both grow up with the Reynolds number increasing. The present study shows that the more stenoses can result in a lower critical Reynolds number that means an earlier occurrence of turbulence for the stenotic flows. Copyright © 2003 John Wiley & Sons, Ltd.     

窦部对称狭窄对颈动脉内流场影响的数值研究

以颈动脉分岔血管为例，采用数值方法研究了窦部环缩狭窄之后的流场分布情况，并和正 常血管情况下的流场分布进行了比较. 结果表明，采用环缩方式给颈动脉分岔血管施加对称 的狭窄改变了颈动脉窦内流场，特别是壁面剪应力的分布规律. 低剪应力区出现在狭窄段之 后的窦内，并且沿整个周向均匀分布. 根据低剪应力和动脉粥样硬化的关系，指出： 若人为地给颈动脉窦内施加对称狭窄，则脂质沉积将在狭窄下游的窦内沿周向轴对称 发展. 为了更真实地反映颈动脉窦内的狭窄，建议根据动脉血管中的实际狭窄情况，采用非 对称的狭窄分布模式.     

Brownian dynamics simulation of rodlike polymers under shear flow

The highly nonlinear behaviors of rodlike polymers in nematic phase under shear flow are studied with Brownian dynamics simulation. The LebwohlLasher nematogen model is taken as the prototype of the simulation and the mean-field approximation is avoided. By considering the nearest-neighbor intermolecular interaction, the spatial orientational correlation is introduced and therefore the spatial inhomogeneity such as the multiple-domain effect can automatically be incorporated. The transient order parameters, birefringence axes, shear stresses and first normal stress differences are calculated. The important finding of this work is that the director wagging and damped oscillation share the same molecular origin as director tumbling. The only difference is that the system is split into micro-domains which tumble with different phase angles in the wagging and damped oscillation regimes. The tumbling of the director of the whole system is suppressed due to the spatial inhomogeneity of director fields and then the damped oscillation of macroscopic stresses becomes predominant. The negative first normal stress difference exists at moderate shear rates, where both elasticity and viscosity play important role. Our simulation results including some dimensionless scaling parameters find good agreement with experimental observations in literature.