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

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

爆震室壁面条件对脉冲爆震发动机爆震性能的影响

在爆震室内快速形成稳定传播的爆轰波是脉冲爆震发动机的关键．本文利用有限速率化学反应模型,考虑粘性、热对流,基于N-S方程对氢气与空气/氧气为反应混合物的爆震发动机爆震室内流场进行计算．从流场压力、速度、涡量、湍流动能等方面研究爆震室壁面条件对燃烧爆轰性能的影响,分析流场爆轰波压力与流场湍动能的关系,讨论可燃气体燃烧转爆轰的机理．结果表明：爆震室内燃烧爆轰机理受到化学反应能量释放、壁面摩擦效应、壁面与外界热交换的影响．在文中讨论的范围内,相比于半圆形和三角形的爆震室装置,矩形的爆震室增强装置能在更短的时间内得到较高的爆轰波压力和湍动能峰值．壁面粗糙层高度（粗糙度）影响爆震室的燃烧爆轰性质．当壁面粗糙度为0.15mm时,粗糙度对爆轰的激励作用大于抑制作用,能较快形成稳定的爆轰波,且推力为35.5N;随着壁面对流换热系数的增大,爆震室壁面的散热加剧．当壁面对流换热系数大于临界值2.6W/(m2·K)时,爆震室内不能形成稳定的爆震波．     

节理岩体的分维特征及其工程地质意义

本文用分形理论,对节理网络分布和粗糙度曲线进行分形分析。结果表明,节理的空间分布和粗糙度曲线具自相似性;分维是一个表征岩体强度、岩体损伤程度、岩体质量及节理粗糙程度的几何参数,并可用于评价岩体稳定性、断裂体系活动性及划分风化卸荷带等。     

无壁面参数低雷诺数非线性涡黏性模式研究

建立了一个低雷诺数的非线性涡黏性湍流模式,该模式的一个显著特性是它不包含壁面参数（如y^ ,n等）,因而特别适用于复杂几何流场的计算,本模式在几种包括回流、分离、激波等典型流动中进行了验证,结果令人满意。     

人工粗糙壁面的水跃特性研究

根据已有文献对密排加糙壁面水跃共轭水深、水跃旋滚长度、水跃长度的试验结果,分析了密排加糙壁面水跃的共轭水深、水跃旋滚长度、水跃长度、壁面平均切应力随弗劳德数、跃前和跃后断面水深、壁面粗糙度的变化规律;给出了人工粗糙壁面水跃共轭水深、水跃旋滚长度、水跃长度、壁面阻力系数、壁面平均切应力的计算公式;通过已有文献的试验结果对公式进行了验证,得到了水跃共轭水深的平均误差为4.06%,水跃旋滚长度和水跃长度的平均误差分别为4.25%和7.16%。研究表明:人工粗糙壁面水跃的共轭水深和水跃长度随着跃前断面弗劳德数的增大而增大,随着壁面粗糙度的增大而减小;壁面平均切应力随着壁面粗糙度和跃前断面弗劳德数的增大而增大,随着共轭水深比的增大而减小。     

Taylor补丁对新型动脉旁路移植流场影响的数值分析

为了研究Taylor补丁对新型(S型)动脉旁路移植术中吻合口处流场的影响,使用数值方法研究了采用Taylor补丁和未采用该补丁的两个S型旁路移植模型内流场的血流动力学差异. 对流速、壁面切应力和切应力梯度等参数进行了比较分析. 结果表明,Taylor补丁对吻合口的流场有显著影响. 采用Taylor补丁的模型其下游吻合口处的流场分布较未采用补丁的模型更均匀,二次流平均流速减小约34.48%,壁面切应力梯度减小约52.22%,从壁面切应力梯度方面分析,这将有助于改善血流动力学分布,抑制动脉粥样硬化. 但从壁面切应力值分析,其动脉底部的壁面低切应力区明显增大,平均壁面切应力值减小30.33%,这又将促使动脉粥样硬化. 因此,Taylor补丁是否对S型搭桥术具有治疗优越性,仅从血流动力学分析尚不能定论,配合数值计算结果进行动物和临床实验研究是十分必要的.      

表面粗糙度对螺旋槽干式气体端面密封性能预测与结构优化的影响

