基于滑移速度壁模型的复杂边界湍流大涡模拟

采用滑移速度壁模型实现了浸入边界方法与壁模型相结合的大涡模拟.本文首先分别采用平衡层模型和非平衡壁模型对周期山状流进行数值模拟,以考查在壁模型中考虑切向压力梯度的作用.数值结果表明,流场的压力对本文所采用的壁模型形式并不敏感,但是考虑切向压力梯度可以显著改进壁面摩擦力的计算结果,并且能够准确的预测强压力梯度区以及分离区内的流动平均统计特性.不考虑压力梯度效应的平衡层模型显著低估了壁面摩擦力的分布,同时无法准确预测分离区内的平均速度剖面.非平衡模型的修正项正比于切向压力梯度和壁面法向距离,因此在强压力梯度区或者网格较粗时,计算得到的平均压力和摩擦力分布以及流动的低阶统计量均与参考的实验和计算结果吻合.在此基础上,通过回转体绕流的大涡模拟考查了该方法用于模拟高雷诺数壁湍流的适用性,非平衡壁模型可以准确地捕捉流动的物理结构并较准确地预测其水动力学特性.结果表明,将浸入边界方法与非平衡滑移速度壁模型相结合的大涡模拟,有望成为数值模拟复杂边界高雷诺数壁湍流的工具.      

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

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

Numerical studies of blood flow in healthy,stenosed, and stented carotid arteries

Numerical analysis of pulsatile blood flow in healthy, stenosed, and stented carotid arteries is performed with the aim of identifying hemodynamic factors in the initiation, growth, and the potential of leading to severe occlusions of a diseased artery. The Immersed Finite Element Method is adopted for this study to conveniently incorporate various geometrical shapes of arteries without remeshing. Our computational results provide detailed quantitative analysis on the blood flow pattern, wall shear stress, particle residence time, and oscillatory shear index. The analysis of these parameters leads to a better understanding of blood clot formation and its localization in a stenosed and a stented carotid artery. A healthy artery is also studied to establish a baseline comparison. This analysis will assist in developing treatments for diseased arteries and novel stent designs to reduce restenosis. Copyright © 2008 John Wiley & Sons, Ltd.     

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

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

Effect of couple stresses on the flow in a constricted annulus

The flow of an incompressible couple stress fluid in an annulus with local constriction at the outer wall is considered. This configuration is intended as a simple model for studying blood flow in a stenosed artery when a catheter is inserted into it. The effects couple stress fluid parameters α and σ, height of the constriction (ε), and ratio of radii (k) on the impedance and wall shear stresses are studied graphically. Graphical results show that the resistance to the flow as well as the wall shear stress increases as the ratio of the radii increases and decreases as the couple stress fluid parameters increases.     

Structure of near wall turbulence downstream of a wall mounted protrusion: an interesting Reynolds stress suppression phenomena

A local suppression in the generation of near wall Reynolds stress is achieved by modifying the buffer region and sublayer (y ⁺ <30) of a turbulent pipe flow with a 16.4 wall unit high wall mounted protrusion. Multi-component, multi-point, time resolved laser Doppler velocimetry measurements are made in the undisturbed and modified ARL/PSU glycerin tunnel pipe flow at a Reynolds number of approximately 10000. A downstream converging flow field is produced by the divergence of the approaching mean flow around the protrusion. A pair of counter-rotating vortices, 15 wall units in diameter with common flow down, are generated by the protrusion and also contribute to the wall directed flow convergence. The convergence region is 15 wall units high and more than 100 wall units long and appears to decouple the near wall region from the outer turbulent wall layer. Locally, turbulent velocity fluctuations in the form of Reynolds stress producing events, sweeps and ejections, are retarded within this region. This results in a reduction in near wall uv Reynolds stress and local wall shear. Interestingly, the counter-rotating vortices act to increase turbulent diffusion in a manner which is uncorrelated with Reynolds stress generation.     

