冠状动脉对称双路移植管中脉动血流的数值模拟

冠状动脉旁路移植管搭桥术后，常会产生血管再狭窄，导致手术失败，。这与移植管的几何结构及血流动力学是密切相关的。作为改进措施，作者提出了采用对称双路搭桥的设想。本文利用有限元分析方法，对冠状动脉搭桥术中对称双路移植管内的生理流动进行了数值仿真。计算了缝合区附近的流场、壁面切应力、压力等血流动力学因素在心动周期内的时空分布情况。计算结果表明，对称双路搭桥具有较好的血流动力学，可以改善血管流场状况和减轻再狭窄发生。这对临床手术计划是很有帮助和指导意义的。     

Numerical simulation of the haemodynamics in end‐to‐side anastomoses

In this paper, the commercial CFD package Ansys workbench 11 was used to analyse the three‐dimensional haemodynamics of a typical stenotic Coronary Artery Bypass Grafting (CABG). Two end‐to‐side CABG configurations with anastomosis angle of 20^° and 40^° and graft–artery diameter ratios of 1/0.6, 1/1 and 1/1.6 were examined. The flow measurements from in vitro Doppler guide wire technique acquired in left interior mammary artery (LIMA) and grafted to the left anterior descending artery (LAD) were used to impose the physiologically flow conditions at proximal and distal CABG inlet and outlet, respectively. The blood flow was considered to be incompressible, pulsatile, Newtonian, and laminar rheology. The main objective was to determine the effect of anastomosis angle and graft–artery diameter ratio on the flow patterns and the long‐term functionality of the graft. In analysing the results, the distributions of temporal and spatial wall shear stress (WSS) gradient and oscillating shear index (OSI) in the critical regions of CABG such as heel, toe and the centre of the junction were presented and the vortex motions and the occurrence of recirculation zones were examined. The findings showed asymmetrically disturbed flows in the localized regions of the proximal and distal host artery for all models considered and the movement of the recirculation zones from heel to toe was found to depend on the time at the cardiac cycle. These regions are known as susceptible sites to thrombosis and re‐stenosis due to their association with low values of WSS. Copyright © 2010 John Wiley & Sons, Ltd.     

Advantageous swirling flow in 45° end-to-side anastomosis

The effects of swirling flow on the flow field in 45° end-to-side anastomosis are experimentally investigated using a particle image velocimetry technique to reveal fluid dynamic advantages of swirling flow in the vascular graft. Non-swirling Poiseuille inlet flow unnecessarily induces pathological hemodynamic features, such as high wall shear stress (WSS) at the ‘bed’ side and large flow separation at the ‘toe’ side. The introduction of swirling flow is found to equalize the asymmetric WSS distribution and reduces the peak magnitude of WSS. In particular, the intermediate swirling intensity of S = 0.45 induces the most uniform axial velocity and WSS distributions compared with weaker or stronger swirling flows, which addresses the importance of proper selection of swirling intensity in the vascular graft to obtain optimum flow fields at the host vessel. In addition, swirling flow reduces the size of flow separation because it disturbs the formation of Dean-type vortices in secondary flow and inhibits secondary flow collision. The beneficial fluid dynamic features of swirling flow obtained in this study are helpful for designing better vascular graft suppressing pathological hemodynamic features in the recipient host vessel.     

Multiphase non‐Newtonian effects on pulsatile hemodynamics in a coronary artery

Hemodynamic stresses are involved in the development and progression of vascular diseases. This study investigates the influence of mechanical factors on the hemodynamics of the curved coronary artery in an attempt to identify critical factors of non‐Newtonian models. Multiphase non‐Newtonian fluid simulations of pulsatile flow were performed and compared with the standard Newtonian fluid models. Different inlet hematocrit levels were used with the simulations to analyze the relationship that hematocrit levels have with red blood cell (RBC) viscosity, shear stress, velocity, and secondary flow. Our results demonstrated that high hematocrit levels induce secondary flow on the inside curvature of the vessel. In addition, RBC viscosity and wall shear stress (WSS) vary as a function of hematocrit level. Low WSS was found to be associated with areas of high hematocrit. These results describe how RBCs interact with the curvature of artery walls. It is concluded that although all models have a good approximation in blood behavior, the multiphase non‐Newtonian viscosity model is optimal to demonstrate effects of changes in hematocrit. They provide a better stimulation of realistic blood flow analysis. Copyright © 2008 John Wiley & Sons, Ltd.     

