弯曲血管中的血液流动和大分子传质

采用数值模拟的方法,研究主动脉弯曲血管中的定常/脉动流动及低密度脂肪蛋白(LDL)和血清白蛋白(Albumin)传质．计算结果表明,对于主动脉弓模型,二次流漩涡的位置随时间变化．在弯曲变化比较剧烈的区域大分子浓度较高,壁面浓度外壁高于内壁．这些流动变化比较剧烈的区域可能是动脉硬化或血栓形成的危险区域．     

冠状动脉狭窄情况下的非牛顿血液流动和大分子传质

针对冠状动脉狭窄的情况,采用数值模拟方法求解了牛顿流体与非牛顿流体(幂次律流体和Casson流体)的定常与脉动的流场。在此基础上,求解了LDL(低密度脂肪蛋白)和Albumin(血清白蛋白)的浓度场。根据计算结果,详细讨论了壁面剪应力、非牛顿流效应、分子大小等因素对大分子传质的影响;并对牛顿流体与非牛顿流体、定常流动与脉动流动的大分子浓度场进行了比较,这些结果对于了解动脉硬化成因与流动特性和大分子传质的联系提供了较为丰富的信息。     

运动水平平板上混合对流的群论分析

本文用变换群理论对运动水平平板混合对流边界层流动的动量、能量和浓度扩散方程进行了分析，得到了与X4/(7-5n)成正比的壁面温度分布和浓度分布，同时壁面运动速度正比于X(3-n)／(7-5n)时存在相似性解．导出了相似性解方程，用四阶Runge－Kutta方法进行了计算，给出了Pr＝0．72　和Sc，K1，K2，K3参数下的速度、温度和浓度分布，得出了各参数对流场、温度场的影响。     

圆截面螺旋管道内非定常流动研究

以血液流动为背景,利用双参数摄动法研究了圆截面螺旋管内低频振荡流动,得到问题的二阶摄动解,分析了轴向速度、二次流、壁面剪应力在不同时刻的特点及随时间和Womersley数的变化情况。研究表明:挠率对圆截面曲线管道内低频振荡流动的影响不可忽略,尤其是轴向压力梯度绝对值很小时,挠率将对二次流动结构起主要影响作用。流函数的剧烈变化只发生在正、负数值发生转变的很小的时间段内,大部分时间段内变化平缓。壁面剪应力随θ的变化也很大。     

三维非定常水的跨血管壁流动与大分子的跨血管壁传质

对水的跨血管壁流动与大分子的跨血管对流-扩散传质提出一个统一的三维非定常模型,研究了内皮细胞在转换过程中形成漏的缝隙后的大分子跨壁传质,其浓度随时间的增长,内皮细胞截面形状对浓度分布的影响及生理参数对浓度分布的影响,采用解析方法求出了水跨壁流动的速度场与压力场的分析解,根据以前学者对血管壁的超微结构研究,将血管壁分为4层,用数值模拟方法求解了低密度脂蛋白（LDL）、血清白蛋白（albumin）以及辣根过氧化酶（HRP）三种大分子在各层的浓度分布,结果表明,随着时间的增长,在不同时刻,低剪切流中内皮细胞截面形状接近圆形时,大分子在漏的缝隙处浓度大小和分布范围,都大于高剪切流中内皮细胞截面被拉长时的情形;弹性层厚度等参数显著地影响传质,浓度分布随时间的增长与轴对称结果类似,即漏的缝隙处有一明显的浓度尖峰,在两侧浓度迅速减少,浓度尖峰值和浓度分布的范围随着时间的增长而降低与变小,这些结果对了解动脉粥样硬化形成的机理有重要意义。     

180°弯曲方管牛顿流体及一种非牛顿流体湍性流动的数值模拟

运用湍流k-ε模式及实测壁面函数分别模拟牛顿流体(清水)及一种非牛顿流体(聚合物稀薄减阻溶液)流经180°弯曲方管的湍性流动,取得与实测速度分布吻合较好的结果.对于湍流模式对存在大涡的复杂流动的适应性,根据计算和试验结果进行了分析和讨论.     

双向流固耦合作用下狭窄左冠状动脉内两相血流分析