The effects of wall roughness on the separated flow over a smoothly contoured ramp

Experimental measurements address the effects on a turbulent boundary layer of wall roughness on a flat plate and a ramp that produces a separation bubble over the ramp trailing edge. A fully rough flow condition is achieved on the upstream flat plate. The main effect of the wall roughness on the outer layer turbulence on a flat plate is to change the friction velocity. The separation region is substantially larger for the rough-wall case. The rough-wall boundary layer turbulence is less sensitive to the onset of an adverse pressure gradient over the ramp, producing substantially smaller Reynolds stress peaks in upstream flat-plate, wall-unit coordinates.     

A reformulated synthetic turbulence generation method for a zonal RANS–LES method and its application to zero-pressure gradient boundary layers

A synthetic turbulence generation (STG) method for subsonic and supersonic flows at low and moderate Reynolds numbers to provide inflow distributions of zonal Reynolds-averaged Navier–Stokes (RANS) – large-eddy simulation (LES) methods is presented. The STG method splits the LES inflow region into three planes where a local velocity signal is decomposed from the turbulent flow properties of the upstream RANS solution. Based on the wall-normal position and the local flow Reynolds number, specific length and velocity scales with different vorticity content are imposed at the inlet plane of the boundary layer. The quality of the STG method for incompressible and compressible zero-pressure gradient boundary layers is shown by comparing the zonal RANS–LES data with pure LES, pure RANS, and direct numerical simulation (DNS) solutions. The distributions of the time and spanwise wall-shear stress, Reynolds stress distributions, and two point correlations of the zonal RANS–LES simulations are smooth in the transition region and in good agreement with the pure LES and reference DNS findings. The STG approach reduces the RANS-to-LES transition length to less than four boundary-layer thicknesses.     

Effect of an adverse pressure gradient on the streamwise Reynolds stress profile maxima in a turbulent boundary layer

It is widely accepted that in a turbulent boundary layer （TBL） with adverse pressure gradient （APG） an outer peak usually appears in the profile of streamwise Reynolds stress. However, the effect of APG on this outer peak is not clearly understood. In this paper, the effect of APG is analysed using the numerical and experimental results in the literature. Because the effect of upstream flow is inherent in the TBL, we first analyse this effect in TBLs with zero pressure gradient on flat plates. Under the individual effect of upstream flow, an outer peak already appears in the profile of streamwise Reynolds stress when the TBL continues developing in the streamwise direction. The APG accelerates the appearance of the outer peak, instead of being a trigger.     

Three-dimensional particle-tracking velocimetry measurement of turbulence statistics and energy budget in a backward-facing step flow

Using a three-dimensional (3-D) particle-tracking velocimeter, detailed turbulent flow measurements were made in a plane channel with a one-sided 50% abrupt expansion, which acted as a backward-facing step. The turbulent channel flow reached a fully developed state well upstream of the step. The Reynolds number based on the upstream centerline velocity and the step height H was 5540. With the mean reattachment point located at 6.51H downstream of the step, the measurement region ranged from −2H upstream to 12H downstream of the step. Various turbulent statistics and the energy budget were calculated from numerous instantaneous vector distributions. As in previous experimental investigations, the Reynolds normal and shear stresses had maximum values upstream of the reattachment. The stress anisotropy tensor revealed a peculiar phenomenon near the reattachment wall, wherein the spanwise normal stress was the largest among the three normal stresses. The triple velocity correlations indicated large values in the separating shear layer, and hence the turbulent diffusion was a major term in the energy budget. Comparison was made between the present results and those of the direct numerical simulation (DNS) of Le et al. (1993), and it was found that the mean and fluctuating velocities, the Reynolds shear stress, and the turbulent energy budget were in excellent agreement, although there was a considerable difference in the inflow conditions.     

An experimental study of turbulent boundary layers approaching separation