不同心排出量下主动脉瓣血流动力学的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}$. 心排出量较低时,主动脉瓣较低的跨瓣流动速度和黏性剪应力易导致血栓形成,研究结果可为主动脉瓣置换术提供理论参考.      

An improved near‐wall modeling for large‐eddy simulation using immersed boundary methods

An improved near‐wall modeling for large‐eddy simulation using the immersed boundary method is proposed. It is shown in this study that the existing near‐wall modeling for the immersed boundary (IB) methods that imposes the velocity boundary condition at the IB node is not sufficient to enforce a correct wall shear stress at the IB node. A new method that imposes a shear stress condition through the modification of the subgrid scale‐eddy viscosity at the IB node is proposed. In this method, the subgrid eddy viscosity at the IB node is modified such that the viscous flux at the face adjacent to the IB node correctly approximates the total shear stress. The method is applied to simulate the fully developed turbulent flows in a plane channel and a circular pipe. It is demonstrated that the new method improves the prediction of the mean velocity and turbulence stresses in comparison with the existing wall modeling based solely on the velocity boundary condition. Copyright © 2015 John Wiley & Sons, Ltd.     

A σ‐coordinate non‐hydrostatic model with embedded Boussinesq‐type‐like equations for modeling deep‐water waves

A σ‐coordinate non‐hydrostatic model, combined with the embedded Boussinesq‐type‐like equations, a reference velocity, and an adapted top‐layer control, is developed to study the evolution of deep‐water waves. The advantage of using the Boussinesq‐type‐like equations with the reference velocity is to provide an analytical‐based non‐hydrostatic pressure distribution at the top‐layer and to optimize wave dispersion property. The σ‐based non‐hydrostatic model naturally tackles the so‐called overshooting issue in the case of non‐linear steep waves. Efficiency and accuracy of this non‐hydrostatic model in terms of wave dispersion and nonlinearity are critically examined. Overall results show that the newly developed model using a few layers is capable of resolving the evolution of non‐linear deep‐water wave groups. Copyright © 2009 John Wiley & Sons, Ltd.     

Dynamic scaling of unsteady shear-thinning non-Newtonian fluid flows in a large-scale model of a distal anastomosis

Dimensional analysis has been applied to an unsteady pulsatile flow of a shear-thinning power-law non-Newtonian liquid. An experiment was then designed in which both Newtonian and non-Newtonian liquids were used to model blood flow through a large-scale (38.5 mm dia.), simplified, rigid arterial junction (a distal anastomosis of a femorodistal bypass). The flow field within the junction was obtained by Particle Imaging Velocimetry and near-wall velocities were used to calculate the wall shear stresses. Dimensionless wall shear stresses were obtained at different points in the cardiac cycle for two different but dynamically similar non-Newtonian fluids; the good agreement between the measured dimensionless wall shear stresses confirm the validity of the dimensional analysis. However, blood exhibits a constant viscosity at high-shear rates and to obtain complete dynamic similarity between large-scale experiments and life-scale flows, the high-shear viscosity also needs to be included in the analysis. How this might be done is discussed in the paper.     

Interference effect of two equal‐sized square cylinders in tandem arrangement: With planar shear flow