将螺旋槽干式气体端面密封(S-DGS)的粗糙表面分成软环端表面、硬环端面开槽底面与非开槽表面3个区域,建立粗糙表面S-DGS性能的有限元分析理论模型,采用该模型研究3个区域表面粗糙度对S-DGS性能参数和端面几何结构参数优化的影响.结果表明:在相同工况和表面综合均方根差的条件下,粗糙表面S-DGS的开启力、气膜刚度和摩擦扭矩均大于光滑表面S-DGS,而泄漏量较小;在有关标准范围之内,硬环端面槽底面和软环端面的表面粗糙度对S-DGS的性能参数预测值产生较大影响,而硬环端面表面粗糙度的影响可以忽略;表面粗糙度对S-DGS端面几何结构参数的优化值没有影响.研究结果对S-DGS性能的正确预测和结构优化设计以及密封的修复工作具有重要的指导意义.     

表面粗糙度对微细管内气体流动特性的影响

采用了表面粗糙度粘性系数模型对微细管内的气体流动进行数值模拟，以研究微管内壁表面粗糙度对微管内气体流动的影响。运用本文改进的表面粗糙度粘性系数模型，数值模拟与实验数据十分吻合。计算结果表明，进出口压力一定时，表面粗糙度对流场的压力、密度及温度分布的影响不大，但是对速度场的影响十分显著，表面粗糙度使气体流动速度减小，并使壁面附近的速度梯度减小，从而使通过管道的气体质量流量减小，在微管内的气体流动中，表面粗糙度的影响是不能被忽略的。     

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

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

粗糙度对大间隙环流偏心转子动特性系数的影响

基于作者建立的大间隙环流中转子运动的理论模型,用摄动法推导了大间隙环流流场非线性控制方程的一阶摄动方程,采用数值方法研究了静子和转子壁面粗糙度对大间隙环流中偏心转子动特性系数的影响。研究结果表明：静子和转子壁面粗糙度对大间隙环流中偏心转子动特性系数有较大影响。所得到的数值结果与已有的解析解和实验结果具有较好的一致性。     

The effects of surface roughness geometry of flow undergoing Poiseuille flow by molecular dynamics simulation

Numerical simulation of Poiseuille flow of liquid Argon in a rough nano-channel using the non-equilibrium molecular dynamics simulation is performed. Density and velocity profiles across the channel are investigated in which roughness is implemented only on the lower wall. The Lennard–Jones potential is used to model the interactions between all particles. The effects of surface roughness geometry, gap between roughness elements (or roughness periodicity), surface roughness height and surface attraction energy on the behavior of the flow undergoing Poiseuille flow are presented. Results show that surface shape and roughness height have a decisive role on the flow behaviors. In fact, by increasing the roughness ratio (height to base ratio), the slip velocity and the maximum velocity in the channel cross section are reduced, and the density fluctuations near the wall increases. Results also show that the maximum density near the wall for a rough surface is less than a smooth wall. Moreover, the simulation results show that the effect of triangle roughness surface on the flow behavior is more than the cylindrical ones.     

Effect of surface roughness on a turbulent wall-attaching offset jet

The flow characteristics of a two-dimensional offset jet discharged parallel to a rough wall is experimentally investigated by using a split film probe. The distributions of the mean velocity and turbulent stresses in the flow field are obtained and compared with those of the wall attaching offset jet on a smooth wall. It is found that the wall-attaching region on the rough wall is longer than on the smooth wall for the same offset height and the jet speed. The normal distance of the maximum velocity point is farther away from the wall than that for the smooth wall case because of the thick wall boundary layer established by the surface roughness. It is also found that the roughness of the wall accelerates the relaxation process towards redeveloped plane wall jet and that it exhibits a quite different turbulent diffusion behavior especially near the wall from that in the wall jet over a smooth surface.     

A plane turbulent wall jet on a fully rough surface

