等宽多孔介质壁面管道中磁流体的流动

研究等宽管道中,磁场、可渗透壁面、Darcy速度和滑动参数,对流体稳定流动的综合影响.假设管道中流动的流体是均匀的、不可压缩的Newton流体.利用Beavers-Joseph滑动边界条件,得到控制方程的解析解.详细地讨论了磁场、可渗透性、Darcy速度和滑动参数对轴向速度、滑动速度和剪应力的影响.可以看出,Hartmann数、Darcy速度、多孔参数和滑动参数,在改变流动方向,进而改变剪应力方面,起着至关重要的作用.     

流变学流体的蠕动传输:食道中食物块的运动模型

研究食道中蠕动传输的流体力学.对任意的波形和任意的管道长度,建立起流变学流体蠕动传输的数学模型.用粘性流体的Ostwald-de Waele幂定律,描述非Newton流体的流动特性.解析公式化模型,详细且精确地给出食物块在食道中蠕动传输相关的一些重要性质.分析中应用了润滑理论,本研究特别适合于Reynolds数不大的情况.将食道看作环形的管道,通过食道壁周期性的收缩来传输食物块.就单个波和周期性收缩一组波的传播,研究与传输过程有关变量的变化,如压力、流速、食物颗粒轨迹以及流量等.局部压力的变化,对流变指数n有着高度的敏感性.研究结果清晰地表明,食物块在食道中蠕动传输时,Newton流体或流变学流体构成的连续流体,以组合波传播比大间隔单波传播,传输效率要高得多.     

具有壁面粗糙度的平行板微管道内三阶流体的电磁流动

研究了具有粗糙壁面的平行板微管道内三阶流体的电磁驱动流.假设两个壁面粗糙度的形状是相位差为0或π的小振幅正弦波形状.将洛伦兹力作为体积力,利用摄动法解析求出了速度和流率的近似解.通过数值计算,结果表明随着波数或非牛顿参数的增加,壁面粗糙度对三阶流体的阻力增加.随着Hartmann数的增加,壁面粗糙度对三阶流体的阻力减小.相位差为0的壁面粗糙度对流动的阻力大于相位差为π的粗糙度对流动的阻力.当波数或Hartmann数充分大时,壁面粗糙度的相位差变得不太重要.     

不对称柔性壁管道内幂律流体蠕动传输的精确解

在不对称管道内,研究了壁面柔曲性对非Newton流体蠕动流的影响.流变学性质由幂律流体本构方程表征.在数学表达中,采用了长波和低Reynolds数近似.得到了流函数和速度的精确解.给出了流线图及其俘获现象.对所讨论的流动,陈列了关键参数的显著特征,并最后给出了主要结论.     

磁场和热辐射对可变表面热通量作用下的竖直圆锥体自然对流的影响

就圆锥体表面受到可变表面热通量作用,计及磁场和热辐射的综合影响,数值研究了流经竖直圆锥体的自然对流及其热交换特点.认为流体是灰色的、吸收-发射的辐射介质,而非散射介质,通过近似变换,将自由对流区中流动的边界层控制方程,简化为无量纲方程.利用Crank-Nicol-son形式的隐式有限差分法(具有收敛快、精度高、无条件稳定的特点),求解了无量纲的控制方程.得到了数值结果,以及空气和水中的速度、温度、局部和平均的壁面剪应力、局部和平均的Nusselt数.将所得到的结果与先前文献报道的结果进行比较,发现两者有着很好的一致性.     

T型微通道内液滴在幂律流体中运动机理的格子Boltzmann方法研究北大核心CSCD