Numerical simulations have been performed for flow past two equal‐sized square cylinders in tandem arrangement subjected to incoming planar shear flow. Effect of L/d ratio and the shear parameter has been studied. The range of L/d ratio (ratio of center‐to‐center distance (L) to cylinder width (d)) is varied from 2 to 7 and the non‐dimensional shear parameter (K) is varied from 0.0 to 0.4 in steps of 0.1. For all the cases the Reynolds number (Re) based on centerline velocity and cylinder width is fixed at 100. The results are compared with that of isolated square cylinder with uniform flow. Strouhal number decreases with increasing shear parameter. There are more than one shedding frequency at high shear parameters and L/d ratios. The mean drag coefficient is decreased with shear parameter and lesser than that of the single cylinder. The root mean square (RMS) value of both lift and drag coefficients is higher for the downstream cylinder for all values of shear parameter. With increasing L/d ratio, for both lift and drag, the RMS value increases and then decreases for upstream cylinder, whereas it continuously increases for the downstream cylinder. The stagnation point is moved towards the top leading edge with increasing shear. The critical L/d ratio, which is defined as the distance between two cylinders, beyond which the vortex shedding from the upstream cylinder occurs, decreases with increasing shear parameter. Copyright © 2007 John Wiley & Sons, Ltd.     

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

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

Non‐Newtonian blood flow dynamics in a right internal carotid artery with a saccular aneurysm

Flow dynamics plays an important role in the pathogenesis and treatment of cerebral aneurysms. The temporal and spatial variations of wall shear stress in the aneurysm are hypothesized to be correlated with its growth and rupture. In addition, the assessment of the velocity field in the aneurysm dome and neck is important for the correct placement of endovascular coils. This work describes the flow dynamics in a patient‐specific model of carotid artery with a saccular aneurysm under Newtonian and non‐Newtonian fluid assumptions. The model was obtained from three‐dimensional rotational angiography image data and blood flow dynamics was studied under physiologically representative waveform of inflow. The three‐dimensional continuity and momentum equations for incompressible and unsteady laminar flow were solved with a commercial software using non‐structured fine grid with 283 115 tetrahedral elements. The intra‐aneurysmal flow shows complex vortex structure that change during one pulsatile cycle. The effect of the non‐Newtonian properties of blood on the wall shear stress was important only in the arterial regions with high velocity gradients, on the aneurysmal wall the predictions with the Newtonian and non‐Newtonian blood models were similar. Copyright © 2005 John Wiley & Sons, Ltd.     

A locally DIV‐free projection scheme for incompressible flows based on non‐conforming finite elements

We present a projection scheme whose end‐of‐step velocity is locally pointwise divergence free, using a continuous ?₁ approximation for the velocity in the momentum equation, a first‐order Crouzeix–Raviart approximation at the projection step, and a ?₀ approximation for the pressure in both steps. The analysis of the scheme is done only for grids that guarantee the existence of a divergence free conforming ?₁ interpolant for the velocity. Optimal estimates for the velocity error in L²‐ and H¹‐norms are deduced. The numerical results demonstrate that these estimates should also hold on grids on which the continuous ?₁ approximation for the velocity locks. Since the end‐of‐step velocity is locally solenoidal, the scheme is recommendable for problems requiring good mass conservation. Copyright © 2005 John Wiley & Sons, Ltd.     

A non‐linear dynamical system approach to finite amplitude Taylor‐Vortex flow of shear‐thinning fluids

The effect of shear thinning on the stability of the Taylor–Couette flow is explored for a Carreau–Bird fluid in the narrow‐gap limit. The Galerkin projection method is used to derive a low‐order dynamical system from the conservation of mass and momentum equations. In comparison with the Newtonian system, the present equations include additional non‐linear coupling in the velocity components through the viscosity. It is found that the critical Taylor number, corresponding to the loss of stability of the circular Couette flow, becomes lower as the shear‐thinning effect increases. That is, shear thinning tends to precipitate the onset of Taylor vortex flow, which coincides with the onset of a supercritical bifurcation. Comparison with existing measurements of the effect of shear thinning on the critical Taylor and wave numbers show good agreement. The Taylor vortex cellular structure loses its stability in turn, as the Taylor number reaches a critical value. At this point, an inverse Hopf bifurcation emerges. In contrast to Newtonian flow, the bifurcation diagrams exhibit a turning point that sharpens with shear‐thinning effect. Copyright © 2004 John Wiley & Sons, Ltd.     