The present paper deals with the experimental analysis of a strong decelerated turbulent boundary layer developed on a flat plate. The aim of the study was to examine the effects of pressure gradient on a non-equilibrium boundary layer while indicating local areas of equilibrium flow. The effect of the Reynolds number on a turbulent boundary layer developed with matching the external pressure gradient conditions was also analysed. The emphasis was on the analysis of mean flow statistics i.e. mean velocity profiles, streamwise Reynolds stress and the effect of large- and small-scale interactions by analysing the skewness factor and energy isocontours maps. The comparative analysis of the external data indicated that the structure of the turbulent boundary layer depends not only on local effects of pressure gradient but also on the upstream history of the flow. For the same condition of pressure gradient, the increased momentum is observed near the wall with the increase of the Reynolds number at the Incipient Detachment, where increased turbulence production is also observed, leading to the failure of the outer scaling methods. Surprisingly, the effect of the Reynolds number decays at the intermittent transitory detachment where similar profiles were observed. The upper inflection point in the mean profile corresponded well with the outer maximum of the Reynolds stress and zero crossing of skewness factor. Position of this point occurs at different locations, depending on the flow history effects. The last observation demonstrates that the inflection points results from large- and small-scale interactions, which led to the increased convection velocity of small scales near the wall.     

不同心排出量下主动脉瓣血流动力学的PIV实验研究

主动脉瓣发生病变时导致心排出量(cardiac output, CO)减少,而心排出量减少与主动脉瓣血流动力学耦合作用, 引发瓣膜继发性疾病.本文基于医学影像数据三维重构带有冠状动脉的主动脉根部,制备高度光滑和透明的主动脉根部实验模型, 构建体外脉动循环模拟系统,利用粒子图像测速技术(particle image velocimetry,PIV)研究冠状动脉存在时心排出量对主动脉瓣速度分布、黏性剪应力(viscous shear stress, VSS)和雷诺剪应力(Reynolds shear stress, RSS)等血流动力学的影响.研究结果表明: 冠状动脉的存在改变了主动脉窦中的涡旋运动和涡度,冠状动脉存在时流体经由冠状动脉流出, 主动脉窦中的涡旋运动逐渐消失,涡度较早开始减小. 峰值期, 中心对称流动两侧区域存在正、负高黏性剪切区域,存在冠状动脉一侧的升主动脉下游存在高雷诺剪应力区域.心排出量显著影响主动脉瓣的速度分布、VSS和RSS等血液流动和受力状况.随着心排出量增大, 冠状动脉存在时峰值期的最大速度、VSS和RSS增大, 即$CO=2.1$, 2.8, 3.5和4.2 l/min时, 最大速度分别为0.98, 1.13, 1.21和1.37 m/s, 最大VSS分别为0.87, 0.95, 0.96和1.02 N/m$^{2}$, 最大RSS分别为103.76, 116.25, 138.68和146.55 N/m$^{2}$. 心排出量较低时,主动脉瓣较低的跨瓣流动速度和黏性剪应力易导致血栓形成,研究结果可为主动脉瓣置换术提供理论参考.      

倾斜吹吸控制下湍流边界层减阻的直接数值模拟

雷诺切应力是壁湍流高摩擦阻力的重要来源, 有理论认为可以通过壁面生成负雷诺应力(数值上为正)的方式来削弱湍流流场中雷诺应力的分布, 以此获得流动减阻. 而通过对雷诺平均运动方程的法向二次积分, 可以发现壁面生成正雷诺应力(数值上为负)对壁面摩擦阻力系数才有负贡献. 文中在湍流边界层流动的控制区域下边界设置一系列倾斜狭缝, 利用该装置通过周期性吹吸的方法产生壁面生成正(负)雷诺应力, 并采用直接数值模拟方法考察和验证上文提到的减阻理论. 文中采用的湍流边界层流动模型, 其流动雷诺数(基于外流速度及动量损失厚度)从300 发展到860. 文中通过多组数值模拟算例, 考察了射流强度和频率对壁面摩擦阻力系数的影响, 并对比了壁面生成正或负雷诺应力对流动的影响. 研究表明, 壁面生成正雷诺应力控制的减阻率能达到3.26, 而壁面生成负雷诺应力控制的减阻效果较壁面生成正雷诺应力控制的要差; 壁面生成的正雷诺应力对壁面摩擦阻力有负贡献, 而壁面生成的负雷诺应力对壁面摩擦阻力有正贡献; 通过考察控制的收支比, 发现控制方案不能获得能量净收益.      