The mean velocity field and skin friction characteristics of a plane turbulent wall jet on a smooth and a fully rough surface were studied using Particle Image Velocimetry. The Reynolds number based on the slot height and the exit velocity of the jet was Re = 13,400 and the nominal size of the roughness was k = 0.44 mm. For this Reynolds number and size of roughness element, the flow was in the fully rough regime. The surface roughness results in a distinct change in the shape of the mean velocity profile when scaled in outer coordinates, i.e. using the maximum velocity and outer half-width as the relevant velocity and length scales, respectively. Using inner coordinates, the mean velocity in the lower region of the inner layer was consistent with a logarithmic profile which characterizes the overlap region of a turbulent boundary layer; for the rough wall case, the velocity profile was shifted downward due to the enhanced wall shear stress. For the fully rough flow, the decay rate of the maximum velocity of the wall jet is increased, and the skin friction coefficient is much larger than for the smooth wall case. The inner layer is also thicker for the rough wall case. The effects of surface roughness were observed to penetrate into the outer layer and slightly enhance the spread rate for the outer half-width, which was not observed in most other studies of transitionally rough wall jet flows.     

Numerical study of particle deposition in a turbulent channel flow with transverse roughness elements on one wall

A numerical study is presented for the effect of wall roughness on the deposition of solid spherical particles in a fully developed turbulent channel flow based on large eddy simulation combined with a Lagrangian particle-tracking scheme. The interest is focused on particles with response times in wall units in the range of 2.5 ≤ τ_p⁺ ≤ 600 depositing onto a vertical rough surface consisting of two-dimensional transverse square bars separated by a rectangular cavity. Predictions of particle deposition rates are obtained for several values of the cavity width to roughness element height ratio and particle response time. It is shown that the accumulation of particles in the near wall region and their preferential concentration in flow areas of low streamwise fluid velocity that occur in turbulent flows at flat channels are significantly affected by the roughness elements. Particle deposition onto the rough wall is considerably increased, exhibiting a subtle dependence on the particle inertia and the spacing between the bars. The observed augmentation of deposition coefficient can be attributed to the flow modifications induced by the roughness elements and to the inertial impaction of particles onto the frontal deposition area of the protruding square bars.     

Structure of the turbulent boundary layer over a rod-roughened wall

Turbulent coherent structures near a rod-roughened wall are scrutinized by analyzing instantaneous flow fields obtained from direct numerical simulations (DNSs) of a turbulent boundary layer (TBL). The roughness elements used are periodically arranged two-dimensional spanwise rods, and the roughness height is k/δ = 0.05 where δ is the boundary layer thickness. The Reynolds number based on the momentum thickness is varied in the range Re_θ = 300–1400. The effect of surface roughness is examined by comparing the characteristics of the TBLs over smooth and rough walls. Although introduction of roughness elements onto the smooth wall affects the Reynolds stresses throughout the entire boundary layer when scaled by the friction velocity, the roughness has little effect on the vorticity fluctuations in the outer layer. Pressure-strain tensors of the transport equation for the Reynolds stresses and quadrant analysis disclose that the redistribution of turbulent kinetic energy of the rough wall is similar to that of the smooth wall, and that the roughness has little effect on the relative contributions of ejection and sweep motions in the outer layer. To elucidate the modifications of the near-wall vortical structure induced by surface roughness, we used two-point correlations, joint weighted probability density function, and linear stochastic estimation. Finally, we demonstrate the existence of coherent structures in the instantaneous flow field over the rod-roughened surface.     

Effects of two-dimensional V-shaped grooves on turbulent channel flow

Detailed Laser Doppler velocimeter (LDV) measurements have been carried out in a turbulent rectangular channel flow with one rough wall. The roughness elements of two-dimensional spanwise 120° V-shaped grooves are periodically arranged with different depths and pitches. The Reynolds number based on the centerline velocity, and the channel half height ranges from 2,740 to 20,000. The comparisons of turbulence statistics over smooth and rough walls indicate that the present roughness leads to a significant change in the turbulence both in the inner and in the outer flow. Particularly, the distribution density of the grooves is a key parameter to evaluate the effect of roughness. The low-Reynolds-number dependence of turbulence statistics is also observed. The rough walls with the same pitch-to-depth ratio exhibit the equivalent roughness function under the corresponding Reynolds numbers. The disagreement of velocity defect profiles between smooth and rough walls challenges the defect universal law. The variations of the turbulence stresses and Reynolds shear stress decomposition in the outer layer suggest that the turbulent motions may be modified by the present grooves. The importance of sweep events for the present groove-roughened walls is reflected by the differences in relative contribution to Reynolds shear stress from each quadrant and the higher-order moments over smooth and rough walls.     

Numerical observations of turbulence structure modification in channel flow over 2D and 3D rough walls

