主动脉瓣倾斜角度血流动力学的 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$,为血栓形成提供了有利环境.研究结果可为临床主动脉瓣置换术参数选择以及继发性瓣膜疾病的避免提供理论依据和技术参考.      

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$,为血栓形成提供了有利环境.研究结果可为临床主动脉瓣置换术参数选择以及继发性瓣膜疾病的避免提供理论依据和技术参考.     

壁湍流多尺度相干结构复涡黏模型的实验研究

在湍流相干结构动力学方程中,非相干结构成分对相干结构贡献的雷诺应力的模型为涡黏性模型,即涡黏性系数乘以相干结构平均速度变形率的形式.基于非相干结构成分对相干结构贡献的雷诺应力与相干结构速度变形率之间存在相位差的事实,在理论上提出了非相干结构成分对相干结构贡献的雷诺应力复涡黏性模型的假设.应用热线测速技术,在低速风洞中对湍流边界层非相干结构成分对相干结构贡献的雷诺应力与相干结构法向速度变形率之间的相位关系进行了实验测量.通过分析湍流相干结构猝发过程中非相干结构成分对相干结构贡献的雷诺应力与相干结构速度变形率之间的相位关系,研究了相干结构雷诺应力分量与流向速度法向梯度之间的相位差沿湍流边界层法向的变化规律,肯定了湍流相干结构复涡黏性系数模型的合理性.     

Experimental research on complex eddy viscosity modeling of multi-scale coherent structures in wall turbulence

在湍流相干结构动力学方程中,非相干结构成分对相干结构贡献的雷诺应力的模型为涡黏性 模型,即涡黏性系数乘以相干结构平均速度变形率的形式. 基于非相干结构成分对相干结构贡 献的雷诺应力与相干结构速度变形率之间存在相位差的事实,在理论上提出了非相干结构成 分对相干结构贡献的雷诺应力复涡黏性模型的假设. 应用热线测速技术,在低速风洞中对湍 流边界层非相干结构成分对相干结构贡献的雷诺应力与相干结构法向速度变形率之间的相位 关系进行了实验测量. 通过分析湍流相干结构猝发过程中非相干结构成分对相干结构贡献的 雷诺应力与相干结构速度变形率之间的相位关系,研究了相干结构雷诺应力分量与流向速度 法向梯度之间的相位差沿湍流边界层法向的变化规律,肯定了湍流相干结构复涡黏性系数模 型的合理性.     

应用PIV技术测试涡旋波流场

涡旋波流动作为一种特殊的流动现象,可以使流体在相对较宽的槽道中产生较强的波动和对流混合,从而在小Re数条件下起到强化传质的效果。本文利用PIV流场显示技术,对振荡流在非对称槽道中所形成的涡旋波的产生机理和发展规律进行了实验研究和定量分析,测得了涡旋波流场的速度矢量图,阐明了涡旋波流场周期性变化的特点。分析了Re数和St数对涡旋波流动的影响,并得出了旋涡涡心位置以及涡心处涡量的动态变化规律。     

k-ω SST湍流模式三维激波分离流动修正

“激波?边界层分离”是航空气动领域的典型湍流非平衡流动问题, 准确模拟激波分离对于跨声速飞行器气动性能评估和优化设计具有重要意义. 然而传统涡黏性湍流模式中涡黏性系数的定义方式并不适用于非平衡流动, k-ω SST湍流模式为此引入的Bradshaw假设在应用于三维强逆压梯度和较大分离流动时反而限制了雷诺应力的生成, 导致包括k-ω SST在内的常用涡黏性湍流模式均无法对此类流动进行准确模拟. 同时, 现有的非线性雷诺应力本构关系也并不能有效提高模拟精度. 为此, 针对k-ω SST模式分别提出了基于Bradshaw假设和基于长度尺度的两种激波分离流动修正方法. 前者通过提高Bradshaw常数的方式放宽了对雷诺应力生成的限制, 后者则从湍流长度尺度概念出发, 利用混合长度理论、湍动能生成/耗散之比和一种新定义的长度尺度之比构造了ω方程耗散项修正函数, 提高了模式在三维激波分离流动中的建模长度尺度. 两种方法对ONERA M6机翼跨声速大攻角流动均能得到较雷诺应力模式更好的模拟结果. 进一步的雷诺应力分析表明, 三维激波分离流动中“主雷诺应力分量”的概念不再成立, 各雷诺应力分量大小接近. 网格收敛性分析、对其他攻角状态的验证以及湍流平板边界层壁面律验证进一步确认了所提出的两种修正方法的合理性、有效性和通用性.      