Solid‐particle erosion in the tube end of the tube sheet of a shell‐and‐tube heat exchanger

Erosion is one of the major problems in many industrial processes, and in particular, in heat exchangers. The effects of flow velocity and sand particle size on the rate of erosion in a typical shell‐and‐tube heat exchanger were investigated numerically using the Lagrangian particle‐tracking method. Erosion and penetration rates were obtained for sand particles of diameters ranging from 10 to 500 µm and for inlet flow velocities ranging from 0.197 to 2.95 m/s. A flow visualization experiment was conducted with the objective of verifying the accuracy of the continuous phase calculation procedure. Comparison with available experimental data of penetration rates was also conducted. These comparisons resulted in a good agreement. The results show that the location and number of eroded tubes depend mainly on the particle size and velocity magnitude at the header inlet. The rate of erosion depends exponentially on the velocity. The particle size shows negligible effect on the erosion rate at high velocity values and the large‐size particles show less erosion rates compared to the small‐size particles at low values of inlet flow velocities. The results indicated that the erosion and penetration rates are insignificant at the lower end of the velocity range. However, these rates were found to increase continuously with the increase of the inlet flow velocity for all particle sizes. The particle size creating the highest erosion rate was found to depend on the flow velocity range. Copyright © 2005 John Wiley & Sons, Ltd.     

Some features of a turbulent wing–body junction vortical flow

Laser-Doppler velocimeter measurements of a wing/body junction flow field made within a plane to the side of the wing/wall junction and perpendicular both to a 3:2 elliptical nose—NACA 0020 tail wing, and a flat wall are presented. Reynolds number of the approach boundary layer was, Re_θ = 5940, and free-stream air velocity was, U_ref = 27.5 m/s. A large vortical structure residing in the outer region redirects the low-turbulence free-stream flow to the vicinity of the wing/wall junction, resulting in thin boundary layers with velocity magnitudes higher than free-stream flow. Lateral pressure gradients result in a three-dimensional separation on the uplifting side of the vortex. Additionally, a high vorticity vortical structure with opposite sense to the outer-layer vortex forms beneath the outer-layer vortex. Normal and shear stresses increase to attain values an order of magnitude larger compared to values measured in a three-dimensional boundary layer just outside the junction vortex. Bimodal histograms of the w fluctuating velocity occur under the outer-layer vortex near the wall due to the time-dependent nature of the horseshoe vortex. In such a flow the shear-stress angle (SSA) highly lags the flow-gradient angle (FGA), and the turbulence diffusion is highly altered due to presence of vortical structures.     

A two‐dimensional vertical non‐hydrostatic σ model with an implicit method for free‐surface flows

An implicit finite difference model in the σ co‐ordinate system is developed for non‐hydrostatic, two‐dimensional vertical plane free‐surface flows. To accurately simulate interaction of free‐surface flows with uneven bottoms, the unsteady Navier–Stokes equations and the free‐surface boundary condition are solved simultaneously in a regular transformed σ domain using a fully implicit method in two steps. First, the vertical velocity and pressure are expressed as functions of horizontal velocity. Second, substituting these relationship into the horizontal momentum equation provides a block tri‐diagonal matrix system with the unknown of horizontal velocity, which can be solved by a direct matrix solver without iteration. A new treatment of non‐hydrostatic pressure condition at the top‐layer cell is developed and found to be important for resolving the phase of wave propagation. Additional terms introduced by the σ co‐ordinate transformation are discretized appropriately in order to obtain accurate and stable numerical results. The developed model has been validated by several tests involving free‐surface flows with strong vertical accelerations and non‐linear waves interacting with uneven bottoms. Comparisons among numerical results, analytical solutions and experimental data show the capability of the model to simulate free‐surface flow problems. Copyright © 2004 John Wiley & Sons, Ltd.     

An efficient and accurate non‐hydrostatic model with embedded Boussinesq‐type like equations for surface wave modeling