The effects of wall roughness on turbulence structure modifications were explored by numerical experiments, carried out using Large Eddy Simulation techniques. The wall geometry was made using an archetypal artificial method, thus to achieve irregular two- and three-dimensional shapes. The proposed roughness shapes are highly irregular and are characterised by high and small peaks, thus it can be considered a practical realistic roughness. Their effects are analysed comparing the turbulence quantities over smooth, 2D and 3D rough walls of fully developed channel flow at relatively low friction Reynolds number ${\mathit{Re}}_{\tau }=395$. Both transitional and fully rough regimes have been investigated. The two rough surfaces were built in such a way that the same mean roughness height and averaged mean deviation is obtained. Despite of this, very different quantitative and qualitative results are generated. The analysis of the mean quantitative statistics and turbulence fluctuations shows that deviations are mainly concentrated in the inner layer. These results support the Townsend’s similarity hypothesis. Among the geometrical parameters, which characterise the wall geometries, roughness slope correlates well with the roughness function $\mathrm{\Delta }{U}^{+}$. Specifically, a logarithmic law is proposed to predict the downward shift of the velocity profile for the transitional regime. Instantaneous view of turbulent organised structures display differences in small-scale structures. The flow field over rough surfaces is populated with coherent structures shorter than those observed over flat planes. The comparative analysis of both streaks and wall-normal vortical structures shows that 2D and 3D irregularities have quite different effects. The results highlight that 3D rough wall are representative of a more realistic surface compared to idealised 2D roughness.     

Effects of irregular two-dimensional and three-dimensional surface roughness in turbulent channel flows

Wall-resolved Large Eddy Simulation of fully developed turbulent channel flows over two different rough surfaces is performed to investigate on the effects of irregular 2D and 3D roughness on the turbulence. The two geometries are obtained through the superimposition of sinusoidal functions having random amplitudes and different wave lengths. In the 2D configuration the irregular shape in the longitudinal direction is replicated in the transverse one, while in the 3D case the sinusoidal functions are generated both in streamwise and spanwise directions. Both channel walls are roughened in such a way as to obtain surfaces with statistically equivalent roughness height, but different shapes. In order to compare the turbulence properties over the two rough walls and to analyse the differences with a smooth wall, the simulations are performed at the same Reynolds number Re_τ = 395. The same mean roughness height h = 0.05δ (δ the half channel height) is used for the rough walls.The roughness function obtained with the 3D roughness is larger than in the 2D case, although the two walls share the same mean height. Thus, the considered irregular 3D roughness is more effective in reducing the flow velocity with respect to the 2D roughness, coherently with the literature results that identified a clear dependence of the roughness function on the effective slope (see Napoli et al. (2008)), higher in the generated 3D rough wall. The analysis of higher-order statistics shows that the effects of the roughness, independently on its two- or three-dimensional shape, are mainly confined in the inner region, supporting the Townsend’s wall similarity hypothesis. The tendency towards the isotropization is investigated through the ratio between the resolved Reynolds stress components, putting in light that the 3D irregular rough wall induces an higher reduction of the anisotropy, with respect to the 2D case.     

Görtler vortices promoted by wall roughness

Numerical experiments are conducted to investigate spatially developing Görtler vortices and the way in which wall roughness promotes their formation and growth. Several different types of walls are examined and their relative merits as vortex promoters assessed. The only disturbances of the flow are due to the rough wall; hence, at each downstream station the local field feels (1) the upstream flow distribution (produced by the upstream wall conditions) and (2) the local forcing at the wall. Rapid vortex formation and growth, like in the case of ribleted walls, can be qualitatively explained by the positive combination of these two effects; when the two influences on the local flow field compete, e.g. for randomly distributed wall roughness, the equations with the boundary conditions filter the disturbances over some streamwise length, function of the roughness amplitude, to create coherent patches of vorticity out of the random noise. These patches can then be amplified by the instability mechanism. If a thin rough strip is aligned along the span of an otherwise smooth wall to trip the boundary layer, the filtering region is shorter and growth of the vortices starts earlier. Also for the case of an isolated three-dimensional hump a rapid disturbance amplification is produced, but in this case the vortices remain confined and a very slow spanwise spreading of the perturbation occurs. In all naturally developing cases, where no specific wavelengths are explicity favored, the average spanwise wavelengths computed are very close to those of largest growth from the linear stability theory.