GENERAL REYNOLDS ANALOGY RELATION ON BLUNT-NOSED BODIES$^{\bf 1)}$

钝头体壁面的摩阻和热流分布规律不同,平板流动中的雷诺比拟关系在钝头体壁面失效. 文章在前期高超声速广义雷诺比拟理论研究工作的基础上,利用数值仿真的方法对不同外形和来流参数条件下的钝头体广义雷诺比拟关系开展进一步研究. 通过建立钝头体绕流边界层的理论分析模型,得到了钝头体壁面雷诺比拟系数的线性分布预示公式. 采用数值求解 N-S 方程的方法,计算了圆柱和幂次体壁面的摩阻和热流以及二者之间的比拟系数. 通过与前期数值和理论结果对比,以及计算收敛性和网格无关性检验,对数值方法进行了验证. 通过在不同雷诺数 ($Re_\infty = 3.98\times 10^2 \sim 1.59\times 10^6$) 和马赫数 ($M_\infty = 3\sim 12$) 条件下的计算结果对比分析雷诺比拟系数的分布,总结了钝头体中广义雷诺比拟关系受外形和来流条件的影响,评估了广义雷诺比拟理论的适用性. 研究发现,在较高雷诺数条件下,离驻点较远的下游 ($\theta > 60^\circ$) 部位,雷诺比拟系数的分布不同程度地偏离理论预示的线性规律. 相比于圆柱外形,幂次体壁面的雷诺比拟系数分布的线性规律相对较好,其分布斜率略低于圆柱壁面的结果. 研究表明,如果针对实际外形和雷诺数进行适当修正,可以提高广义雷诺比拟关系的预示精度.     

等速上仰翼型动态失速现象研究

翼型大迎角绕流的静态失速将造成升力突降和气动性能急剧恶化，但利用非定常运动所产生 的动态失速效应，可以大大地延缓气流分离和失速现象的发生. 采用Rogers发 展的双时间步Roe格式，求解拟压缩性修正不可压N-S方程. 数值模拟了低雷诺数 ($Re=4.8 \times 10^{4}$)条件下NACA0015翼型作等速上仰($\alpha =0^{\circ} \sim 60^{\circ}$)的动态失速过程，同Walker的试验结果比 较，验证了计算结果的正确性. 研究了该过程中主涡、二次涡和三次涡的发展，升 力系数随攻角变化，以及不同上仰速度对动态失速效应所造成的影响.     

研究了雷诺数Re=200, 1000, 线速度比$\alpha =0.5$,

研究了雷诺数Re=200, 1000, 线速度比$\alpha =0.5$, 2.0, 4.0, 强迫振荡频率$f_{s}=0.1\sim 2.0$情况下的旋转振荡圆柱绕流问题. 通 过基于非结构同位网格有限体积法对Navier-Stokes方程进行数值求解. 对流项、扩 散项和非恒定项的离散格式均具有二阶精度，利用SIMPLE算法处理压力-速度耦合. 计算得到了作用力系数随不同控制参数的变化规律. 通过对升力系数的频谱分析得到 自然脱落频率和强迫振荡频率下的作用力振幅. 通过对不同频率作用力幅值的分析， 得到频率之间的竞争关系，进而定量地给出了不同尾迹涡脱落模式的分区图.     

低雷诺数沟槽表面湍流/非湍流界面特性的实验研究

湍流/非湍流界面是流动中湍流和无旋流的边界,其相关研究在加深对湍流与无旋流之间的物质、动量和能量交换的理解有重要意义.本文采用时间解析的二维粒子图像测速技术,分别对零压梯度光滑、顺流向锯齿形沟槽表面平板在不同雷诺数下对湍流/非湍流界面的几何特征及动力学特性进行了实验研究.实验雷诺数为$Re_{\tau } =400\sim1000$.本文采用了湍动能准则对湍流/非湍流界面进行了识别,并分析界面高度分布、分形特征及界面附近的条件平均速度和涡量.结果表明在不同雷诺数下, 无论是光滑壁面还是沟槽壁面,界面平均高度在0.8 $\sim$ 0.9$\delta_{99} $附近. 对于沟槽壁面而言,减阻时对应的界面高度的概率密度分布与光滑壁面基本一致, 均遵循正态分布,而当阻力增大时, 界面高度分布偏离正态分布出现正的偏度. 在本实验情况下,界面分形维度、跨界面速度跳变均会随着雷诺数增大而增大. 此外,不同壁面情况下无量纲条件平均涡量在界面附近的分布相近,而界面附近无量纲速度梯度最大值近似为常数.      