Some Reynolds number effects on two- and three-dimensional turbulent boundary layers

Experimental data for a two-dimensional (2-D) turbulent boundary layer (TBL) flow and a three-dimensional (3-D) pressure-driven TBL flow outside of a wing/body junction were obtained for an approach Reynolds number based on momentum thickness of Re _θ =23,200. The wing shape had a 3:2 elliptical nose, NACA 0020 profiled tail, and was mounted on a flat wall. Some Reynolds number effects are examined using fine spatial resolution (Δy ⁺=1.8) three-velocity-component laser-Doppler velocimeter measurements of mean velocities and Reynolds stresses at nine stations for Re _θ =23,200 and previously reported data for a much thinner boundary layer at Re _θ =5,940 for the same wing shape. In the 3-D boundary layers, while the stress profiles vary considerably along the flow due to deceleration, acceleration, and skewing, profiles of the parameter correlate well and over available Reynolds numbers. The measured static pressure variations on the flat wall are similar for the two Reynolds numbers, so the vorticity flux and the measured mean velocities scaled on wall variables agree closely near the wall. The stresses vary similarly for both cases, but with higher values in the outer region of the higher Re _θ case. The outer layer turbulence in the thicker high Reynolds number case behaves similarly to a rapid distortion of the flow, since stream-wise vortical effects from the wall have not diffused completely through the boundary layer at all measurement stations. Received: 9 June 2000/Accepted: 26 January 2001     

Turbulence production in flow over a wavy wall

Measurements of the spatial and time variation of two components of the velocity have been made over a sinusoidal solid wavy boundary with a height to length ratio of 2a/λ = 0.10 and with a dimensionless wave number of α⁺ = (2π/λ)(v/u ^?) = 0.02. For these conditions, both intermittent and time-mean flow reversals are observed near the troughs of the waves. Statistical quantities that are determined are the mean streamwise and normal velocities, the root-meansquare of the fluctuations of the streamwise and normal velocities, and the Reynolds shear stresses. Turbulence production is calculated from these measurements. The flow is characterized by an outer flow and by an inner flow extending to a distance of about α^?1 from the mean level of the surface. Turbulence production in the inner region is fundamentally different from flow over a flat surface in that it is mainly associated with a shear layer that separates from the back of the wave. Flow close to the surface is best described by an interaction between the shear layer and the wall, which produces a retarded zone and a boundary-layer with large wall shear stresses. Measurements of the outer flow compare favorably with measurements over a flat wall if velocities are made dimensionless by a friction velocity defined with a shear stress obtained by extrapolating measurements of the Reynolds stress to the mean levels of the surface (rather than from the drag on the wall).     

PIV EXPERIMENTAL STUDY ON THE HEMODYNAMICS OF AORTIC VALVE UNDER VARIED TILTED ANGLES 1)