求解复杂湍流的非线性涡黏性系数模型和代数应力模型

非线性涡黏性系数模型和代数应力模型联系了线性涡黏性系数湍流模型和完整的微分 雷诺应力模型.随着它们受到日益关注,其形式也越来越多样化.本篇综述的目的是对这些模 型加以总结并比较它们之间的共同点及不同之处,指出它们与完整微分雷诺应力模型之间的 关系,以及相对于线性涡黏性系数模型而言它们在预报流场上所具有的优势.     

非对称槽道中涡旋波的特性研究

利用PIV流场显示技术，对振荡流体在非对称槽道中涡旋波的产生、发展和消失的规律进 行了实验研究和分析，测得了涡旋波流场的速度矢量图，阐明了涡旋波流场周期性变化的特 点. 结合涡动力学方程，深入分析并揭示了做周期性运动的流体能在槽道中产生波的特性这 一规律，从中发现：流体周期变化的非定常性和不对称的槽道结构是形成涡旋波流动的主要 因素. 本文对涡旋波流场中各个旋涡的速度分布和涡量进行了测量和计算，分析了涡旋波 强化传质的机理，并研究了Re数对涡旋波流动的影响     

对二维分离流涡黏性系数非线性分布的新认识

以弱非线性涡黏性模型为出发点,对Delery分离流动实验结果进行分析并获得了非平衡态分离区涡黏性系数与形状因子J之间的非线性关系. 该非线性关系显示在分离起始阶段,涡黏性系数较平衡态先减小,后增大;再附阶段,涡黏性系数较平衡态数值逐渐增大,并在再附点位置接近最大,而后又逐渐减小,恢复到平衡态水平. 总结涡黏性系数的这种非线性发展数学关系式,并将它应用于BL模型,在不添加微分方程的情况下发展出一种适用于分离流动的改进代数湍流模型. 对低速平板流动,跨声速,超声速以及高超声速分离流动的计算结果表明,该改进湍流模型可以较准确地模拟各类复杂分离流动,计算精度明显优于传统代数模型以及一些两方程模型,而计算工作量仍与BL模型相当. 这表明所提出的涡黏性系数非线性发展规律是正确的,且应用在二维分离流动中具有一定的普适性.      

基于涡旋识别方法的高速列车尾涡结构的讨论

利用改进型延迟分离涡模拟方法对缩尺比例1:30的高速列车简化模型的绕流流场进行数值计算,主要针对近尾流区的涡旋结构展开具体讨论. 通过不同的涡旋识别方法,发现在尾涡结构中,高涡量的强涡旋主要聚集于尾车附近,而涡量较低但处于相对稳定状态的涡旋分布在大部分尾流空间中. 对此,主要基于最新提出的涡旋定义及其物理意义认为,由于边界层在尾部发生的流动分离,剪切变形以及高涡量的扩散对强涡旋的形成发挥着重要的作用,而涡旋会被较强的剪切旋转拉伸,使得局部复杂的流动表现出突出的湍流特性;另一方面,尽管涡强度明显下降,但是在强剪切应变迅速衰减的情况下,流向涡核中的涡旋涡量是主要的,此时,在较接近地面的情况下,流体微团以涡核为中心的旋转运动使得涡旋与地面之间的相互作用成为主导的流动机制. 虽然涡强度会相对缓慢地衰减,但是从湍流能量产生的角度,该机制对涡旋的自维持发挥重要的作用,从而使尾涡结构能够相对稳定地存在于尾流流动中.      

基于等温吸附的页岩水分传输特征研究