EXPERIMENTAL STUDY ON PROPERTIES OF TURBULENT/NON-TURBULENT INTERFACE OVER RIBLETS SURFACES AT LOW REYNOLDS NUMBERS1)

湍流/非湍流界面是流动中湍流和无旋流的边界,其相关研究在加深对湍流与无旋流之间的物质、动量和能量交换的理解有重要意义.本文采用时间解析的二维粒子图像测速技术,分别对零压梯度光滑、顺流向锯齿形沟槽表面平板在不同雷诺数下对湍流/非湍流界面的几何特征及动力学特性进行了实验研究.实验雷诺数为$Re_{\tau } =400\sim1000$.本文采用了湍动能准则对湍流/非湍流界面进行了识别,并分析界面高度分布、分形特征及界面附近的条件平均速度和涡量.结果表明在不同雷诺数下, 无论是光滑壁面还是沟槽壁面,界面平均高度在0.8 $\sim$ 0.9$\delta_{99} $附近. 对于沟槽壁面而言,减阻时对应的界面高度的概率密度分布与光滑壁面基本一致, 均遵循正态分布,而当阻力增大时, 界面高度分布偏离正态分布出现正的偏度. 在本实验情况下,界面分形维度、跨界面速度跳变均会随着雷诺数增大而增大. 此外,不同壁面情况下无量纲条件平均涡量在界面附近的分布相近,而界面附近无量纲速度梯度最大值近似为常数.     

垂直壁面附近上升单气泡的弹跳动力学研究

采用高速摄影技术结合阴影法,对静止水中垂直壁面附近上升单气泡运动进行实验研究,对比气泡尺度及气泡喷嘴与壁面之间的初始无量纲距离 ($S^{\ast}$)对气泡上升运动特性的影响,分析气泡与壁面碰撞前后,壁面效应与气泡动力学机制及能量变化规律.结果表明,对于雷诺数$Re \approx 580 \sim 1100$,无量纲距离$S^{\ast } <2 \sim3$时,气泡与壁面碰撞且气泡轨迹由无约束条件下的三维螺旋转变成二维之字形周期运动;当$S^{\ast } >2 \sim3$时,壁面效应减弱,有壁面约束的气泡运动与无约束气泡运动特性趋于一致.气泡与壁面碰撞前后,壁面效应导致横向速度峰值下降为原峰值的70%,垂直速度下降50%;气泡与壁面碰撞前,通过气泡中心与壁面距离($x/R$)和修正的斯托克斯数相关式可预测垂直速度的变化规律.上升气泡与壁面碰撞过程中,气泡表面变形能量单向传输给气泡横向动能,使得可变形气泡能够保持相对恒定的弹跳运动.提出了气泡在与壁面反复弹跳时的平均阻力系数的预测模型,能够很好地描述实验数据反映出的对雷诺数${Re}$、韦伯数${We}$和奥特沃斯数${Eo}$等各无量纲参数的标度规律.      

不同壁面取向下超疏水平面直轨道上的气泡滑移

利用特定几何分布的超疏水表面实现气泡定向输运在矿物浮选和生物孵化等领域具有广阔的应用前景, 对平面直线超疏水轨道而言, 其壁面取向是相关工程结构的关键参数, 但超疏水壁面取向对倾斜壁面气泡滑移的影响尚不明确. 本文采用高速阴影成像系统研究了不同壁面取向($-90^\circ\leqslant \beta \leqslant 90^\circ$)及轨道倾角($45^\circ\leqslant \alpha \leqslant 75^\circ$)下, 气泡($D_{eq}=2.4$ mm, $Re=500$ $\sim$ 700, $We=7$ $\sim$ 13)在轨道宽度为2 mm的超疏水直线轨道上的运动特性. 气泡在轨道上的滑移近似为匀速, 形状为具有多脊的半子弹型. 根据气液界面波动程度的不同, 滑移气泡可分为波动型和稳定型, 稳定型气泡只在较小倾角且较大方位角时出现($45^\circ\leqslant \alpha < 70^\circ$, $| \beta | \geqslant 45^\circ$). 根据倾角不同, 滑移速度关于$\beta $有2种变化规律: 当$\alpha \leqslant 65^\circ$, 气泡滑移速度近似为关于$\beta =0^\circ$ 的单峰分布($\beta =0^\circ$时, 气泡滑移速度最大); 当$\alpha \geqslant 70^\circ$, 气泡滑移速度在不同的方位角下基本保持稳定. 气泡的最大滑移速度可达0.66 m/s ($\beta =0^\circ$, $\alpha =70^\circ$), 远大于相同尺度的自由上升气泡($\approx0.25$ m/s), 这主要是壁面浸润性分布和惯性力的耦合效应所致. 轨道取向(方位角$\beta )$及轨道倾角($\alpha )$通过改变气泡沿轨道方向的驱动力和气泡迎风面积影响气泡的滑移速度和气液界面稳定性.      

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