瓣叶血栓是主动脉瓣置换术后典型的继发性瓣膜疾病,血流动力学特征异常在其发展过程中至关重要.本文利用粒子图像测速 (particle image velocimetry,PIV) 系统,实验研究了主动脉瓣开口纵向轴线与升主动脉纵向轴线之间倾斜角度 ($\alpha =0^\circ$, $\alpha=5^\circ$,$\alpha =10^\circ$ 和 $\alpha =15^\circ$) 对速度、涡度和黏性剪应力分布等血流动力学特性的影响.研究结果表明：当 $\alpha =0^\circ$ 时,主动脉根部跨瓣血液流动为中心对称流动,而 $\alpha =5^\circ$,$\alpha=10^\circ$ 和 $\alpha =15^\circ$ 时跨瓣血液流动向升主动脉的左冠状动脉一侧倾斜.随着倾斜角度增大,跨瓣血液流动方向倾斜程度增加,血液流动冲击升主动脉壁,损伤内皮细胞导致血栓形成.主动脉瓣倾斜时主动脉窦血液流动速度增大,涡旋也更向主动脉窦底部运动,不利于血液从冠状动脉口流出向心肌供血.同时,主动脉根部的高涡度和高黏性剪应力区域也向升主动脉的左冠状动脉一侧倾斜,主动脉窦的高涡度区域位于主动脉窦底部、高黏性剪应力区域分布于主动脉窦壁面处.主动脉瓣存在倾斜角度时,涡度和黏性剪应力较大,特别是 $\alpha =10^\circ$ 和 $\alpha=15^\circ$,为血栓形成提供了有利环境.研究结果可为临床主动脉瓣置换术参数选择以及继发性瓣膜疾病的避免提供理论依据和技术参考.     

主动脉瓣倾斜角度血流动力学的 PIV 实验研究

瓣叶血栓是主动脉瓣置换术后典型的继发性瓣膜疾病,血流动力学特征异常在其发展过程中至关重要.本文利用粒子图像测速 (particle image velocimetry,PIV) 系统,实验研究了主动脉瓣开口纵向轴线与升主动脉纵向轴线之间倾斜角度 ($\alpha =0^\circ$, $\alpha=5^\circ$,$\alpha =10^\circ$ 和 $\alpha =15^\circ$) 对速度、涡度和黏性剪应力分布等血流动力学特性的影响.研究结果表明:当 $\alpha =0^\circ$ 时,主动脉根部跨瓣血液流动为中心对称流动,而 $\alpha =5^\circ$,$\alpha=10^\circ$ 和 $\alpha =15^\circ$ 时跨瓣血液流动向升主动脉的左冠状动脉一侧倾斜.随着倾斜角度增大,跨瓣血液流动方向倾斜程度增加,血液流动冲击升主动脉壁,损伤内皮细胞导致血栓形成.主动脉瓣倾斜时主动脉窦血液流动速度增大,涡旋也更向主动脉窦底部运动,不利于血液从冠状动脉口流出向心肌供血.同时,主动脉根部的高涡度和高黏性剪应力区域也向升主动脉的左冠状动脉一侧倾斜,主动脉窦的高涡度区域位于主动脉窦底部、高黏性剪应力区域分布于主动脉窦壁面处.主动脉瓣存在倾斜角度时,涡度和黏性剪应力较大,特别是 $\alpha =10^\circ$ 和 $\alpha=15^\circ$,为血栓形成提供了有利环境.研究结果可为临床主动脉瓣置换术参数选择以及继发性瓣膜疾病的避免提供理论依据和技术参考.      

血管残余应力的一种确定方法

动脉无载荷状态下存在残余应力,分析血管壁在载荷状态下的应力分布必须计及血管壁的残余应力。本文通过检测动脉不同轴向伸长比条件下的压力—容积(p—V)数据,根据其p—V实验曲线的“S”型特征,提出了一种用反正切函数关系对试验数据点进行拟合,进而导得动脉管壁周向残余应力沿壁厚分布的方法。文章对10只正常大鼠颈动脉的分析结果表明,大鼠颈动脉周向残余应力沿血管壁厚呈单调上升趋势。残余应力在内壁处为负呈压缩状态,在外壁处为正呈拉伸状态,并且内外壁处应力绝对值大小基本相等,大约在2．5kPa左右。文章结果对进一步分析动脉管壁在载荷作用下的应力分布提供必要的信息。     

BLOOD FLOW AND MACROMOLECULAR TRANSPORT IN CURVED BLOOD VESSELS

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.