研究页岩的水分传输特征至关重要,不仅有助于认识页岩的物理化学性质,而且也有助于评价页岩气的吸附扩散和流动能力.本文设计了页岩的水分传输实验装置,采用美国伍德福德和中国南方龙马溪组页岩为研究对象,开展了不同温度、不同湿度下页岩的水分传输实验,研究了页岩的水分传输特征和影响因素.结果表明,页岩的水分吸附属于II型曲线,包含着单分子层吸附、多分子层吸附和毛细凝聚的过程,GAB模型可用于描述页岩的水分吸附过程;水分吸附随着相对压力的增大而增强,有机碳含量和温度对页岩水分吸附起着增强作用,而方解石会抑制页岩的水分吸附;随着相对压力的增大,页岩的水分扩散系数呈现先增大后减小随后增加的趋势,其系数大约在8.73$\times$10$^{ - 9}\sim $5.95$\times $10$^{ - 8 }$m$^{2}$/s之间;伍德福德页岩的等量吸附热均大于龙马溪页岩的等量吸附热,这与其页岩的成熟度有关.研究结果为认识页岩的物理化学性质和力学性能以及评价页岩气的吸附流动能力提供参考依据.      

Spatiotemporal complexity of the aortic sinus vortex

The aortic sinus vortex is a classical flow structure of significant importance to aortic valve dynamics and the initiation and progression of calcific aortic valve disease. We characterize the spatiotemporal characteristics of aortic sinus vortex dynamics in relation to the viscosity of blood analog solution as well as heart rate. High-resolution time-resolved (2 kHz) particle image velocimetry was conducted to capture 2D particle streak videos and 2D instantaneous velocity and streamlines along the sinus midplane using a physiological but rigid aorta model fitted with a porcine bioprosthetic heart valve. Blood analog fluids used include a water–glycerin mixture and saline to elucidate the sensitivity of vortex dynamics to viscosity. Experiments were conducted to record 10 heart beats for each combination of blood analog and heart rate condition. Results show that the topological characteristics of the velocity field vary in timescales as revealed using time bin-averaged vectors and corresponding instantaneous streamlines. There exist small timescale vortices and a large timescale main vortex. A key flow structure observed is the counter vortex at the upstream end of the sinus adjacent to the base (lower half) of the leaflet. The spatiotemporal complexity of vortex dynamics is shown to be profoundly influenced by strong leaflet flutter during systole with a peak frequency of 200 Hz and peak amplitude of 4 mm observed in the saline case. While fluid viscosity influences the length and timescales as well as the introduction of leaflet flutter, heart rate influences the formation of counter vortex at the upstream end of the sinus. Higher heart rates are shown to reduce the strength of the counter vortex that can greatly influence the directionality and strength of shear stresses along the base of the leaflet. This study demonstrates the impact of heart rate and blood analog viscosity on aortic sinus hemodynamics.     

Three-dimensional vorticity in a turbulent cylinder wake

The three components of the vorticity vector in the intermediate region of a turbulent cylinder wake were measured simultaneously using a multi-hot-wire probe. This probe has an improved spatial resolution compared with those reported in the literature. The behavior of the instantaneous velocity and vorticity signals is examined. Both coherent and incoherent vorticity fields are investigated using a phase-averaged technique. The iso-contours of the phase-averaged longitudinal and lateral vorticity variances, and , wrap around the spanwise structures of opposite sign and run through the saddle point along the diverging separatrix. The observation conforms to the previous reports of the occurrence of the longitudinal structures based on flow visualizations and numerical simulations. The magnitude of these contours is about the same as that of the maximum coherent spanwise vorticity at the vortex center, indicating that the strength of the longitudinal structures is comparable to that of the spanwise vortices. Furthermore, and exhibit maximum concentration away from the vortex center, probably because of a combined effect of the large-scale spanwise vortices and the intermediate-scale longitudinal structures. Coherent structures contribute about 36% to the spanwise vorticity variance at x/d=10. The contribution decreases rapidly to about 5% at x/d=40. The present results suggest that vorticity largely reside in relatively small-scale structures.     

Vortex formation process in gravity-driven starting jets