基于血流与血管壁间双向流固耦合作用,将血液设为两相流体,运用计算流体力学方法对左冠状动脉内血流进行瞬态数值模拟.研究了一个心动周期内典型时刻下左冠状动脉内血流分布特性,并与Newton(牛顿)血液和两相血液模型对比,分析了两相血液和流固耦合作用对血流特性的影响.结果表明,左冠状动脉左回旋支远段和钝缘支近心端外侧分布了低速涡流区,该区域内壁面切应力和红细胞体积分数均较小,为动脉粥样硬化的形成与发展提供了合适的生理环境.左冠状动脉分叉处管壁形变量较大,引起管壁内膜功能发生紊乱,促进了粥样硬化斑块的形成.3种血液模型对比可知,红细胞的流动特性对血流速度及壁面切应力等血流动力学特性影响较大,双向流固耦合模型更符合真实的血液流动情况.     

低Mach数翼型绕流数值模拟的低耗散预处理格式

低Mach数流动中,基于可压缩流动的数值模拟算法存在严重的刚性问题,预处理方法可以有效地解决这一问题,但其计算结果不稳定．基于原有的预处理Roe格式,引入可调节参数,得到一种新的低耗散格式．该格式可以减弱边界层以及极低速区域的过度耗散,使得整个流场计算稳定．低Mach数、低Reynolds数定常圆柱绕流和低Mach数、高Reynolds数翼型(NACA0012和S809)绕流3个验证算例表明,带可调节参数的低耗散预处理方法正确可靠,是低速流动数值模拟的有效方法．     

精确的磁流体动力学汇流解析解

就一个特殊的磁流体动力学(MHD)流动,即速度幂指数为-1时的汇流,得到著名的Falkner-Skan方程精确的解析解.解析解是封闭的,并有多重解分支.分析了磁场参数和壁面伸长参数的影响.发现了有趣的速度分布现象:即使壁面固定,回流区域依然出现.在一个罕见的FalknerSkan MHD流动中,得到了一组解,以精确封闭的解析公式表示,极大地丰富了著名的Falkner-Skan方程的解析解,也加深了对这重要又有趣方程的理解.     

磁场和Hall电流对狭窄动脉中血液流动的影响

对一个水平向不对称、竖直向对称,带有轻微狭窄的动脉,提出了血液磁流体动力学流动的微极模型.为了估计狭窄形状的影响,几何上加以适当地考虑,通过选取不同的参数(称为形状参数),很方便地改变水平向的狭窄情况.在不同形状参数、Hartmann数和Hall参数下,计算了流动参数,例如流速、流动阻力(阻力阻抗)、狭窄区域血管壁面剪应力分布以及狭窄最大凸起高度位置处(狭窄喉部)的壁面剪应力大小.结果表明,流动阻力随着确定狭窄情况参数值和Hall参数值的增大而减小,并随着Hartmann数的增大而增大.对任意给定的Hartmann数和Hall参数,血管壁面剪应力和血管狭窄部位凸起最大高度处的管壁剪应力,具有与流动阻力相反的特征.最后,给出了Hartnmrm数和Hall参数对水平速度的影响.     

Effects of pulsatile flow and straight length on low-density lipoprotein transport in a curved arterial wall

Abnormal accumulation of macromolecules such as low-density lipoproteins (LDLs) in the arterial wall causes narrowing and blockage of vessels, which leads to atherosclerosis. Effects of pulsatile nature of blood flows as well as the initial length on transport of the LDL species in the arterial boundary layer region are analyzed numerically in the present work. The set of governing equations consisting of continuity, Navier-Stokes, and species transport is solved using a projection method based on the second-order central difference discretization. The obtained results are in excellent agreement with the pertinent data. The computational results imply that the flow field and concentration distribution are time dependent but the variation of the filtration velocity can be ignored. The LDL concentration boundary layer thickness decreases in the outer part and increases in the inner part for both with or without straight length. Presence of initial straight length generates about 26% growth in the boundary layer thickness, although its effect on the LDL surface concentration (LSC) is negligible. The maximum LSC is related to the regions with minimum wall shear stress in the inner part of the curved artery, which have more potential for formation of atherosclerosis. A new numerical correlation between the LSC and boundary layer thickness is proposed and examined.     

Momentum and heat transfer of the Falkner–Skan flow with algebraic decay: An analytical solution