高温后高强高性能混凝土双轴压力学性能

利用大型静动真三轴试验机,进行了常温20${^\circ}$C以及200${^\circ}$C$\sim $ 600${^\circ}$C\,6个温度等级高温后高强高性能混凝土在7种应力比双轴压应力状态下的强度与变形试验.测得了双轴主压方向的静态强度、峰值应变与应力应变曲线,剖析了温度和应力比对单、双轴压强度与峰值应变发展趋势的影响规律性以及试件破坏形态. 试验结果表明:随温度的升高,高强高性能混凝土的单轴压减摩强度并不一定降低;双轴压强度相对于单轴压强度的提高倍数取决于应力比、不同温度等级后的高强高性能混凝土``脆硬性'. 提出了带有温度和应力比参数的Kupfer-Gerstle破坏准则公式.      

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

两不等直径圆柱高雷诺数时均风力场的数值模拟

圆柱间气动干扰研究具有重要的理论和现实意义.尽管国内外开展了圆柱组合风效应的风洞试验研究,但主要针对等直径圆柱,并且雷诺数多为10~5以下.考虑到工程结构风场的高雷诺数特征,采用数值模拟方法,模拟两不等直径圆柱在串列、并列及交错排列下的高雷诺数(Re=4.5×10~5)时均绕流场.通过改变组合的间距和风向,分析探讨两柱阻力、升力及总风力的变化规律.     

Newtonian variational problems for nonshallow bodies

On the basis of the analysis of variational problems in which the Newton formula is used to calculate the pressure, it is shown that: 1) the laws of variation with length of the areas of the cross sections of axisymmetric bodies and bodies whose cross sections have the form of regular polygons are identical; 2) the law of variation of the radius of the body of minimum drag along its length depends on the magnitude of the specified base radius or body volume. In spite of the indeterminacy of the applicability conditions, the Newton law for the pressure on the surface of a body moving at hypersonic speed

$$p = \rho _\infty V_\infty ^2 , \sin ^2 \alpha $$     

压电体椭圆孔边的力学分析

基于复变函数的方法,以PZT-4材料为例,分别采用精确电边界条件和非导通电边界条件进行了远场均匀载荷作用下的横观各向同性压电体椭圆孔的力学分析并与相关结果进行对比。结果表明当椭圆孔退化为圆孔时,无论在远场作用力载荷或电载荷,两种电边界条件下的结果均能完全吻合。随着椭圆孔的愈加尖锐化,非导通电边界条件逐渐不能适用。     

Asymptotic solution of the two-dimensional equations for a thin viscous rarefied shock layer on a cold surface

Hypersonic rarefied flow past blunt bodies is studied in the continuum-free-molecular transition regime. On the basis of an asymptotic analysis three rarefied gas flow patterns are established depending on the relation between the relevant parameters of the problem. In the first regime corresponding to a cold surface asymptotic solutions of the equations of a thin viscous shock layer are derived at low Reynolds numbers in the axisymmetric and plane cases. Simple analytical expressions for the pressure and the heat transfer and friction coefficients are obtained as functions of the freestream parameters and the body geometry. With decrease in the Reynolds number the coefficients approach the values corresponding to free-molecular flow. In this regime a similarity parameter for the hypersonic rarefied flow past bodies is determined. The asymptotic solutions are compared with numerical solutions and the results of direct statistical simulation by the Monte Carlo method.     

Active control of flow past an elliptic cylinder using an artificial neural network trained by deep reinforcement learning