A novel approach that embeds the Boussinesq‐type like equations into an implicit non‐hydrostatic model (NHM) is developed. Instead of using an integration approach, Boussinesq‐type like equations with a reference velocity under a virtual grid system are introduced to analytically obtain an analytical form of pressure distribution at the top layer. To determine the size of vertical layers in the model, a top‐layer control technique is proposed and the reference location is employed to optimize linear wave dispersion property. The efficiency and accuracy of this NHM with Boussinesq‐type like equations (NHM‐BTE) are critically examined through four free‐surface wave examples. Overall model results show that NHM‐BTE using only two vertical layers is capable of accurately simulating highly dispersive wave motion and wave transformation over irregular bathymetry. Copyright © 2008 John Wiley & Sons, Ltd.     

Simulations of high‐aspect‐ratio jets

There are many practical situations when jets are emanating from non‐axis‐symmetric apertures, yet numerical simulations of such three‐dimensional jets are scarce and most of them have failed to reproduce some of the unique flow features. Examples of this type of jets are gas leaks from flanges. These can be treated as jets issuing from high aspect ratio rectangular orifices. The present work consists of a series of large eddy simulations typifying such jets using different inflow boundary conditions. Good agreement with available experiments was observed provided appropriate boundary conditions were present. A discrete method for formulating turbulence data with a known energy spectrum for an inflow condition is outlined and evaluated with three other inflow conditions–a steady uniform profile, a steady parabolic profile and pseudo‐random noise. The implementation of the new inlet condition results in a more realistic centreline velocity decay where the division between the end of the potential core region and the start of the characteristic decay region is clearly visible. Large velocity oscillations are also observed in the final quarter of the domain (15–20 slot widths downstream). Similar oscillations have been observed in real jets. Off‐centre mean velocity peaks are present along the major axis 10 slot widths downstream of the exit in all the simulations. The peaks are approximately 3% of the centreline velocity. The presence of the off‐centre peaks are proved to be independent of jet inflow boundary conditions and an explanation for the mechanism causing the off‐centre peaks is given. Copyright © 2002 John Wiley & Sons, Ltd.     

Numerical simulation of two‐dimensional laminar slot‐jet impingement flows confined by a parallel wall

Two‐dimensional laminar incompressible impinging slot‐jet is simulated numerically to gain insight into flow characteristics.Computations are done for vertically downward‐directed slot‐jets impinging on a plate at the bottom and confined by a parallel surface on top. The behaviour of the jet with respect to aspect ratio (AR) and Reynolds number (Re) are described in detail. The computed flow patterns for various AR (2–5) and for a range of jet‐exit Reynolds numbers (100–500) are analysed to understand the flow characteristics. The transient development of the flow is also simulated for AR = 4 and Re = 300. It is found that the reattachment length is dependent on both AR and Reynolds number for the range considered. The correlation for reattachment length is suggested. The maximum resultant velocity V_rmax and its trajectory is reported. A detailed study of horizontal velocity profile at different downstream locations is reported. It is found that the effect of Reynolds number and AR is significant to the bottom wall vorticity in the impingement and wall jet regions. Copyright © 2007 John Wiley & Sons, Ltd.     

Generation of turbulent inflow and initial conditions based on multi‐correlated random fields

Direct or large eddy simulation of a turbulent flow field is strongly influenced by its initial or inflow boundary condition. This paper presents a new stochastic approach to generate an artificial turbulent velocity field for initial or inflow boundary condition based on digital filtering. Each velocity component of the artificial turbulent velocity field is generated by linear combination of individual uncorrelated random fields. These uncorrelated random fields are obtained by filtering random white‐noise fields. Using common elements in these linear combinations results in multi‐correlation among different velocity components. The generated velocity field reproduces locally desired Reynolds stress components and integral length scales including cross‐integral length scales. The method appears to be simple, flexible and more accurate in comparison with previously developed methods. The accuracy and performance of the method are demonstrated by numerical simulation of a homogeneous turbulent shear flow with high and low shear rates. To assess the accuracy and performance of the method, simulation results are compared with a reference simulation. Copyright © 2007 John Wiley & Sons, Ltd.