采用非Newton不可压两相流格子Boltzmann模型研究了T型微通道内Newton液滴在非Newton幂律流体中的运动过程.研究了非Newton流体幂律指数n、主管道毛细数Ca、两相流量比Q、两相黏度比M以及主管道壁面润湿性θ对液滴在T型微通道内的形成尺寸、形成时间和变形参数(DI)的影响.研究结果表明:首先,主管道流体幂律指数n从0.4增加到1.6时,液滴的形成尺寸近似呈线性减小,而液滴的形成时间和变形参数先快速减小,然后缓慢减小;其次,黏度比对液滴形成尺寸、液滴形成以及变形参数的影响与幂律指数的影响基本一致;再者,随着Ca和主管道壁面润湿性的增加,形成液滴的尺寸近似呈线性减小,形成液滴的时间和变形参数先快速减小然后缓慢减小,且减小趋势随幂律指数的增加而减缓;最后,研究结果还表明主管道和子管道的流量比Q越大,液滴形成时间越长,液滴形成尺寸和变形参数越小.     

壁面结构对三维可压缩气泡群影响的数值模拟研究

基于流体体积(VOF)法追踪自由液面,研究了壁面结构对三维可压缩气泡群流动的影响.通过在待测壁面上设置不同形状的壁面结构(长方体、椭球体和圆锥体)并改变它们各自的几何参数(位置和长度),来研究壁面结构对壁面附近的气泡群流动的影响,该影响表现为气泡群对壁面的空间平均压力.研究发现,壁面结构对气泡群的拓扑结构的影响会造成壁...     

Maxwell流体在有限长管道中作不稳定的蠕动传输:食道吞咽进程分析

解析地研究了有限长管道中Maxwell流体的不稳定蠕动传输.管壁受到不超过静止边界的收缩波作用.对无量纲形式的方程,应用长波长近似进行分析.导出了轴向速度和径向速度的表达式,评估了沿波长和管道长度方向的压力.讨论了回流现象,确定了回流极限区域.对食道中咀嚼食物(如面包、蛋白等)传输的数学公式给出了物理上的解释.可以看出,与Newton流体相比,Maxwell流体有利于在食道中的流动.与Takahashi等[Rheology,1999,27:169-172]的实验结果相符合.进一步揭示了松弛时间既不影响剪应力,也不影响回流极限.发现了压力的峰值,对整数值波列是相同的,而对非整数值波列是不同的.     

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

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

多孔饱和矩形管中粘性随温度变化对熵产、热和流体流动的影响

研究了充填流体-饱和多孔介质的矩形管中,随温度变化的粘性对充分发展强迫对流的影响.采用Darcy流动模型并假设粘性-温度为倒线性关系.管壁视为均匀热通量,即Kays和Craw-ford称为的H边界条件.当流体粘性随温度升高而降低时,管壁的Nusselt数增大.求解速度和温度分布时,利用热力学第二定律求解了局部平均熵产率.根据Brinkman数、Péclet数、粘性变化数、无量纲管壁热通量和管道截面宽高比,给出了熵产率、Bejan数、传热不可逆性和流体流动不可逆性的表达式.这些表达式是该类问题参数研究的基础.可以看出,当管道截面宽高比的增大使熵产率减小时,方形管中流动产生的熵大于矩形管,这类似于Ratts和Raut研究的明流(clear flow)情况.     

One-dimensional viscoelastic fluid model where viscosity and normal stress coefficients depend on the shear rate

We study the unsteady motion of a viscoelastic fluid modeled by a second-order fluid where normal stress coefficients and viscosity depend on the shear rate by using a power-law model. To study this problem, we use the one-dimensional nine-director Cosserat theory approach which reduces the exact three-dimensional equations to a system depending only on time and on a single spatial variable. Integrating the equation of conservation of linear momentum over the tube cross-section, with the velocity field approximated by the Cosserat theory, we obtain a one-dimensional system. The velocity field approximation satisfies both the incompressibility condition and the kinematic boundary condition exactly. From this one-dimensional system we obtain the relationship between average pressure and volume flow rate over a finite section of the tube with constant and variable radius. Also, we obtain the correspondent equation for the wall shear stress which enters directly in the formulation as a dependent variable. Attention is focused on some numerical simulation of unsteady/steady flows for average pressure, wall shear stress and on the analysis of perturbed flows.     

Nonlinear mathematical analysis for blood flow in a constricted artery under periodic body acceleration