In this paper, an analytical solution of the Falkner–Skan equation with mass transfer and wall movement is obtained for a special case, namely a velocity power index of ?1/3, with an algebraically decaying velocity profile. The solution is given in a closed form. Under different values of the mass transfer parameter, the wall can be fixed, moving in the same direction as the free stream, or opposite to the free stream (reversal flow near the wall). The thermal boundary layer solution is also presented with a closed form for a prescribed power-law wall temperature, which is expressed by the confluent hypergeometric function of the second kind. The temperature profiles and the wall temperature gradients are plotted. Interesting but complicated variation trends for certain combinations of controlling parameters are observed. Under certain parameter combinations, there exist singular points or poles for the wall temperature gradients, namely wall heat flux. With poles, the temperature profiles can cross the zero temperature line and become negative. From the results, it is also found empirically that under certain given values of the Prandtl number (Pr) and flow controlling parameter (b), the number of times for the temperature profiles crossing the zero line is the same as the number of poles when the wall temperature power index varies between zero and a given value n. The current result provides a new analytical solution for the Falkner–Skan equation with algebraic decay and greatly enriches the understanding of this important equation as well as the heat transfer characteristics for this flow.     

Computational modelling of slug flow in a capillary microreactor

The benefits of slug flow capillary microreactor exhibit the ability to adjust two individual transport mechanisms, i.e., convection inside the slug and diffusion between two consecutive slugs. The mass transfer rate is enhanced by internal circulation, which arises due to the shear between slug axis and continuous phase or capillary wall. The knowledge of circulation patterns within the slug plays an important role in the design of a capillary microreactor. Apart from this, well defined slug flow generation is a key activity in the development of methodology to study hydrodynamics and mass transfer. In the present paper we discuss computational fluid dynamics (CFD) modelling aspects of internal circulations (single phase) and slug flow generation (two-phase).     

Mathematical Study of an Inflammatory Model for Atherosclerosis: A Nonlinear Renewal Equation

In this work we study the coupling of a nonlinear renewal equation to an ordinary differential equation. We start with existence and uniqueness issues for the coupled equations and, in particular cases, we study the long-time behaviour. The novelty here is the nonlinearity in the renewal equation. This model arises in the context of atherosclerosis. The renewal part accounts for the inflammatory process: leucocyte recruitment in the arterial wall, differentiation when internalizing low-density lipoprotein (LDL) and death. The ordinary differential equation describes the LDL dynamics in the arterial wall, leucocyte absorption and release in the blood.

     

Factors Influencing the Disturbed Flow Patterns Downstream of Curved Atherosclerotic Arteries

Pulsatile blood flows in curved atherosclerotic arteries are studied by com- puter simulations.Computations are carried out with various values of physiological parameters to examine the effects of flow parameters on the disturbed flow patterns downstream of a curved artery with a stenosis at the inner wall.The numerical re- suits indicate a strong dependence of flow pattern on the blood viscosity and inlet flow rate,while the influence of the inlet flow profile to the flow pattem in downstream is negligible.     

Effects of variable viscosity on pulsatile flow of blood in a tapered stenotic flexible artery

The objective of the present study is to investigate the effects of variable viscosity on incompressible laminar pulsatile flow of blood through an overlapping doubly constricted tapered artery. To mimic the realistic situation, wall of the artery is taken to be flexible, and physiologically relevant pulsatile flow is introduced. The governing equations of blood flow are made dimensionless. A coordinate transformation is used to make the overlapping doubly constricted wall geometry of tube to a straight tube. Taking advantage of the Stream function–Vorticity formulation, the system of partial differential equations is then solved numerically by finite difference approximations. Effects of Reynolds number, Strouhal number, degree of contraction, tapering angle, and viscosity parameters are presented graphically and analyzed. The results show that formation of stenosis and tapering disturb the flow field significantly, and degree of stenosis is more important in influencing blood flow compared with tapering.     

Viscous flow over a shrinking sheet with an arbitrary surface velocity

In this paper, an analytical solution in a closed form for the boundary layer flow over a shrinking sheet is presented when arbitrary velocity distributions are applied on the shrinking sheet. The solutions with seven typical velocity profiles are derived based on a general closed form expression. Such flow is usually not self-similar and the solution can only be implemented when the mass transfer at the wall is prescribed and determined by the moving velocity of the wall. The characteristics of the flows with the typical velocity distributions are discussed and compared with previous similarity solutions. The flow is observed to have quite different behavior from that of the self-similar flow reported in the literature and the results demonstrate distinctive momentum and energy transport characteristics. Some plots of the stream functions are also illustrated to show the difference in flow field between the shrinking sheet and the stretching sheet. An integral approach to solve boundary layer flow over a shrinking or stretching sheet with uncoupled arbitrary surface velocity and wall mass transfer velocity is outlined and the effectiveness of this approach is discussed.