The active control of flow past an elliptical cylinder using the deep reinforcement learning (DRL) method is conducted. The axis ratio of the elliptical cylinder $\Gamma$ varies from 1.2 to 2.0, and four angles of attack $\alpha=0^\circ, 15^\circ, 30^\circ$, and $45^\circ$ are taken into consideration for a fixed Reynolds number $Re=100$. The mass flow rates of two synthetic jets imposed on different positions of the cylinder $\theta_1$ and $\theta_2$ are trained to control the flow. The optimal jet placement that achieves the highest drag reduction is determined for each case. For a low axis ratio ellipse, i.e., $\Gamma=1.2$, the controlled results at $\alpha=0^\circ$ are similar to those for a circular cylinder with control jets applied at $\theta_1=90^\circ$ and $\theta_2=270^\circ$. It is found that either applying the jets asymmetrically or increasing the angle of attack can achieve a higher drag reduction rate, which, however, is accompanied by increased fluctuation. The control jets elongate the vortex shedding, and reduce the pressure drop. Meanwhile, the flow topology is modified at a high angle of attack. For an ellipse with a relatively higher axis ratio, i.e., $\Gamma\ge1.6$, the drag reduction is achieved for all the angles of attack studied. The larger the angle of attack is, the higher the drag reduction ratio is. The increased fluctuation in the drag coefficient under control is encountered, regardless of the position of the control jets. The control jets modify the flow topology by inducing an external vortex near the wall, causing the drag reduction. The results suggest that the DRL can learn an active control strategy for the present configuration.     

Asymptotic solutions of the two-dimensional equations for a thin viscous shock layer in a rarefied gas

Two-dimensional hypersonic rarefied gas flow around blunt bodies is investigated for the continuum to free-molecular transition regime. In [1], as a result of an asymptotic analysis, three rarefied gas flow regimes, depending on the relationship between the problem parameters, were detected and one of these regimes was investigated. In the present study, asymptotic solutions of the thin viscous shock layer equations at small Reynolds numbers are obtained for the other two flow regimes. Analytical expressions for the heat transfer, friction and pressure coefficients are obtained as functions of the incident flow parameters and the body geometry and temperature. As the Reynolds number tends to zero, the values of these coefficients approach their values in free-molecular flow. The scaling parameters of hypersonic rarefied gas flow around bodies are determined for different regimes. The asymptotic solutions are compared with the results of direct Monte Carlo simulation.     

Theoretical amplitude and period of precursor solition generation in two-layer flow

An fKdV equation of two-layer flow and an averaged fKdV equation (AfKdV equation) with respect to phase are derived to determine the theoretical amplitude and period of the precursor solitons in the present paper. In terms of the AfKdV equation derived by the authors, a new theory on the precursor soliton generation based on Lee et al.'s concept is presented. Concepts of asymptotic mean hydraulic fall and level are introduced in our analysis, and the theoretical amplitude and period both depend on the asymptotic mean levels and stratified parameters. From the present theoretical results, it is obtained that when the moving velocity of the topography is at the resonant points, there exist two general relations: (1) amplitude relation Å=2F, (2) period relation $\mathop \tau \limits^ \circ = - 8m_1 m_3^{ - 1} \sqrt {6m_4 m_3^{ - 1} } \mathcal{F}$ , in which Å and $\tau $ are the amplitude and period of the precursor solitons at the resonant points respectively,m ₁,m ₃ andm ₄ are coefficients of the fKdV equation, andF is an asymptotic mean half-hydraulic fall at subcritical cutoff points. The theoretical results of this paper are compared with experiments and numerical calculations of two-layer flow over a semicircular topography and all these results are in good agreement. Due to the canonical character of the coefficients of fKdV equations, this theory also holds for any two-dimensional system, which can be reduced to fKdV equations.     

Investigation on the mechanism of convective mass transfer from a horizontal rotating cylinder without air jet flow

The local and mean convective mass transfer coefficients from the surface of a large-diameter horizontal circular rotating cylinder without air jet flow were investigated by measuring the concentration gradient. The results indicate that rotation performs different effects on the convective mass transfer at different regions. Based on the experimental data, the correlation equations of the mean convective mass transfer Sherwood number Sh and the critical Reynolds number Re _r,cri have been formulated as follows: $ Sh = 0.32[(8.5Re_{\text{r}}^{2} + Gr) \cdot Sc]^{1/3} $ and $ Re_{\text{r,cri}} = 0.44(Gr \cdot Sc)^{1/2} $ .