This study analyses the pulsatile flow of blood through mild stenosed narrow arteries, treating the blood in the core region as a Casson fluid and the plasma in the peripheral layer as a Newtonian fluid. Perturbation method is employed to solve the resulting coupled implicit system of non-linear partial differential equations. The expressions for shear stress, velocity, wall shear stress, plug core radius, flow rate and longitudinal impedance to flow are obtained. The effects of pulsatility, stenosis depth, peripheral layer thickness, body acceleration and non-Newtonian behavior of blood on these flow quantities are discussed. It is noted that the plug core radius, wall shear stress and longitudinal impedance to flow increase as the yield stress and stenosis depth increase and they decrease with the increase of the body acceleration, pressure gradient, width of the peripheral layer thickness. It is observed that the plug flow velocity and flow rate increase with the increase of the pulsatile Reynolds number, body acceleration, pressure gradient and the width of the peripheral layer thickness and the reverse behavior is found when the yield stress, stenosis depth and lead angle increase. It is also recorded that the wall shear stress and longitudinal impedance to flow are considerably lower for the two-fluid Casson model than that of the single-fluid Casson model. It is found that the presence of body acceleration and peripheral layer influences the mean flow rate and mean velocity by increasing their magnitude significantly in the arteries.     

Pulsatile flow of Herschel–Bulkley fluid through catheterized arteries – A mathematical model

The pulsatile flow of blood through catheterized artery has been studied in this paper by modeling blood as Herschel–Bulkley fluid and the catheter and artery as rigid coaxial circular cylinders. The Herschel–Bulkley fluid has two parameters, the yield stress θ and the power index n. Perturbation method is used to solve the resulting quasi-steady nonlinear coupled implicit system of differential equations. The effects of catheterization and non-Newtonian nature of blood on yield plane locations, velocity, flow rate, wall shear stress and longitudinal impedance of the artery are discussed. The existence of two yield plane locations is investigated and their dependence on yield stress θ, amplitude A, and time t are analyzed. The width of the plug core region increases with increasing value of yield stress at any time. The velocity and flow rate decrease, whereas wall shear stress and longitudinal impedance increase for increasing value of yield stress with other parameters held fixed. On the other hand, the velocity, flow rate and wall shear stress decrease but resistance to flow increases as the catheter radius ratio (ratio of catheter radius to vessel radius) increases with other parameters fixed. The results for power law fluid, Newtonian fluid and Bingham fluid are obtained as special cases from this model.     

Mathematical modeling of pulsatile flow of non-Newtonian fluid in stenosed arteries

The pulsatile flow of blood through mild stenosed artery is studied. The effects of pulsatility, stenosis and non-Newtonian behavior of blood, treating the blood as Herschel–Bulkley fluid, are simultaneously considered. A perturbation method is used to analyze the flow. The expressions for the shear stress, velocity, flow rate, wall shear stress, longitudinal impedance and the plug core radius have been obtained. The variations of these flow quantities with different parameters of the fluid have been analyzed. It is found that, the plug core radius, pressure drop and wall shear stress increase with the increase of yield stress or the stenosis height. The velocity and the wall shear stress increase considerably with the increase in the amplitude of the pressure drop. It is clear that for a given value of stenosis height and for the increasing values of the stenosis shape parameter from 3 to 6, there is a sharp increase in the impedance of the flow and also the plots are skewed to the right-hand side. It is observed that the estimates of the increase in the longitudinal impedance increase with the increase of the axial distance or with the increase of the stenosis height. The present study also brings out the effects of asymmetric of the stenosis on the flow quantities.     

弹性动脉血管中血液流动特性的模拟和分析

研究了两个不同的非牛顿血液流动模型:低粘性剪切简单幂律模型和低粘性剪切及粘弹性振荡流的广义Maxwell模型．同时利用这两个非牛顿模型和牛顿模型,研究了磁场中刚性和弹性直血管中血液的正弦型脉动．在生理学条件下,大动脉中血液的弹性对其流动性态似乎并不产生影响,单纯低粘性剪切模型可以逼真地模拟这种血液流动．利用高剪切幂律模型模拟弹性血管中的正弦型脉动流,发现在同一压力梯度下,与牛顿流体相比较,幂律流体的平均流率和流率变化幅度都更小．控制方程用Crank-Niclson方法求解．弹性动脉中血液受磁场作用是产生此结果的直观原因．在主动脉生物流的模拟中,与牛顿流体模型比较,发现在匹配流率曲线上,幂律模型的平均壁面剪切应力增大,峰值壁面剪切应力减小．讨论了弹性血管横切磁场时的血液流动,评估了血管形状和表面不规则等因素的影响．     