Experiments were conducted to investigate the vortex formation process in gravity-driven starting jets at three Reynolds numbers ( 3,528 and 4,716, where D is the nozzle diameter, the average discharging velocity and ν the kinematic viscosity). Planar laser induced fluorescence (PLIF) was used to visualize the flow while particle image velocimetry (PIV) was used to quantify the 2D velocity and vorticity fields. Vortex leapfrogging was found to appear at Re = 2,358, while pinch-off for the leading vortex was observed in the other two cases. Conditions for the vortex pinch-off appeared to be different from those found in the piston-driven starting jets. Although the leading vortex attained the maximum circulation and energy level very quickly after the jet was initiated, its detachment from the main stem may not be achieved necessarily. Attempt had been made to elaborate the different flow characteristics for respective cases from the related distribution of circulation and energy level via the measured velocity field by PIV.     

Drag reduction and thrust generation by tangential surface motion in flow past a cylinder

Sensitivity of drag to tangential surface motion is calculated in flow past a circular cylinder in both two- and three-dimensional conditions at Reynolds number \(\textit{Re} \le 1000\). The magnitude of the sensitivity maximises in the region slightly upstream of the separation points where the contour lines of spanwise vorticity are normal to the cylinder surface. A control to reduce drag can be obtained by (negatively) scaling the sensitivity. The high correlation of sensitivities of controlled and uncontrolled flow indicates that the scaled sensitivity is a good approximation of the nonlinear optimal control. It is validated through direct numerical simulations that the linear range of the steady control is much higher than the unsteady control, which synchronises the vortex shedding and induces lock-in effects. The steady control injects angular momentum into the separating boundary layer, stabilises the flow and increases the base pressure significantly. At \(\textit{Re}=100\), when the maximum tangential motion reaches 50% of the free-stream velocity, the vortex shedding, boundary-layer separation and recirculation bubbles are eliminated and 32% of the drag is reduced. When the maximum tangential motion reaches 2.5 times of the free-stream velocity, thrust is generated and the power savings ratio, defined as the ratio of the reduced drag power to the control input power, reaches 19.6. The mechanism of drag reduction is attributed to the change of the radial gradient of spanwise vorticity (\(\partial _{r} \hat{\zeta }\)) and the subsequent accelerated pressure recovery from the uncontrolled separation points to the rear stagnation point.     

Semianalytical Solution for \text{ CO}_{2} Plume Shape and Pressure Evolution During \text{ CO}_{2} Injection in Deep Saline Formations

The injection of supercritical carbon dioxide ( $\text{ CO}_{2})$ in deep saline aquifers leads to the formation of a $\text{ CO}_{2}$ rich phase plume that tends to float over the resident brine. As pressure builds up, $\text{ CO}_{2}$ density will increase because of its high compressibility. Current analytical solutions do not account for $\text{ CO}_{2}$ compressibility and consider a volumetric injection rate that is uniformly distributed along the whole thickness of the aquifer, which is unrealistic. Furthermore, the slope of the $\text{ CO}_{2}$ pressure with respect to the logarithm of distance obtained from these solutions differs from that of numerical solutions. We develop a semianalytical solution for the $\text{ CO}_{2}$ plume geometry and fluid pressure evolution, accounting for $\text{ CO}_{2}$ compressibility and buoyancy effects in the injection well, so $\text{ CO}_{2}$ is not uniformly injected along the aquifer thickness. We formulate the problem in terms of a $\text{ CO}_{2}$ potential that facilitates solution in horizontal layers, with which we discretize the aquifer. Capillary pressure is considered at the interface between the $\text{ CO}_{2}$ rich phase and the aqueous phase. When a prescribed $\text{ CO}_{2}$ mass flow rate is injected, $\text{ CO}_{2}$ advances initially through the top portion of the aquifer. As $\text{ CO}_{2}$ is being injected, the $\text{ CO}_{2}$ plume advances not only laterally, but also vertically downwards. However, the $\text{ CO}_{2}$ plume does not necessarily occupy the whole thickness of the aquifer. We found that even in the cases in which the $\text{ CO}_{2}$ plume reaches the bottom of the aquifer, most of the injected $\text{ CO}_{2}$ enters the aquifer through the layers at the top. Both $\text{ CO}_{2}$ plume position and fluid pressure compare well with numerical simulations. This solution permits quick evaluations of the $\text{ CO}_{2}$ plume position and fluid pressure distribution when injecting supercritical $\text{ CO}_{2}$ in a deep saline aquifer.     

CH4/O2/CO2预混体系爆炸动力学研究

为探究甲烷在富氧条件下的火焰动力学规律,以CH₄/O₂/CO₂预混体系为研究对象,在小尺度方形透明管道中进行了一系列爆炸实验,探讨了初始环境温度波动对爆炸参数的影响,并对预混体系的燃烧机理进行分析。结果表明:在273 K的环境温度下,化学当量比φ=0.8~1.0且氧气相对比γ<0.30和φ=1.2且γ<0.35的预混体系不能被点燃,而其他预混体系均可被点燃,最终产生郁金香与非郁金香两种火焰类型,并且根据郁金香火焰独特的演变特征,又划分为T形郁金香火焰和不对称郁金香火焰;随着γ的增大,无量纲火焰传播速度v/(S_Lσ)的变化趋势由“两升两降”转变为“一升一降”。初始环境温度的升高并未对火焰传播速度和爆炸超压的变化趋势产生影响,但是会导致最大爆炸超压p_max和最大火焰传播速度降低。值得注意的是,初始环境温度对爆炸强度的影响随化学当量比的减小而增强。另外,与最大爆炸超压相比,最大火焰传播速度与层流燃烧速度之间的关系更紧密。从敏感性分析中可知:层流燃烧速度对自由基链式反应R38(即H+O₂=O+OH)表现出最大的正敏感度,对R52(即H+CH₃(+M)=CH₄(+M))表现出最大的负敏感度,并且对自由基OH的生成速率最敏感,当初始环境温度升高至303 K时,层流燃烧速度对R38(正)和R52(负)的敏感度降低;H、O和OH自由基总摩尔分数的增大会削弱热扩散的不稳定性,增强流体力学的不稳定性。     

水中高压脉动气泡与浮体流固耦合特性研究

本文针对水中放电气泡与水面浮体流固耦合作用开展实验和数值研究, 采用边界积分法对气泡运动进行数值模拟, 利用辅助函数法提高非线性流固耦合问题的计算精度, 同时运用双节点法保证气-液-固三相交界线的计算稳定性. 实验中, 采用水下放电技术生成气泡, 使用高速摄影捕捉气泡动力学行为与浮体运动响应. 首先对比数值与实验结果, 二者吻合良好, 验证了数值计算模型的有效性和正确性. 然后通过对气泡与浮体的无量纲距离$\gamma_{s} $ (气泡最大半径为特征长度)进行系统研究发现: (1) $\gamma_{s} $从0.2增大至2时, 气泡在坍塌阶段分别形成了颈缩型环状射流(本文针对水中放电气泡与水面浮体流固耦合作用开展实验和数值研究,采用边界积分法对气泡运动进行数值模拟,利用辅助函数法提高非线性流固耦合问题的计算精度,同时运用双节点法保证气-液-固三相交界线的计算稳定性.实验中,采用水下放电技术生成气泡,使用高速摄影捕捉气泡动力学行为与浮体运动响应.首先对比数值与实验结果,二者吻合良好,验证了数值计算模型的有效性和正确性.然后通过对气泡与浮体的无量纲距离γ_s(气泡最大半径为特征长度)进行系统研究发现:(1)γ_s从0.2增大至2时,气泡在坍塌阶段分别形成了颈缩型环状射流(0.2≤γ_s≤0.3)、接触射流(0.4≤γ_s≤0.6)、非接触射流(0.7≤γ_s≤1)、对射流(1.1≤γ_s≤1.3)和反射流(1.4≤γ_s≤2)等5种典型射流模式;(2)正射流速度随γ_s先增大后减小再增大,并且当0.7≤γ_s≤0.9时,速度可达约1000 m/s;反射流速度随γ_s增大而增大;(3)在本文实验条件下,γ_s1.5时浮体对气泡的Bjerknes吸引力强于自由液面的Bjerknes排斥力导致气泡在坍塌阶段向浮体迁移;当γ_s≥1.5时自由液面对气泡的排斥作用更强,气泡在坍塌阶段远离自由液面.