Two-fluid Casson model for pulsatile blood flow through stenosed arteries: A theoretical model

Pulsatile flow of blood through mild stenosed narrow arteries is analyzed by treating the blood in the core region as a Casson fluid and the plasma in the peripheral layer as a Newtonian fluid. Perturbation method is used to solve the coupled implicit system of non-linear differential equations. The expressions for velocity, wall shear stress, plug core radius, flow rate and resistance to flow are obtained. The effects of pulsatility, stenosis, peripheral layer and non-Newtonian behavior of blood on these flow quantities are discussed. It is found that the pressure drop, plug core radius, wall shear stress and resistance to flow increase with the increase of the yield stress or stenosis size while all other parameters held constant. The percentage of increase in the resistance to flow over the uniform diameter tube is considerably very low for the present two-fluid model compared with those of the single-fluid model.     

Numerical computation of pulsatile flow through a locally constricted channel

This paper deals with the numerical solution of a pulsatile laminar flow through a locally constricted channel. A finite difference technique has been employed to solve the governing equations. The effects of the flow parameters such as Reynolds number, flow pulsation in terms of Strouhal number, constriction height and length on the flow behaviour have been studied. It is found that the peak value of the wall shear stress has significantly changed with the variation of Reynolds numbers and constriction heights. It is also noted that the Strouhal number and constriction length have little effect on the peak value of the wall shear stress. The flow computation reveals that the peak value of the wall shear stress at maximum flow rate time in pulsatile flow situation is much larger than that due to steady flow. The constriction and the flow pulsation produce flow disturbances at the vicinity of the constriction of the channel in the downstream direction.     

主动脉弓及分支血管内非稳态血流分析

运用流体力学中的三维非定常Navier-Stokes方程作为血液流动的控制方程,并采用计算流体力学方法对人体主动脉弓及分支血管内非Newton(牛顿)血液黏度模型下血流进行瞬态数值模拟.分析了一个心动周期内不同时刻血流动力学特征参数的分布对动脉粥样硬化斑块形成的影响,并与Newton血液黏度模型下的血管壁面压力和壁面切应力特征参数进行对比.结果表明:与Newton血液模型相比,非Newton血液模型下血流分布更符合真实血流特性;在心动收缩期,分支血管外侧壁附近存在面积较大的低速涡流区,该区域内血管壁面压力与壁面切应力具有较大的变化量,血液中的血小板、脂质和纤维蛋白等易沉积,血管内壁易疲劳损伤并发生血管重构,促使动脉粥样硬化斑块形成;而在心动舒张期,分支血管内血流速度分布均匀,血管壁面压力与壁面切应力变化量较小,血管壁受到较小的应力作用,对动脉粥样硬化斑块形成的作用较小.     

Falling film flow of a viscoelastic fluid along a wall

In this paper, we study the heat transfer in the fully developed flow of a viscoelastic fluid, a slag layer, down a vertical wall. A new constitutive relation for the stress tensor of this fluid is proposed, where the viscosity depends on the volume fraction, temperature, and shear rate. For the heat flux vector, we assume the Fourier's law of conduction with a constant thermal conductivity. The model is also capable of exhibiting normal stress effects. The governing equations are non‐dimensionalized and numerically solved to study the effects of various dimensionless parameters on the velocity, temperature, and volume fraction. The effect of the exponent in the Reynolds viscosity model is also discussed. The different cases of shear‐thinning and shear‐thickening, cooling and heating, are compared and discussed. The results indicate that the viscous dissipation and radiation (at the free surface) cause the temperature to be higher inside the flow domain. Copyright © 2013 John Wiley & Sons, Ltd.